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PREFACE 


The FIBER-TEX 1991 proceedings contain the papers jointly sponsored by the National Aeronau- 
tics and Space Administration, the Center for Advanced Engineering Fibers at Clemson University, 
the Department of Defense, the North Carolina State University, and Drexel University. The con- 
ference was held in Raleigh, North Carolina on October 15-17, 1991 to create a forum to encourage 
an interrelationship of the various disciplines involved in the fabrication of materials, the types of 
equipment, and the processes used in the production of advanced composite structures. Topics 
discussed were Reinforcing Fibers, Matrix Materials, Mechanics of Woven Materials, Structural 
Fabric Production, Pultrusion and Composite Processes, and Structures and Applications. 

Certain materials and processes are identified in this publication in order to adequately specify 
procedures. In no case does such identification imply recommendation or endorsement by the 
government, nor does it imply that the materials or processes are the only or best ones available 
for the purpose. 


John D. Buckley 

NASA Langley Research Center 
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An Evaluation of Composites 
Fabricated from Powder Epoxy 

Towprea 


J. Timothy Hartness and Tim Greene 
BASF Structural Materials 


Abstract 

BASF has developed a unique process for applying powdered resin systems to 
continuous reinforcement fibers in order to produce flexible towpreg 
material. Evaluation of three powder epoxy resins by BASF using this 
towpregging process is in progress under NASA contract NASI-18834. Shell 
RSS-1952, Dow CET-3, and 3M PR500 powder epoxy systems have been 
successfully towpregged with G30-500 6K carbon fiber. Both neat resin and 
basic unidirectional composite properties have been developed to compare 
performance. Cure cycles for each system have also been developed for 
repeatable fabrication of high-quality composite laminates. Evaluations of 
the powder towpreg material for use in textiles processes such as weaving 
and braiding are underway. Traditional 8-harness weaving has been 
successfully performed with one system (PR500/G30-500) to date, with 
some basic composite properties generated. Ongoing work will demonstrate 
scaleup of the towpregging process for higher throughput, as well as 
evaluation of the powder towpreg material in advanced preforming processes 
such as 3-D braiding and weaving. 
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POWDER EPOXY TOWPREG DEVELOPMENT 

• NASA Contract/Materials Characterization 

• Primary Focus Subsonic Applications (180°F Service) 

• Evaluate Mechanical Performance 

• Evaluate Use in Textiles Preforming 

• Evaluate Manufacturing/Processing Methods 

ADVANTAGES OF "DRY" MATERIAL FORMS 

• Conformability 

• Textile Yarn Form/Textile Technology 
Applications 

• No Refrigeration Required 
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POWDER EPOXY TOW ADVANTAGES 


. Predetermined Fiber Volume 
"Predetermined" Chemistry 
Good Fiber/Resin Distribution 
. No Solvents 

• Room Temperature Storage 


CANDIDATE POWDER EPOXY SYSTEMS 


. PR-500 (3M) 

• RSS-1952 (SHELL CHEMICAL) 


CET-3 (DOW CHEMICAL) 



PHYSICAL PROPERTIES. NEAT RESIN 



EB-500 (3M) 

HSS-1952 (SHELL) 

CET-3P (DOW) 

Tg (DSC, °C/°F) 

205/401 

219/426 

1 64/327 

TENSILE STRENGTH (KSI) 

8.3 


13.0 

MODULUS (KSI) 

507 


4 1 0 

ELONGATION (%) 

1.9 


5.0 

FLEXURAL STRENGTH (KSI) 

18.4 

16.9 

21.0 

MODULUS (KSI) 

504 

426 

45 0 

STRAIN (%) 

4.2 

5.1 

7.0 

DENSITY (gm/cc) 

1 .25 

1.15 

1.27 

MOISTURE ABSORPTION (% WT.) 1.56 

1.1 

1 .35 


B47B3-51 177 DEG C CURING POWDER 



PR-500 Powder, Viscosity Sweep 
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Powder Coated Tow, 1000X 
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CURE CYCLE. FUSED TOWPREG 


PR-500 (3M): 350°F/2 hours 

RSS-1952 (Shell): 300°F/2 hours, ramp to 400°F/4 hours 

CET-3P (Dow): 300°F/4 hours; 400°F/4 hour post-cure 


COMPOSITE PROPERTIES. UNIDIRECTIONAL TOW 


FIBER VOLUME 


PR-500/G30-500 RSS- 1 952/G30-50n 
55% 63% 


0° 3 PT. FLEXURE (RT. 32-11 
STRENGTH (KSI) 

MODULUS (MSI) 

0° 4 PT. SHEAR (ffOM) 
STRENGTH (KSI) 


90° 3 PT. FLEXURE 1RT1 
STRENGTH (KSI) 
MODULUS (MSI) 


242 
1 6 


12.2 


1 1.0 
1.2 


320 
1 9 


10.0 


9.0 

1.18 


PHYSICAL PROPERTIES. 8-HARNESS FABRIC 


RESIN: PR-500 (3M) 

FIBER VOLUME: 56% 

3 PT. FLEXURE (RT. 32-1.1 

STRENGTH (KSI) 102.0 

MODULUS (MSI) 7.1 

4 PT. SHEAR (RT. 16-1) 

STRENGTH (KSI) 6.0 



Photomicrograph, Cured 8 Harness Laminate 


CONCLUSIONS/ACCOMPLISHMENTS 


• Fused Towpreg Approach Is Viable With Powder Epoxy 

• 5-Ends Successfully Demonstrated 

• Fused Epoxy Tow Is Weaveable Using Standard Techniques 

• Good Fiber/Resin Distribution and Wet-Out Demonstrated 

• Initial Unidirectional and 8-Harness Data Generated 

FUTURE PLANS 

• Generate RSS-1952 and CET-3 8-Harness Fabric Data 

• Complete RT and Hot/Wet Mechanical Testing 

• Determine Processing Window 

• Verify Preforming Feasibility (2-D and 3-D) 

• Develop De-bulklng and Part-Manufacturing Methods 
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Frank L. Scardino 

I FAI Professor of Industrial Fabrics 
Philadelphia College of Textiles and Science 
School House Lane and Henry Avenue 
Phi ladelphia , Pennsylvania 

INTRODUCTION 

In the design of textile composites, the selection of materials 
and constructional techniques must be matched with product performance, 
productivity, and cost requirements. Constructional techniques may 
vary from slow and expensive, hand laid-up batch manufacturing (one 
unit at a time) to very quick and cost effective continuous pull- 
through processing. No single textile material, resin system, textile 
assemblage, or constructional technique can be considered optimal for 
all applications. 

A classification of various textile composite systems is given in 
Table I. In general, the chopped fiber system (Type I) is not suitable 
for structural composite applications because of fiber discontinuity, 
uncontrolled fiber orientation and a lack of fiber integration or 
entanglement. Linear filament yarn systems (Type II) are quite 
acceptable for structural components which are exposed to simple 
tension in their applications. To qualify for more general use as 
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structural components, filament yarn systems must be multi- 
directionally positioned. With the most sophisticated filament winding 
and laying techniques, however, the Type II systems have limited 
potential for general load-bearing applications because of a lack of 
filament integration or entanglement, which means vulnerability to 
splitting and delamination among filament layers. 


TABLE I. TEXTILE COMPOSITE SYSTEMS 



Re i nf orcemen t Textile 

Fiber 

Fiber 

Fiber 

IXE1 

System 

Construction 

Leny th 

Orientation 

En tang lement 

l 

Suspended 

Chopped Fiber 

Di scont i nuous 

Uncontrol led 

None 

II 

Linear 

Filament Yarn 

Continuous 

Linear 

None 

[ii 

Laminar 

Simple Fabric 

Continuous 

Planar 

Planar 

IV 

I n teg rated 

Advanced Fabric 

Continuous 

3-D 

3-D 


The laminar systems (Type III) represented by a variety of simple 
fabrics (woven, knitted, braided and nonwoven) are especially suitable 
for load-bearing panels in flat form and for beams in a rolled up or 
wound form. The main features of simple fabric systems are fiber 
continuity, planar fiber orientation and planar fiber entanglement or 
integration, in general. The major vulnerability of simple fabric 
laminate systems is delamination between layers of the fabrics which 
tends to be more critical in flat panels than in rolled up tubular or 
rectangular configurations. 

The totally integrated, advanced fabric systems (Type IV) are 
thought to be the most reliable for general load-bearing applications 
because of fiber continuity and because of controlled multiaxial fiber 


orientation and entanglement. Consequently, the risk of splitting and 
delamination is minimized and practically omitted. Type IV systems can 
be woven, knitted, braided or stitched through with very special 
equipment. 

In general, multiaxial fabrics are classified as Type IV in Table 
I. A practical advantage of multiaxial fabrics is the elimination of 
much of the hand lay-up work in composite manufacturing which is so 
labor intensive and time consuming. Also, multiaxial fabrics are 
easier to handle because the various yarn orientations are held in a 
fixed position during manipulation. (1) 

MULTIAXIAL FABRIC TECHNOLOG IBS 

Several alternate technologies are commercially available today 
for the conversion of yarn into multiaxial constructions for a variety 
of industrial fabric applications, but particularly for flexible and 
rigid composites. These multiaxial fabric technologies include 
adhesive bonding, triaxial weaving, triaxial braiding, weft insertion 
warp knitting, bias ply stitch bonding and bias warp knitting. Each 
multiaxial fabric technology has unique attributes and limitations, and 
accordingly, will find a place in the industrial fabric market on a 
cost/ pe rf ormance/ava ilability/processabili ty/machineabi 1 i ty/ join- 
ability/maintainability basis. 

Adhes ively-bonded multiaxial systems are easily produced by 
combining several layers of yarn or fabric at various angles (including 
skewed WIWK fabrics, bias woven, etc.). While the productivity is 
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attractive, adhesive bonding does not provide resistance to 
delamination or to crack propagation. Another major limitation of 
adhesively bonded yarn and fabric systems is a lack of through-the- 
thickness strength. Adhesively-bonded constructions are classified as 
Type II if made directly from yarn or Type ill if made by combining 
layers of fabric. By contrast, all of the following multiaxial fabric 
systems are classified as Type IV. 

In triaxial weaving technology, three systems of yarns are 
interlaced at sixty degree angles (2,3). With two warps and one 
filling, the resulting orientation is + 30 O /-30 O /90° . While triaxial 
weaving technology has been developed for quite some time, triaxial 
woven fabric is commercially available in limited styles from Sakase 
Textile Company in Japan. 

Triaxial braids incorporate one longitudinal yarn system with two 
helical or diagonal yarn systems. Consequently, the resulting yarn 
orientation can vary from 0 O /+30°/-30° to 0°/+60°/-60° , depending upon 
processing variables. Triaxial braids are particularly well suited for 
tubular applications .( 4 ) 

Weft insertion warp k nit (WIWK) technology is well developed and 
established Cor a variety of industrial fabric applications. Biaxial 
WIWK fabrics are made by knitting 0° and 90° straight- laid yarn systems 
together. It is possible to produce quasi-quadraxial WIWK fabrics by 
adding two additional laid-in yarn systems in a counter zig-zag fashion 
in the machine direction, wherein each yarn continues for up to two 
inches in the +45° or -45° direction before reversing its path. (5, 6) 


Several companies are capable of producing quas i-quadrax ial WIWK 
fabrics on LIBA or Mayer WIWK machines. 

Multiaxial stitch bonded fabric preforms can be made by stitching 
together several layers of yarn at various angles or plies of skewed 
fabric. The stitching processes provides through- the- thickness 
strength and integrity. Stitching can be performed efficiently and 
with high productivity on any suitable multi-needle process capable of 
through-the-thickness penetration including sewing machines, Mali 
machines, warp knitting stitch- through machines, etc. The major 
advantages of multiaxial st itched— through fabric preforms are high 
density, control of yarn orientation in each layer, and the integration 
of the yarn/fabric layers by stitching. The major disadvantage is the 
impaling and localized dislocation of yarn during stitching which leads 
to reduced strength and to poor structural cons is tency . ( 7 ) 

Several companies are presently involved in multiaxial stitch- 
through technology including Hexcel in Seguin, Texas; Bean Fiberglass 
Company in Jaffrey, New Hampshire; Brunswick Technologies in Brunswick, 
Maine; Bay Mills in Ontario, Canada; Chomarat in France; and others in 
Europe; each of whom lias devised their own unique multiaxial orienta- 
tion system. Commercial stitch-through machines without bias 
orientation capability are available from LIBA, Karl Mayer and Textima. 
Also, LIBA is presently marketing a machine called Copcentra-Mult iaxial 
which coordinates multiaxial orientation with stitching through on a 
weft insertion warp knit type of machine. Multiaxial stitch-through 
fabrics are available in several plies of yarn systems oriented in the 
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warp, weft and bias directions. Light nonwoven webs can be added, as a 
P-iy, it desired, and fabric widths of tour to eight feet are possible. 

All of the laid- i n yarns are linear, are in a specific plane, and 
are continuous between fabric edges. Obviously, multiaxial stitch- 
through fabrics can be engineered for specific directional properties. 
Currently, most fabrics are composed of high modulus laid- in tow/yarn 
systems with a fine polyester stitching yarn system. Fabric con- 
structions of the various ply orientations can be made up to 
approximately 1/4— inch thickness. Multiaxial stitch— through fabric 
weights vary according to construction but are often in the range 
between 12 and 48 ounces per square yard. 

Mult iaxial warp k ni t fabric technology has been developed recently 
by the Karl Mayer Textile Machine Company in West Germany. The process 
is officially known as multiaxial magazine weft insertion and 
informally referred to as the bias machine. As an alternate to 


multiaxial stitch-through technology, the Mayer bias machine precisely 
knits rather than punches the stitching yarn through the various layers 
of laid-in yarns. Consequently, no yarns are impaled and the texture 
is quite uniform throughout the fabric, leading to higher translation 
efficiencies. The major design limitation of the Mayer bias machine is 
a maximum of four yarn layers (0°, 90°, +45°, -45°) plus a fiber web, 
if desired. Also, fabrics made from the Mayer bias machine tend to be 
somewhat more voluminous than comparable stitch-through multiaxial 
fabrics, which provides for easier resin penetration but slightly lower 
fiber volume fractions, potentially. ( 6 ) Milliken is the only company 




in the U.S. with a Mayer bias machine. (8) Nine other Mayer bias 
machines are in operation in Europe. 

MULT I AXIAL FABRIC BEHAVIOR 

The behavior of various fabric constructions under uniaxial stress 
in the machine (0°), crosswise (90°) and bias 45°) directions are 
illustrated in Figures .1,2, and 3, respectively. In-plane shear 
resistance is parallel to the direction indicated in each figure. With 
the variety of fabric constructional forms available, it is possible to 
have textile composite preforms behave according to any position 
desired on the performance maps illustrated in figures 1 through 3, in 
general . ( 9 ) 

When made from high modulus or fully drawn yarns, most multiaxial 
fabric constructions tend to have excellent dimensional stability and 
outstanding in-plane shear resistance in all directions (machine, 
crosswise, and bias) compared with regular woven, knitted, and braided 
constructions. Normally, dimensional stability and conformability are 
mutually exclusive. However, when made from partially oriented yarns 
(POY), multiaxial fabric constructions can be designed for directional 
conformability during 3-D draw-molding, 3-D draw-stamping or 3-D 
draw-wrapping with’ exce llent omni-directional dimensional stability and 
shear resistance in the final 3-D configuration. 

In many composite applications, isotropy (for performance) and 
porosity (for ease of resin impregnation) are quite important in 
composition preforms. A performance map plotting in-plane isotropy 
versus porosity for a variety of fabric constructional forms is 
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illustrated in Figure 4. It can be readily seen that warp knitting 
technology (regular, weft insertion, s t i tched- through and multiaxial) 
offers the greatest range of design possibilities. 
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EXPLORATIONS IN STITCHING TECHNOLOGY 

A number of multiaxial sti tch- through samples have been made from 
a variety of high modulus yarns in recent years by hand orienting and 
stitching to simulate machine-made constructions. Typically, 0°/-45°/ 
+45°/90° orientations have been made in 4, 8, 12, and 16 layers on 
tension frames of various dimensions and chain stitched at 2 or 3 wales 
per inch simulating a machine stitch-through operation. Yarn impaLing/ 
filament damage during hand stitching remains a problem especially 
with high density multiaxial orientations. 

Most of the hand-made multiaxial s t i tch- through samples have been 
made with chain stitching in the 0° (simulated machine) direction as 
straight and as uniformly as possible. Such an approach is appropriate 
for making large, thick, dense multiaxial preforms for planar 
applications requiring substantial, through- the- thickness strength and 
in-plane shear resistance. Currently, irregular stitching patterns 
with various stitching materials in various paths and angles are being 
explored in order to improve the conformabi lity of large, thick, dense 
multiaxial s t i tch- through samples in adapting to non-planar applica- 
tions . 

In the analysis of 3-D textile composites, Norris Dow showed 
clearly that improvements in a given property are usually accomplished 
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He concluded that 


at the sacrifice of some other properties (10). 
there is no magic configuration. Just as traue-offs are required to 
minimize the penalties for enhanced through-the- thickness properties, 
trade-offs are also required to minimize the penalties for enhanced 
conformability. Net shape textile preforms have received a great deal 
of attention in recent years, and rightly so. perhaps advanced 
stitching technology will make as great an impact on net shape textile 
preforms in the near future as it has on planar multiaxial preforms in 
the recent past. 
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MACHINE DIRECTIONAL BEHAVIOR Or FABRIC 



IN-PLANE SHEAR RESISTANCE 


Figure 1. Behavior of Various Fabric Structural Forms 
Under Uniaxial Stress In the Machine (0*) Direction 


CROSSWISE DIRECTIONAL BEHAVIOR OF FABRIC 



IN-PLANE SHEAR RESISTANCE 


Figure 2. Behavior of Various Fabric Structural Forms 
Under Uniaxial Stress In the Crosswise (90*) Direction 
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Flow Behavior in Liquid Molding 

D. Hunston, F. Phelan, and R. Pamas 
National Institute of Standards and Technology 
Polymers Division 
Gaithersburg, MD 


ABSTRACT 

The liquid molding (LM) process for manufacturing polymer composites with structural 
properties has the potential to significantly lower fabrication costs and increase production rates. 
LM includes both resin transfer molding and structural reaction injection molding. To achieve 
this potential, however, the underlying science base must be improved to facilitate effective 
process optimization and implementation of on-line process control. The National Institute of 
Standards and Technology (NIST) has a major program in LM that includes materials 
characterization, process simulation models, on-line process monitoring and control, and the 
fabrication of test specimens. The results of this program are applied to real parts through 
cooperative projects with industry. The key feature in the effort is a comprehensive and 
integrated approach to the processing science aspects of LM. This paper briefly outlines the 
NIST program and uses several examples to illustrate the work. 


INTRODUCTION 

Polymer based composite materials have much to offer in a wide range of products, 
particularly in structural applications. Despite the advantages, however, composites have realized 
only a small fraction of their potential. The most important barrier to their expanded use is the 
inability to make them rapidly, reliably, and cost effectively. Rejection rates in manufacturing 
are often too high, and the labor intensive manufacturing methods now used are slow and costly. 
The solution to this problem is the implementation of more effective processing methods and the 
development of the scientific understanding necessary to optimize their use. 

The process known as liquid molding (LM) has great potential in this regard. LM 
includes both resin transfer molding (RTM) and structural reaction injection molding (SRIM), 
which is a high speed cousin to RTM. LM combines some of the speed advantages of simple 
injection molding with the ability to make the high performance parts associated with continuous 
fiber reinforcement. Flexibility in the type, amount, and orientation of the reinforcement in every 
section of the mold enables the generation of very complex parts. Moreover, since the resin 
viscosity is low, the parts can be quite large and three dimensional. 

Because of these potential advantages, LM is the leading candidate for the fabrication of 
structural composite parts in the automotive industry. Likewise, the drive for more cost effective 
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manufacturing has caused a variety of other industries including aerospace and marine to consider 
LM. Effective utilization of this processing method, however, requires developments in both 
preform fabrication and the technology required for process optimization and on-line process 
control. Because this technique has such great potential, the National Institute of Standards and 
Technology (NIST) has initiated a major program to address the scientific issues associated with 
process optimization and control. This paper will give a general overview of the NIST program 
and then illustrate the work with two examples from the program. 


NIST PROGRAM OVERVIEW 

The LM program at NIST has four tasks which are summarized in Table I. Each will be 
briefly discussed in the material that follows. 


Table Is 

NIST Program on Liquid Molding 


Materials Characterization 
Permeability 

Thermal and Cure Behavior 
Process Simulation 

Macroscopic Models 
Microscopic Models 
Process Monitoring / Control 
Flow Behavior 

Thermal and Cure Monitoring 
Sample Preparation and Application to Parts 


Materials Characterization 

The first task focuses on characterization of the material properties associated with 
processing, i.e. preform permeability, thermal conductivity, and cure behavior. Preform 
permeability measures the resistance o ffered by the preform to th e flow of a resin. Since the 
reinforcement can have very different resistances to flow in different directions, the permeability, 
K, is a tensor quantity. This means a variety of measurements both in the plane and through the 
thickness of the preform material must be conducted to evaluate K. The resistance to flow is 
also very sensitive to the fiber volume fraction so this dependence must be determined. The most 
direct procedure to measure permeability is unidirectional flow experiments, and a variety of 
special molds have been developed at NIST to minimize edge effects and facilitate measurements 
at a variety of fiber volume fractions. For those cases where the part of interest has a shell-like 
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structure which permits the use of two dimensional flow models, it may be adequate to determine 
only the in-plane components of the permeability, and this significantly simplifies the 
characterization. 

This task also has the goal of developing the technology to predict the permeability from 
a knowledge of preform microstructure and fiber surface treatment This would be a very 
significant advance because the preform could then be designed to optimize both performance 
and processibility. Although this technology is well beyond current capabilities, it is an 
important long range objective. An example of research in this area is given later in this paper. 

The measurement of thermal and cure properties for the materials involved relies heavily 
on the process monitoring facilities established at NIST over the last five years. The objective 
in developing these facilities was to take advantage of NIST’s position as an outstanding 
measurement laboratory by assembling a wide range of process measurement techniques 
including virtually all the commonly used methods. Table II lists the ten different techniques that 
have been adapted to process monitoring in this program. This capability permits the 
examination of the chemical and physical changes that occur during processing at size scales 
ranging from individual chemical bonds where the chemistry occurs up to bulk properties such 
as viscosity and viscoelasticity. By applying these techniques both individually, and in 
combinations, a detailed picture of the changes can be achieved. 


Table n: PROCESS MONITORING TECHNIQUES 


Size Scale 


Chemical Bond 

Molecules 

Bulk 

Optical Spectroscopy 
Calorimetry 

Conductance 
Dielectric Spectroscopy 
Fluorescence Spectroscopy 
Chromatography 
SANS 1 

Ultrasonics 
Dynamic Mechanical 
Viscosity 

Measurement Methods 


‘Small Angle Neutron Scattering 


Process Simulation 

The second task in the NIST program is the development of process simulation models. 
The models fall into two somewhat arbitrary categories, macroscopic and microscopic. The 
macroscopic models employ the so called volume averaging approach. This approach focuses 
on volume elements that are large enough so that variations in local features such as the 
arrangement of individual fibers in space, interactions between fluid and fibers, etc. average out 
and thus only the average properties need to be considered. Such treatments are generally quite 



good for analyzing macroscopic events such as mold filling. 

The macroscopic models are less useful, however, for predicting other important events 
like void formation since these events depend on local features. To deal with this, microscopic 
models are needed. Such models include some or all of the local features. In principle, it is 
possible to simulate the entire part at this level of detail, but the computational time required is 
prohibitive. As a result, events dependent on local features are simulated by first using 
macroscopic models of the part to predict boundary conditions on the local area and then 
simulating this area using a microscopic model. 

To analyze complex geometries, the process simulation models are applied using finite 
element methods. In the NIST research effort, a variety of mold filling programs have been 
developed (1). The simplest is a 2D program, which is appropriate for flat, thin structures. The 
most general is a fully 3D program which is used for complex parts with thick sections. 
Bridging the gap between these is what is called a 2.5D model. This is for flow in shell-like 
structures in which locally there is always a plane with a thickness dimension that is small 
compared to the dimensions in the local plane. In these parts, it is assumed that there is no flow 
in the thickness direction, and thus, the flow is treated as two dimensional even though the 
structure is three dimensional. Both the 2D and 2.5D programs have significant computational 
speed advantages relative to the fully three dimensional program, which is the primary advantage 
of their use. The 2.5D and 3D simulations both account for the effect of gravity on the flow. 

Process simulations must analyze many factors including resin flow, heat transfer, 
chemical reactions, etc., as well as the interactions between these factors. Chemical reactions, for 
example, generate heat which must be considered in the heat transfer relationships. In certain 
cases, it is possible to separate some of these factors and thereby simplify the modeling. Often, 
the mold filling may be completed before the chemical reactions produce significant effects. 
When possible, the use of such simplifications can significantly accelerate the analysis. 


Process Monitoring / Control 

The third task in the NIST program involves on-line process monitoring. The purpose 
is to develop this technology and then use it to address two areas. First, the experiments seek 
to test and refine the process simulation models. Second, the technology for on-line process 
control is explored for those situations where it is appropriate. For testing the simulation models, 
the first step is an examination of flow behavior in mold filling. The principle tool is flow 
visualization experiments. One drawback of this approach, however, is that visualization which 
is performed with molds having one or more clear sides provides information only about what 
is happening at the surface. Consequently, sensors such as fiber optic probes are being developed 
to measure flow front position inside the preform. Such sensors may also have the potential to 
measure flow velocity and pressure if properly optimized for LM experiments. The area of cure 
monitoring is also an important part of the program. Cure monitoring has been a major focus 
of NIST’s research for some time. A detailed examination of the monitoring techniques listed 
in Table II was conducted to identify those that have the potential for on-line measurement. One 
result of this effort was a recent report (2) which provides an assessment of the state-of-the-art 
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process monitoring sensors for polymer composites. The preparation of this report was 
greatly facilitated by the fact that NIST had direct experience with all of the techniques involved. 
This capability will now be exploited through application to liquid molding. 


Sample Preparation and Application to Parts 

The final task involves the fabrication of samples for testing in another NIST program 
which considers the performance and durability of polymer composites. In addition, this task is 
applying the developments in the LM program to the study of real parts through cooperative 
efforts with industry. The best example is a program between NIST and the Automotive 
Composites Consortium (ACC). The ACC is a joint effort between Ford, Chrysler, and General 
Motors and was formed to conduct precompetitive research that promotes the use of composites 
for structural applications in automobiles. The ACC is demonstrating the results of their research 
by fabricating a series of parts. The first part chosen was the front end structure of the Ford 
Escort, and the fabrication method of choice was SRIM. 

NIST is cooperating with the ACC by conducting process simulations to help optimize 
the fabrication of this part. In addition, the permeabilities for the reinforcement materials used 
in the part are being measured at NIST. The ACC benefits by the development of the technology 
to optimize part design and manufacturing. NIST benefits by having the opportunity to test 
research results with application to a real part. The initial process simulation focuses on mold 
filling and uses a 2.5-dimensional model with gravity since the part is a shell-like structure. The 
complexity of the reinforcement, however, may make it necessary to employ a fully three 
dimensional flow model eventually. To better visualize the mold filling, a computer animation 
program has been developed. It allows examination of the mold filling from any angle and 
permits halting the flow temporarily at any point for closer examination of the details of the 
filling pattern. Important events such as the closing off of an unfilled area and the formation of 
knit lines can be seen and studied much more easily with the animation than with the normal 
output from process simulations. 


EXAMPLES FROM THE NIST PROGRAM 

To illustrate the NIST program, two examples will be briefly outlined. More details on 
these particular studies can be found in references 3 and 4. The first example looks at the 
prediction of permeability and illustrates the use of microscopic flow models and how they can 
help to examine the effects on flow behavior of localized features such as the distribution of 
fibers in space. The second example examines two techniques that are commonly used to 
measure permeability and compares the results. This work illustrates the research in both 
materials characterization and flow visualization. 


Permeability Prediction 

Macroscopic models for fluid flow in LM are usually based on Darcy’s law 
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( 1 ) 


<v> = - 


jrv<p> 


where <v> and <P> are the volume averaged velocity and pressure, respectively, K is the 
permeability of the porous medium, and |i is the viscosity. These models are the basis for most 
process simulation programs including those developed at NIST (I). One of the most important 
factors affecting the agreement between the actual filling of a part and the predictions of such 
process simulation programs is the accuracy of the permeability data available. At the present 
time, the best and essentially the only method for determining the permeability is through 
experimental measurements (3,5-7). While this is satisfactory in the sense that various techniques 
give consistent results (3), and hence, appear to be reliable, it is also desirable to have available 
theoretical models which can accomplish this. Such models would not only permit optimization 
of preform structure but also greatly increase the speed with which LM processes could be 
implemented and/or modified by eliminating (or reducing) the large number of carefully 
controlled experiments that are currently required for determining the necessary permeabilities. 


A number of theoretical models for predicting permeability as a function of structure in 
fibrous porous media are available (8-10). The most widely used relation in composite 
applications is the Carman-Kozeny equation which for fibrous porous media is given by 


K _ e 3 d 2 
(l-e) 2 16c, 


( 2 ) 


where e is the porosity, d is the diameter of the fibers which make up the media, and c K is the 
Kozeny constant which has a theoretical value of 2. While some studies have reported success 
with this equation (11,12), generally such success is over a limited porosity range with a value 
of the Kozeny constant much different from its theoretical value (13).* Such limited success is 
somewhat typical. Therefore, at this time, there is not available a general relation which 
describes the permeability of fibrous porous media over the entire range of porosities that are of 
interest 


In order to aid the development of theoretical models for permeability, various numerical 
studies have been undertaken with microscopic models (14-17*). In such studies, flow in ideal 
periodic arrays of cylinders (or some other appropriate structure) is modeled, from which 
predictions for permeability as a function of structure are calculated. Corresponding experiments 
on systems which closely resemble the idealized arrangements have also been conducted. A 
review of the relevant experimental and numerical results for fiber porous media has recently 
been completed (13). An interesting result for the case of axial flow through unidirectionally 
aligned fiber beds is that experimental measurements (11,18) consistently yield higher 
permeability values than those predicted by simulation (15). The same result holds true for the 
case of transverse flow (13). Thus, there is difficulty in matching experimental results not only 
with theoretical relations such as Carman-Kozeny but also with detailed simulation. 

One possible reason for this discrepancy between simulation and experiment is that the 
materials used in LM preforms are heterogeneous in the sense that the network of fiber bundles 

* personal communication. 
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(tows) which make up the global porous medium are themselves porous. Thus, in actual 
materials, there is flow not only around but through the tows. The effect of this heterogeneous 
structure on unsteady flow transients, void formation and ultimate properties has recently been 
discussed (19,20) and shown to be quite important. Therefore, it seems reasonable to investigate 
the effect of this property on the permeability as well. 

The numerical study of the microscopic flow in fibrous porous media is complicated by 
the fact that this problem involves both open regions (the space between tows) and porous media 
(the tows themselves). The flow in the open media can be modeled using the Navier-Stokes 
equation. The flow in the porous media might be modeled using Darcy’s law, but with this 
approach it is not possible to properly handle the interface between the two regions because the 
order of Darcy’s equation does not allow for the proper specification of boundary conditions (i.e., 
continuity of velocity and stress) at interfaces. The approach used in this study to achieve the 
correct degree of continuity is to model the flow in the porous region using the Brinkman (21,22) 
form of the volume averaged momentum equation (denoted hereafter as BVAME) which is given 
by 

pV^v) - pJC'Mv) = V(p) (3) 


It is interesting to note that if only the first term on the left were present, the result would 
be mathematically equivalent to Stokes equation while if only the second term on the left were 
present, the relationship would be identical to Darcy’s law. The most important feature of 
Brinkman’s equations, however, is the presence of second order terms. When this relationship 
is used to describe the flow in the porous region, the second order terms make it possible to 
satisfy the boundary conditions at the interface with the open region where Stokes law is used 
to describe the flow. A variety of problems involving heterogeneous media have been 
successfully treated using this equation (23,24). While more general expressions than Eq. (3) are 
available (11,15,25,26) Eq. (3) represents the simplest and most easily applied form. 

Simulation 


The BVAME was used to develop a microscopic model for flow to elucidate the effect 
that tow permeability had on the overall permeability in fibrous porous mediums. The axial flow in 
a periodic, square array of cylinders acting under the influence of a constant pressure gradient 
of value -G was simulated for the case of both solid and porous cylinders. The computational 
domain for these geometries is depicted in Figure 1. The problem was scaled using the cylinder 
diameter, d, as the characteristic length, and a characteristic velocity of d 2 G/ji. The relevant 
component of Eq. (3) is then 


&U f &U a u= ± 
dX 2 dY 2 K* 


( 4 ) 


where X, Y, U, and K* are the dimensionless coordinates, velocity, and permeability, and a is a 
computational parameter equal to 1 in porous media, and 0 in open media. Solutions to this set 
of equations using the appropriate boundary conditions were generated using a Galerkin finite 
element scheme. Details of the boundary conditions and numerical formulation are discussed in 
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detail elsewhere (4). 

Calculations were performed for the case of solid cylinders first to verify the numerical 
scheme by comparing with other calculations, and also, to provide a reference point for looking 
at porous cylinder results. Figure 2 summarizes the computations by showing the dimensionless 
effective permeability (i.e., K/d 2 ) as a function of the porosity. Although not shown in the 
Figure, the results from the calculations performed here agree quite well with the computations 
by Sparrow and Loeffler (27). For the sake of comparison, the predictions of the Carman- 
Kozeny equation and experimental data of Gutowski et al. (1 1) are also plotted. The Carman- 
Kozeny results match the numerical calculations only at intermediate to low porosities, and 
deviate strongly at the extremes. This is similar to what is generally found when comparing 
Carman-Kozeny to experimental data for systems with solid fillers (13). The most significant 
observation from the Figure, however, is that the experimentally determined values of 
dimensionless permeability with porous tows are much higher than any of the numerically 
calculated and theoretical curves for solid cylinders. 

Calculations for porous cylinders were performed with cylinder permeabilities of 0. 1, 0.01, 
0.001, and 0.0002. Iso- velocity contours (at the same contour values) for solid and permeable 
cylinders at a cylinder spacing corresponding to an overall porosity of 0.65 are compared in 
Figure 3 for the case of a cylinder permeability equal to 0.001. The Figure shows that the effect 
of the cylinder permeability is to shift the contour values inward towards the axis of the cylinder. 
This shifting results in an increase in die overall flow rate throughout the system. 


i 

■ 

I 

i 

I 

| 

| 

I 

1 

1 





Axial Flow, Iso- Velocity Contours 



(a) Solid 


(b) Porous 


Porosity = 65% 


Figure 3 : Iso-velocity contours for axial flow through periodic, square arrays of (a) so l i d and 
(b) porous cylinders. There are ten equally spaced contours. The maximum velocity for both 
cases occurs at the upper right-hand comer. The permeability of the porous cylinder is 0.001. 



permeability, was computed for cylinder spacings corresponding to porosities of 0.9 and 0.65. 
This is shown in Figure 4. As was expected from inspection of the velocity contours, the results 
show that there is a significant flow rate enhancement due to cylinder permeability. The level 
of enhance ment prope rly converges towards zero as the inner cylinder permeability decreases. 
Perhaps the most interesting feature of Figure 4 is that it shows that the flow rate enhancement 
due to the permeability of the cylinder becomes more important as the cylinder volume fraction 
increases. Thus, even at a cylinder permeability of 0.0002, the flow enhancement is still over 
13% (see Table III). 

When examining these results, it is important to note that although there is flow through 
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Fiber Tow Permeability 


Figure 4: The percentage increase in the flowrate outside the porous cylinder with 
respect to the case of solid cylinders as a function of tow permeability. The effect 
becomes more pronounced as the overall porosity decreases. 
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the porous cylinders, this flow is quite small, and in fact, almost negligible compared to the 
overall flow rate. The real mechanism for the flow increase is the change in boundary conditions 
at the cylinder surface which shifts all of the flow contours in Figure 3 toward the center of the 
porous cylinders relative to the case for solid cylinders. This causes an increase in velocity 
throughout the entire flow field. Note particularly the increase in the size of the cross-sectional 
area at the highest flow velocity (upper right comer in Figure 3). Since the flow rate is the 
integral of the velocity taken over the entire flow domain, the overall change in this quantity is 
quite substantial. 


Table III: Flow Rate Enhancement 

vs Fiber Tow Permeability 

Fiber Tow 

Flow Rate Enhancement (%) 

Permeability 

Porosity = 0.90 

Porosity = 0.65 




0.1 

225.0 

749.6 

r 0.01 

43.7 

128.8 

0.001 

12.1 

32.7 

0.0002 

5.1 

13.6 


Summary 

The calculations presented here indicate that the heterogeneous nature of the LM preform 
material can have a significant effect on the permeability. The calculations predict that the flow 
rate is substantially higher for porous cylinders than for solid cylinders at equal cylinder spacing 
and that this increase becomes more pronounced as fiber volume fraction increases. This is 
consistent with results in the literature where experimentally determined values of permeability 
in such model networks are generally found to be higher than their theoretical counterparts (13). 
The results here show that in such situations, the permeability of the tows provides a potential 
mechanism for flow rate enhancement, and hence, an increased permeability. At this time, 
calculation of an effective permeability for the overall medium for the case of porous cylinders 
has not been attempted. This will be done in the future, but requires intra-tow permeability vs. 
porosity data which is not presently available. Nevertheless, the results obtained here suggest 
that in order to develop theoretical relationships for permeability as a function of structure in 
fibrous porous media, it is important to take into account the intra-tow properties of the media. 


Permeability Determination 


The second example from the NIST program concerns the measurement of permeabilities. 
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Although all components of the permeability tensor can be important, the discussion here will 
focus on the in-plane components. There are two techniques commonly used to measure the in- 
plane permeabilities of preform materials. The first involves unidirectional flow while the second 
focuses on radial flow. As part of a major effort on permeability characterization, the NIST 
program is working in cooperation with General Electric’s Central Research and Development 
Laboratory to conduct a detailed comparison of results from the two methods. This tests both 
the experiments and the flow model used to compare the results. In addition, the data are being 
examined in terms of the structure of the preform material in an effort to understand the observed 
behavior. The initial experiments involved measurements on three woven glass reinforcements: 
J.P. Stevens 8-Hamess satin, and Carolina Narrow Fabrics 8-Hamess satin and Crowfoot 
weaves. 


Characterization of the in-plane permeability requires the determination of three 
parameters: the two principle permeabilities, K n and K 22 , and the principal angle, 0, i.e the angle 
between the warp direction of the preform material and the direction of the larger of the two 
principal permeabilities. Moreover, the values obtained for the permeabilities are different for 
measurements of flow into a dry preform and for flow into a preform completely saturated with 
the fluid. Although explanations for this difference have been offered (19), experimental 
measurements for both are still required for a complete characterization. 

One-Dimensional Flow Experiments 

The one-dimensional flow experiments were conducted in a mold with a 152.5 mm x 
152.5 mm (6 in. x 6 in.) sample section where the thickness of the test specimen could be 
adjusted between 3.175 mm and 12.7 mm (0.125 in. and 0.5 in.) (see Figure 5a). The top of 
the mold was 12.7 mm (0.5 in.) thick plexiglass to permit observation of the flow front, and the 
flatness of the plexiglass surface was checked during the experiments to detect and prevent 
problems due to pressure driven deflection. Fluid was forced through the mold from a pressure- 
pot, and the flow was measured as it entered the mold with a calibrated rotameter and was also 
monitored at the mold exit with volumetric measurements. Pressures at the mold inlet were 
measured by transducers that were mounted at two positions on the back face of the mold located 
at the lower boundary of the reinforcement sample. 

The results from the one-dimensional flow experiments were accumulated by measuring, 
at preset pressure drops, the flow of com syrup/water solutions through the woven glass fabric, 
compressed to 12.7 mm (0.5 in.) thickness and previously saturated with the test fluid. As the 
pressure drop was increased, the flow rate increased linearly (see Figure 6), as expected, for the 
Newtonian test fluid. Assuming Darcy’s law to be valid, the effective permeability of the 
preform in the flow direction was obtained from the slope of the flow-pressure curve. The linear 
relationship of the flow and pressure drop qualitatively indicates that pressure induced bending 
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Schematic diagrams of the 


equipment and (b) radial flow equipment. 


experimental equipment: (a) one-dimensional flow 


I 


of the upper plexiglass plate did not occur. This provides a very accurate way to determine I 

permeability in saturated flow. Since the top of the mold is plexiglass, experiments can also be 
performed to measure the progression of the flow front into a dry preform. Tlius, although the m 

experiments here focused on saturated flow, the permeabilities for both saturated and unsaturated = 

flow can be obtained with this technique. - 

Radial Flow Experiments S 

The radial flow experiments were conducted at General Electric’s Central Research 
Laboratory using the apparatus shown in Figure 5b (3). A radial planar flow was achieved by - 

injecting the fluid through a central 1.524 mm (0.060 in.) diameter gate into a 203.2 mm x 203.2 | 

mm (8 in. x 8 in.) region between two parallel plates containing the reinforcement. Each layer " 

of reinforcement had a 4.762 mm (0.1875 in.) diameter hole centered over the gate to permit the ■ 

fluid to penetrate evenly through the thickness of the sample. A 25.4 mm (1.0 in.) thick upper | 

glass plate permitted visual observation of the flow. A central hydraulic system supplied both 1 

the pressure to clamp the reinforcement sample and pump the test fluid from a separate injection 1 

cylinder. I 

After an initial transient, the flow front maintained a shape and orientation that appeared 1 

qualitatively unchanged as it grew. Traces of the flow front were made manually, with a typical K 

example shown in Figure 7 for the J.P. Stevens 8-Hamess satin weave fabric. Although B 

relatively crude, this presentation technique illustrates important general features of the flow S 

front The shape of the flow fronts were analyzed by digitizing the traces. For a preform that 1 
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Flow, ml/s 



Figure 6 : Typical results from a one-dimensional flow experiment showing flow vs. pressure 
drop for two fabric orientations. 


is isotropic in the plane, a circular pattern would be expected. If the permeability in the plane 
is anisotropic but orthotopic, an elliptical shape should fit the data. For the materials examined 
here, the flow front was indeed elliptical in shape, and this confirms orthotopic behavior for the 
high volume fractions tested. The direction of the principal axis of the ellipse relative to the 
warp fibers in the preform provides a measure of the principal angle, while the square of the ratio 
of the major to minor axis lengths in the ellipse is equal to the anisotropy which is the ratio of 
the two principal permeabilities, ¥^JK n . 

As can be seen, the determination of the principal angle and the anisotropy in the radial 
flow experiment is quite easy. Although it is also possible to determine the quantitative values 
of the permeabilities, this is more difficult. Moreover, the experiments described above measured 


37 




Fill 


Figure 7: The elliptical flow fronts observed in J.P. Stevens 8-Hamess satin fabric. 


unsaturated flow properties. Although permeability measurements for saturated flow are possible 
in the radial flow experiments, they are considerably more difficult. Consequently, the work here 
focused only on the determination of the principal angle and anisotropy for unsaturated flows and 
compared the results with values for these parameters obtained from saturated, unidirectional flow 
experiments. 

Comparison of Results 

To compare the one-dimensional flow results with the radial flow results, it was necessary 
to calculate the principal permeabilities and the principal angle from the unidirectional flow 
experiments. By using a two-dimensional version of Darcy’s Law, the effective permeability, 
K e , in a direction which is at an angle, <J), relative to the warp direction in the preform can be 
written in terms of the principal permeabilities and the principal angle. 

K% = JC^cos 2 <<|> — 8) • 


(l + ^i tan 2 («j>-e)) 

K n 


(*22/^1 - 1 ) 2 
*2 2 /*ll , 1 

tan 2 (<J>-0) 


(5) 
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Since this equation involves three unknowns, measurements at a minimum of three 
different orientations must be made. To improve reliability, four orientations were tested in this 
work. The results of measurements conducted at two orientations are shown in Figure 6. That 
figure indicates that the permeability is anisotropic which agrees with the radial flow experiments. 
From the four experiments on each material, the values for the principal permeabilities and the 
principal angles were calculated for the three materials tested. The results are given in Table IV 
along with the data obtained in the radial flow experiments. Additional details of the analysis 
are given in reference 3. 

The consistency of the radial flow experiments performed at GE and the one-dimensional 
flow experiments conducted at NIST is illustrated in Table IV. Even though two different 
laboratories conducted the experiments, the agreement is quite good. As the experimental 
techniques are further refined, it will be important to see if the agreement gets even better. 


Table IV: Comparison of 1-D and Radial Flow Experiments 

Material 

Laboratory 

Principle Angle 

W K„ 

J.P. Stevens 
8-H satin 

GE 

12 

0.64 

NIST 

13 

0.72 

Carolina Narrow 
Fabric 
8-H satin 

GE 

-17 

0.77 

NIST 

-15 

0.73 

Carolina Narrow 
Fabric 
Crowfoot 

GE 

3 

0.69 

NIST 

4 

0.61 


It is interesting to note that the principal direction of the permeability is not coincident 
with either the warp or fill yam directions. An hypothesis for this behavior has been developed 
based upon a closer examination of the structure of the woven fabrics. The crossing of the fiber 
tows creates a series of valleys in a woven material known as the crimp. Figure 8 illustrates this 
structure for the J.P. Stevens 8-Hamess satin fabric, and the Carolina Narrow 8-Hamess and 
Crowfoot fabrics used in the experiments reported above. The close proximity of neighboring 
crossing points creates a nearly continuous trough that may act as a conduit of low resistance to 
flow. 


Note that the orientation of the illustrated crimps in Figure 8 match the orientation of the 
principal axes of permeability observed in both the radial and one-dimensional flow experiments 
with the two 8-Hamess satin weave fabrics. A secondary, nonorthogonal crimp direction does 
not appear to influence the flow because the tow crossing points are not on adjacent tows. 
However, a crimp direction in the Crowfoot weave could not be identified for comparison with 
the principal axis of permeability in that material. An additional unresolved point in this 
discussion is that the major crimp directions on opposite sides of the fabrics are not oriented 
identically, creating a complex interface between fabric layers which must be explored further. 
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Figure 8 : Structure of woven fabrics showing crimp directions. The angle between the crimp 
direction and the warp fiber direction is indicated. Note that a crimp direction is not apparent 
in the Crowfoot weave. 


The data developed in this project represent a first step in exploring not only the test 
methods for permeability measurement but also the models used to analyze the data. Constitutive 
models, such as Darcy’s Law, have largely been used without detailed verification. The limited 
results obtained so far not only indicate that the unidirectional flow and radial flow experiments 
are consistent but also support the use of Darcy’s Law in flow simulations. Since the radial 
flow experiments quickly provide a measure of flow anisotropy and principal orientation while 
the one-dimensional flow experiments easily yield quantitative measurements of the permeability, 
the results also suggest that the two experiments could be used in combination to simplify the 
characterization of preform materials. 
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CONCLUSIONS 


The program described here seeks to advance the scientific basis for liquid molding 
through a coordinated effort involving four tasks: (1) characterization of the material properties 
associated with processing, (2) development of process simulation models at both the 
macroscopic and microscopic levels, (3) application of on-line measurement techniques to test 
and refine the simulation models and to explore the technology for on-line process control, and 
(4) production of test samples and the application of the technology developed in the program 
to real parts through cooperative interactions with industry. Such a comprehensive and 
interdisciplinary approach offers the best hope for achieving the full potential of the liquid 
molding process. 
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ABSTRACT 


This paper encompasses an overview of the history of 3-D 
Braiding and an in-depth survey of the most recent, technological 
advances in machine design and implementation. Its purpose is to 
review the major efforts of university and industry research and 
development into the successful machining of this textile process. 


INTRODUCTION 


The technology of 3-D Braiding epitomizes the evolution of 
ancient art and custom into the science of engineering a new, 
innovative textile process. Unlike other textile techniques like 
weaving and knitting, 3-D Braiding involves the continuous 
intertwining of two or more yarn groups in a bias direction, and the 
structure formed possesses remarkable structural integrity, high 
damage tolerance, torsional stability and characteristic handling ease 
[16]*. As this technology expands from university research to 
industrial manufacturing and development, careful investigation into the 
design and automation of successful braiding machines is imperative. 
Furthermore, the future growth of this technology depends on the 
ability of investigators to improve upon and modify processing 
parameters while reducing any process limitations. Accordingly, this 
paper addresses current advances in the area of 3-D braiding 
machines, focusing on the three main design areas: machine bed 

configuration, carrier design, and beat-up mechanisms. These three 
components are instrumental to the overall feasibility of a design for 
a modern 3-D braiding apparatus. 

* (References 1-24 are cited in the text) 
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HISTORY 


The origin of 3-D braiding can be traced to knot-tying, which is 
one of the oldest technologies known to man [17]. Since the beginning 
of time, man has ritualistically twisted and plaited hair for 
convenience, adornment and customs. This decorative styling of the 
hair developed into the practice of intertwining grass and natural 
fibers to form ropes and ornate embellishments in artwork and 
clothing design. In the ancient Orient, skilled artisans crafted 
intricate braids called Kara Kumi for ornamental purposes [21]. 
During the 17th century, English sailors devised spool braiders that 
formed braids for practical rather than artistic use. The spool braider 
or sinnet made braids that served a useful purpose as durable ropes 
[21]. The formation of braids for scientific rather than artistic means 
underwent a slow evolution until the mid-twentieth century. 

During the next phase of development of braiding as an 
advanced technology, the practice of braiding was no longer 
considered only a crafting art. In the 1940's, W. T. Harburger 
authored a paper on the engineering structure of braids [17]. His paper 
addressed geometric factors relating to the concept of braiding. In 
the next decade, D. Brunnschweiler discussed the configuration, 
geometry, and tensile properties of tubular braids [17]. In the 
1960's a plethora of new concepts involving the mechanisms of the 
3-D braiding process evolved. An informative book by W. A. Douglas 
dealing with braiding machinery and processing parameters became 
a cornerstone for the growing technology [17]. Comprehensive 
patents by Bluck and General Electric established 3-D braiding as a 
technical process in the field of science and engineering. Emerging 
from this period of growth and innovation were two fundamental 
classifications for the braiding process, four-step and two step. Both 
classifications share two distinct machine types, regular and circular 
[13]. The former produced braids having right-angled corners, 
whereas the latter produced round cross-sectional shaped braids. In 
1965, General Electric Reentry Systems Divisions invented the 
Omniweave braiding process for the production of multidirectional 
braids for aerospace composites. The 4-step machine incorporated 
sequential, discrete carrier motion as opposed to continuous motion 
[17]. GE claimed that the Omniweave would vary fiber orientation 
while braiding stiff, brittle fibers without delamination failure [13], 
In 1969 the Bluck patent on "High-Speed Bias Weaving" detailed a 



circular 4-step braiding machine for the production of hollow braids 
[13]. On the machine spools sitting beneath the braiding plane 
supplied yarns. During the process, these yarns moved through a 
network of holes in a guide nest which circumrotated the machine. 
This process was faster than its predecessors; however, entanglement 
of feeding yarns due to guide nest movement was a major limitation 
[13]. Consequently, the braided product had a fixed length and yarn 
delivery was intermittent. 

In 1973 Maistre patented the first automated 3-D braiding 
machine. He promoted the new braiding process as SCOUDID, an 
acronym for Structure Composite Unidirectionelle Indelaminable 
(undelaminable, unidirectional composite structure) [13]. A regular- 
type, 4-step braiding machine, the SCOUDID had the capability to 

braid over 3200 yarns in simultaneous fashion. The yarn supply 
rose vertically and did not traverse the cross-section of the braid. 
The braided products had fixed lengths and two layers of yarn with 
adjoining edges [13]. 

In 1982, Florentine developed a versatile braiding machine for 
the production of complex-shaped 3-D preforms having integrated 
structural geometry and varying lengths. The Magnaweave, his 4- 
step regular braiding machine had a system comprised of 4 x 2 

motion due to row/column shifts actuated by pneumatic cylinders, 
and manufactured braids possessing "SCOUDID" texture geometry 
[13]. Also, a combing or compacting movement followed an 
orthogonal shedding mechanism providing a high degree of freedom 
in material orientation control [9]. The Magnaweave had a 

counterpart, the Magnaswirl. This counterpart was a circular braider 
producing braids with circular cross-sections by three distinct 
motions: (1) ring or angular, (2) spoke or radial and (3) vertical or 
combing [16]. The design of the machine consisted of radial rows 
where carriers sat and a base track for tangential motion of arcuate 
members. The partial occupancy of the radially positioned rows by 
carriers produced tight braided structures [13]. 

Contributing further to the evolution of the 3-D braiding 
process was the invention of Through - the - Thickness 1 m 

technology by Atlantic Research Corporation (ARC), which is 
considered in industry to be the leading authority on 3-D or Through 
-the - Thickness Braiding. The circular braider manufactured had 
concentric grooved rings and carriers capable of alternate row and 

column position shifts. During the process of Through - the - 
Thickness 101 braiding, a technique for manufacturing 3-D seamless 
patterns by continuous intertwining of fibers, the carriers situated on 
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the machine bed moved simultaneously [6]. The process was a major 
advancement exhibiting a novel textile technology of arranging 
fibers in irregular geometries for arbitrary thickness [8]. The 
braided products of this process had varying lengths, composite material 
toughness, and adequate delamination resistance. 

As the present decade approaches the 22nd century, the 
advance of 3-D Braiding is inevitable. The current trend is toward 
computer-aided design and manufacture of the braiding process. Ko 
and his associates devised a computer simulation of the 4-step 3-D 
braiding procedure [17J. This system screened braid geometries for 
their respective cross-sections. Then, Brown and Harman used 
computer aided design to simulate the operation of a braiding 
machine and to trace the resulting yarn paths of the braid [18]. 
Developing areas of machine design and exploring multidirectional 
braiding techniques for near-net shapes of preforms through the 
extensive use of CAD/CAM technology are the future of composite 
engineering. 


MACHINE BED DESIGN 


The machine bed, a major component of any 3-D braiding 
apparatus, supports the precise movement of fiber-carrying devices 
(carriers or braiders). Accordingly, the purpose of the machine bed 
is to facilitate and expedite carrier placement along predetermined 
paths within the braiding plane. In recent years many 

configurations for machine beds of regular, circular, and variations of 
both types have been designed and tested. 

The 3-D braiding group at North Carolina State University's 
(NCSU) College of Textiles designed and constructed a 2-step regular 
braiding machine. Its machine bed warrants attention due to its 
potential as a marketable advance. Consisting of an assembly of unit 
tiles made of aluminum plates embedded with orthogonal T-grooves, 

the 2-step braider utilizes a dual arrangement (Figure 1) [18]. The 
first prototype exhibited machine arrangement one, having machine 
sides parallel and perpendicular to slots. The prototype consisted of 
four (4) tiles. Each unit tile has dimensions of 12" x 12" x 1" and a 
distance of 3" between grooves (Figure 2) [19]. Each tile has a hole 

for axial inset tion from spools underneath the braiding plane. Next, 

for the second arrangement, which reduces machine space, machine 
sides are parallel and perpendicular to the preform sides, being at 45 
angles to the slots. 



The actuation system embodies an array of stiffened timing 
gear-belts that push the carriers. While inside the machine bed 
grooves, the belt behaves like a "rigid rod" [18]. The flexible joint 
action occurs once the belt emerges from the groove and drops under 
its own weight. Also, the driving system utilizes stepping motors 
[20]. Four of the eight stepping motors drive the timing gears to 
control carrier stroke. The other four motors transport the driving 

units carriages across the perimeter of the machine bed. 

The 2-step regular machine incorporates a computer terminal 
and controller as its control system. Utilizing a special motor control 
language, the braiding programs are written. The programs are 
transmitted to the CPU of the controller, while the APPCOR IMC-8 
controller guides the action of the eight stepping motors (Figures 3 & 
4) [19]. 

This prototype has many advantages. First, it can be easily 
manufactured and expanded as a result of its tile assembly 
construction of the machine bed. Second, the machine cost can be 
reduced because only eight stepping motors are needed regardless of 
machine size. Lastly, the mechanism of a horizontal rod driven by 
two air cylinders over the machine bed acts as a convergence device 
which eliminates the "over jamming" effect [18]. 

An emerging technology is in the area of multi-ply braiding. 
This concept involves the formation of braid with more than 2 layers 
(plies). Brown and Ratliff of ARC invented a technique to move the 
machine bed that was comprised of a sequenced braider motion for the 
multi-ply technology [7], Patented in 1986, this method includes a 
machine bed configuration of a row/column carrier arrangement. 
When an intermediate row moves to block the motion of a column, a 
tamping force causes the column to shift to one side of the selected 
row [7]. This sequence of action makes carrier alignment in columns 
simpler. 

This braiding innovation claims the following advantages: (1) 
machine jams are eliminated, (2) unjammed operation of a multi-ply 
braiding system is achievable and (3) lower cost, low tolerance 
multi-ply braiders may be used [7]. 

For circular braiders, row movement is caused by the shift of 
concentric rings. Column (radial) movement occurs when discreet 
carriers are shifted. A difficulty in machining of circular braiders 
arises due to the expense of producing properly fitting concentric 
rings [15]. Often, a circular braider is quite large, requiring ample 
space. Also jamming can occur when a large number of carriers shift 
radially. 


Consequently an easily expandable 2-step braiding machine in 
circular form having no concentric rings was next designed and built 
by El-Shiekh and his colleagues at NCSU. The machine bed has 12 
tiles arranged together to form a circular configuration (Figures 5 and 
6) 1 13]. Each tile, made of aluminum, consists of axial tubes, magnets, 
stoppers, and braiding carriers. The yarn supply for the axials passes 
from a point above the plane of braiding down through the axial 
tubes where they are suspended by weights. The spooled carriers 
maintain a supply of up to fifteen feet of yarn. 

On the machine bed a magnet attracts the carriers against yarn 
tensions. The action of the stoppers halts carrier movement 
maintaining adequate yarn tension. The bed has 48 slots in clockwise 

direction. Every slot is formed in the shape of an inverted "T" shape, 

which facilitates the in and out movement of the carriers [13]. 
Elimination of the jamming effect is a major advantage of this 
assembly. The carrier arrangement determines the required 
preform shape. 

An ARC circular braiding machine was patented in 1988 by 

Richard Brown. The apparatus bed consists of a plurality of 
interchangeable rings having the same diameter (Figure 7) [4]. The 
rings are situated side by side in an axial arrangement. These ring 
members hold carriers that move axially and rotate about the central 
axis of the machine. During the braiding process, the Brown 

apparatus intertwines yarn by shifting adjacent carrier rows in 
opposing directions while spinning adjacent ring members in 
opposite directions. Notably, an actuating mechanism causes the ring 
rotation and propels the carrier rows in axial fashion to form the 
braid. This mechanism comprises a "manually or power driven 
actuator" that is connected to the rings through gears [4], The 
actuation system for the carriers is composed of "slidable" rods or 
pistons which are situated at opposite carrier-row ends. These 
pistons may be driven pneumatically or mechanically. 

The cylindrically designed circular braiding machine is a 
noteworthy advancement due to the following distinctions [4]: 

(1) ring members have same dimensions for ease in 
constructing the machine bed 

(2) ring members are both expandable and interchangeable 

(3) the machine requires little space and 

(4) the machine can produce axi-symmetric and cartesian 
braids. 

ARC also patented a continuous circular braider comprised of a 
multitude of "flexible annular" members or belts. These belts 
support the carriers and are arranged in side by side fashion, axially 



aligned [10]. The belts glide circumferentially relative to a central 
axis. The carriers, mounted on the belts move as well. An actuation 
device drives the belts and carrier rows along axial routes 
intertwining the fiber into a braided structure [10]. During the 
braiding process, the motion of adjacent carrier rows in opposing 
axial directions and the movement of adjacent belts in opposite 
circumferential directions initiate braid formation. 

On this apparatus the mechanism of actuation consists of 
sensors mounted on the belts which activate engaging slots found on 
a sector-by-sector basis. After belt shifting occurs this actuating 

mechanism allows the sector positioning of sensors [10]. Next shoes, 
sliding mechanisms, are exchanged with belts in that one sector. This 
sequence of sensor positioning and shoe sliding repeats for each 
sector until the original sector is reached. 

The Culp patent makes the following progressive claims [10]: 

(1) a multitude of belts are arranged axially to support carriers for 
axial motion, 

(2) the machine track supports and guides the belts along a 
circular path 

(3) a group of carrier rows are situated on the belts 

(4) an actuation device moves the belts circumferentially and 
carrier rows axially. 

In 1989 Spain of Airfoil Textron Incorporated (ATI) patented a 
regular, 3-D braiding apparatus with a machine bed designed for the 
integration of braider and axial yarns. The machine surface consists 
of a grooved X-Y grid equipped with axial guides for the yarn supply 
lying beneath the grid on creels [23]. Inside the grooves sit the 
carriers whose movement is controlled by means of solenoids or fluid 
cylinders at opposing ends for each column and row (Figure 8) [23]. 

This patented design is noteworthy due to its inclusion of an 
innovative technique of fugitive (non-permanent) braider or carrier 
yarn removal following braiding. This procedure facilitates matrix 
impregnation, spacing control between non-fugitive (permanent) 
braider yarns and engineering of preform physical and mechanical 
properties [23]. In another embodiment of the machine bed for this 
design both axials and braiders are permanent, creating a resultant 
structure with improved mechanical properties [23]. Therefore, the 
Spain apparatus varies machine configurations in order to diversify 
the design potential for the resultant braided products. 

In 1990 Ivsan and his associates of ATI patented a braiding 
apparatus having a cylindrical surface unlike flat, planar machine 
beds. This cylindrical or quadratic carrier surface serves to minimize 
the distance between carriers and the convergence point of the 


forming braid (Figure 9) [15]. At opposite ends of the quadratic 
surface lie actuators of hydraulic or pneumatic form, which propel 
the carriers along support members or curved tracks comprising 
the bed [15]. The support members move sequentially as carriers 
transit adjacent tracks. A series of tubes on the tracks guide 'axial 
sluffer fiber strands' in a radial fashion [15]. 

The basic structure of the machine bed also includes stationary, 
independent base members that house spindles and guide rods [15]. 
In contrast to carrier movement the cylindrical bed does not allow 
independent axial movement relative to the support members [15]. 
The axial direction of motion forms a circular arc with respect to 
corresponding support members. The shifting of carrier paths across 
adjacent support members causes the intertwining of carrier and 
axial yarns to maintain proper yarn tensioning and convergence of 
the braid. 

Fabrication of a multi-layer interlocked braided preform 
requires a machine bed which unlike track and column machines 
does not limit the insertion of 0° orientation yarns. This process is 
faster, allowing complex carrier systems movement and reducing the 
mechanical efficiency necessary to automate the process [3]. During 
the multi-layer interlock braiding process, carrier movement is both 
smooth and continuous. David Brookstein of Albany International 
Research Corporation (AIRC) commented, "the primary structural 
characteristics of multi-layer interlock braids are the nature of 
interconnectivity of adjacent laminae" [2]. The resultant braid 
consists of yarn that traverses from an original braided layer to an 
adjacent layer and back to the original layer. 

AIRC has manufactured a multi-layer interlocking tubular 
braiding machine, equipped with five layered machine configuration. 
Each layer has the capability of braiding 48 yarns. A network of 
counter rotating horn gears propels the carriers along diagonal paths 
moving the yarns layer to layer of the preform [2]. The machine 
provides for up to 48 axials per layer. Thus, 480 yarns may be 
simultaneously braided. Another system developed by the AIRC 
group was designed to manufacture braided preforms of varying 
shapes. Cross-sections in the form of C, I, J, L, X,and Z are possible to 
construct on the system. The proposed system consists of a group of 
4x2 track and gear modules [2]. These modules have "extra- 
modular" interlock positions which can be assembled to produce a 
given shape. The braided product of this system differs from the 
product of either the 4-step or 2-step process due to the yarn 
configuration patterns [2], For instance, the yarns in one exterior 
layer of the interlocked structure need not pass through the adjacent 
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layer. Also, in contrast to the 2-step process, no axial yarns are 
needed for structural integrity [2J. 


CARRIER DESIGN 

All 3-D braiding machines need sufficient yarn supplies 
delivered from movable housings, called carriers,, which follow 
predetermined paths forming viable braided structures. Since the 
notion of moving a large yarn supply is problematical, the 
manufacture of small, compact and efficient carrier assemblies is 
crucial. Small carrier assemblies transit machine beds quickly and 
precisely. Yet, on flat, planar machine beds, due to their change of 
distance at the convergence point of the braid, these assemblies must 
retract yarn to prevent yarn slackening [13]. Thus the mechanism of 
the carrier serves to maintain a properly tensioned yarn, to allow 
feeding or retraction of the yarns (if necessary), and to direct the 
movement of yarn paths during the braiding process. 

In 1987, William Heine of ARC filed a patent on an apparatus 

equipped with a main fiber spool which sat atop the fiber carrier 

[14]. The yarn follows a path from the main supply around a guide 
and through an aperture positioned within a take-up spool. The 
rotational axis of the main spool is perpendicular to the rotational 

axis of the take-up spool. This take-up spool is equipped with a 
tensioning device in the form of a spring assembly [14]. The carrier 
can be equipped with wheels for ease of movement. 

Fulfilling the need for a carrier assembly that is compact, 
reliable and easily constructed, the Heine carrier also employs an 
adequate tensioning mechanism that has horizontal orientation above 
the yarn supply. This mechanism consists of a ratchet within the 
take-up spool that releases a pretensioned spring giving rotational 

force for the desired fiber tension [14]. While the braiding 
progresses, the take-up spool rewinds slack in the yarn preventing 
sagging. Also, constructed of a suitable material like metal or plastic, 
the carrier is economical. 

The major claims of the patent provided a basis for future 
research and invention in the area of carrier design (Figure 10) [14]: 

(1) an elongated housing with an open end 

(2) a cap member mounted on the open end of the housing which is 
removable and 

(3) a take-up spool rotatably held on the cap member with 
perpendicular rotation to the axis of rotation of the supply spool. 


The design is suitable for use on a multi-ply braiding machine which 
uses a matrix array of carriers undergoing alternate track and 
column shifts. 

At the 3-D braiding laboratory of North Carolina State 
University, College of Textiles, El-Shiekh and his associates devised a 
suitable carrier device for the 2-step braiding process [20]. For the 
2-step process, which requires an expansive yarn supply, the group 
designed and built a prototype having a large rewinding length. The 
prototype meets the requirements for successful carrier design by 
having continuous yarn supply and adequate tension level [18]. 

The NCSU carrier assembly had several notable features. First 
it is driven by a mechanical motor which feeds yarn, rewinds yarn 
and tensions the yarn (Figure 11) [20]. The action of the compact 
carrier is economical, furnishing a continual yarn supply with 
predetermined tension. The carrier made of steel has a rewinding 

length of fifteen feet [18]. The reliability of the carrier action 
compounded by its simplicity enhances the marketability of this 
prototype. 

In 1990 Brown of ARC patented a fiber spool apparatus having 
a rewinding capability of smaller proportions as compared to the 
NCSU carrier assembly. Equipped with a motor housing, which 
consists of a coiled spring, axle and supply spool, the compact device 
advances yarn to the braiding apparatus (Figure 12) [5], The carrier 
also has a tensioning mechanism similar to the NCSU design operation 
to adjust and maintain suitable levels of strand tension during 
braiding. 

A fully automated braider equipped with well-tensioned 
carriers for the 2-step process was suggested by Du and his 

associates at the University of Delaware. The braider and its carriers 
not only fulfill the process requirements of effective speed, but also 
are cost effective. The system comprises "motorized" carriers 
mounted on a reformable track [11]. The system regulates carrier 

position at the start of each cycle (2-Step) due to the varying 

distance moved by each carrier. Each carrier assembly has a "fixed 
parking station" [11]. The carriers move from these stations 

simultaneously traversing in two steps and stopping at their next 
respective station. 

The motorized carriers have the following components: (1) a 
small DC motor with gear head, (2) driving traction, (3) a bobbin 
holder permitting yarn retraction by an electronic tensioning device, 
(4) power contact - brushes and (5) an off/on control microswitch 
[11]. Thus, the Du motorized carrier system utilizes electronic and 
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mechanical devices that insure the precise carrier movement needed 
for the 2-Step process. 

Yokoyama and his colleagues at Kyoto Institute of Technology, 
Japan, devised a computer-controlled system for driving the motion 
of carriers on a 3-D braiding machine. The system facilitates the 
change of the carrier track configuration [24]. Also, the system creates 
and modifies the carrier track for certain braid specifications. 

This robot driven system is a noteworthy advancement. A 

robot answers optical sensors that cause the system to follow a 
guide-tape track plan [24], A computer program dictates the 

relationship between carrier motion and yarn path. All processing 
parameters for carrier motion are programmed through the 

computer. The automated (self-driven) simulation controls both the 
carrier speeds and track positions for the construction of particular 
braids [24]'. Also the computer simulates the product being formed 
by the carrier movement while activating the self-driven system. 

The aforementioned ATI cylindrical braiding apparatus 
advanced not only the design of the machine bed, but served to reduce 
the complexity of carrier design. Its inventors realized that it is 

"desirable to minimize the difference between distances from a 
carrier member to the consolidation point..." [15]. The carrier 
members consist of T-shaped, platform bases for ease of sliding. The 
platform of the base has a spindle with a yarn supply. This platform 
includes a guide tube that carries yarn from a supply spool. No 

tensioning and retraction mechanism is needed due to the side by 
side carrier arrangement on the cylindrical machine bed. Next, for 
the manufacture of an axially reinforced braid, Spain and Bailey also 
designed carriers including "free-wheeling pulley assemblies [15]. 
These assemblies mounted on support rods include a yarn supply 
and connect to a bias movable pulley along the guide rod axis. 


BEAT-UP MECHANISMS 

A necessary element for the production of 3-D 

braided preforms is the convergence of the structure to prevent fiber 
entanglement and to produce a uniform structure. With cartesian, or 
4-step braided structures, the convergence of the structure is 
obtained by the use of a beat-up or combing action. In most cases 
this beat-up action is performed manually, using rods, or by hand. 
Manual beat-up is not only labor intensive, but is slow and 

manufactures inconsistent products. In an effort to improve the 
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quality of the preforms, researchers are automating the beat-up 
process. 

One of the earliest beat-up mechanisms is found in Florentine's 
MAGNAWEAVE process. In MAGNAWEAVE braiding, a combing 
action follows the intertwining of the yarns, controlling the 
orientation and density of the material [12]. 

In 1990, Ashton and Patterson of ARC developed a beat-up 
mechanism for multi-ply braiders. Monofilaments are mounted on 
one side of the braider with a rigid connection at one end and a 
flexible connection (such as a spring) at the other (Figure 13) [1]. 
Pusher arms take the lines to a central point between the fibers. At 
this point, gripper arms come in and grip or engage the lines, pulling 
them through to the other side of the braider. The gripper arms 
travel down support posts, bringing the interlacings down to the 
point of convergence [1], The motion of the pusher and gripper arms 
can be controlled by hydraulic piston and cylinder devices. 

At NCSU, El-Shiekh and his graduate students have developed a 
robotic arm to replace the beat-up motion usually performed 
manually. Utilizing a computer, a controller, two stepping motors, and 
a series of guide rods, a mechanical arm is used to bring the 
interlacings of fibers to the convergence point of the structure [18]. 
The beat up arm, after entering the area between carrier columns, 
swings up toward the convergence point while the beat-up unit 
travels parallel to the braider bed. This beat-up action is repeated 
for a designated number of positions along the braider bed. By 
positioning the beat-up arm between the carriers and moving the 
arm up and across the braider bed simultaneously, the contact 
between the fibers and the beat-up arm is minimized, thus reducing 
fiber damage due to the beat-up process [18]. In addition, the 
mechanized beat-up action produces a repeatable, uniform structure. 
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ABSTRACT 

The development of prototype braiding machines for the production of generalized braid 
patterns is described. Mechanical operating principles and control strategies are presented loi two 
prototype machines which have been fabricated and evaluated. Both machines represent advances 
over current fabrication techniques for composite materials by enabling nearly ideal control of fiber 
orientations within preform structures. They permit optimum design of parts that might be subjected 
to complex loads or that have complex forms. Further, they overcome both the lack of general 
control of produced fiber architectures and the complexity ol other weaving processes that have 
been proposed for the same purpose. One prototype, die Farley braider, consists ol an array ol 
turntables that can be made to oscillate in 90 degree steps. Yarn ends are transported about the 
surface formed by the turntables by motorized tractors which arc conti oiled thiough an optical link 
with the turntables and powered through electrical contact with the turntables. The necessary 
relative motions are produced by a series ol linear tractor moves combined with a series ol turntable 
rotations. As the tractors move about, they weave the yarn ends into the desired pattern. The second 
device, the shuttle plate braider, consists of a braiding surface lormed by an array of stationary 
square sections, each separated from its neighbors by a gap. A plate beneath this surlace is caused 
to reciprocate in two perpendicular directions, first in one direction and then in the other. This 
movement is made possible by openings in the plate that clear short columns supporting the sui lace 
segments. Yarn ends are moved about the surface and interwoven by shuttles which engage the 
reciprocating plate as needed to yield the desired movements. Power and control signals are 
transmitted to the shuttles through electrical contact with the braiding surface. The shuttle plate is 
a passively driven prime mover that supplies the power to move all shuttles and the shuttles are 
simple devices that employ only a solenoid to engage the shuttle plate on command. Each shuttle 
is assigned a unique identity and is controlled independently. When compared to each other, the 
Farley braider is felt to have the advantage of speed and the shuttle plate braider, the advantage of 
simplicity. 


INTRODUCTION. 

The work reported here was part of a study of the potential for practical 3-D braiding machines 
having the ability to generate any prescribed relative motions of the braidingyarns. 1 1-8] Such a braider 
would afford complete freedom in defining braid structures and would facilitate the design of high- 
performance composite components. An ideal braider would possess only the mechanical com- 
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plcxity needed to control the braiding pattern, yet be capable of producing generally variable 
patterns. Most 3-D braiding schemes either achieve simplicity by limiting flexibility or seek 
flexibility at the expense of complexity. For example, most braiders yield structures having 
characteristics inherently linked to the process and that cannot be changed, and they therefore have 
no flexibility at all. Examples are: conventional 2-D mechanical braiders [2,6,7], the Florentine 
Magnaweave scheme [3,4], and the two-step braider [2,6,7], all of which produce braid patterns 
that are intrinsic to the process. On the other hand, methods such as the AYPEX [8] procedure 
and the Bluck [1] and Fakuta [5] processes possess the necessary flexibility but suffer from 
complexity in their implementation. This complexity becomes overwhelming when the processes 
are scaled up to produce large sections with full flexibility. Even when the size of the product is 
modest, the flexibility required to produce a variety of structures requires a great deal of redundant 
capacity. 

Two different approaches were examined in detail as part of this study and small prototypes 
were built. The first machine was based on an idea developed prior to the start of the study by one 
of the authors (Farley) and is referred to here as the Farley braider. The other was originated during 
the course of the study and is referred to as the shuttle plate braider. Four general requirements 
were set initially and were principal influences on the detailed design of the prototypes. These were 
as follows: 

• A completely general braiding capability was to be attained which would permit any 
paiticulai yarn end to be moved from any position on the braiding surface to any other 
position by any prescribed path. 

• The mechanical construction and control requirements had to be practically feasible in 
machines of large size. 

• A large number of non-braiding, axial yarns were to be accommodated. 

• The physical dimensions ol the braiding surface were to be minimized, ideally no greater 
than needed to allow the use of yarn packages of one inch diameter. 

THE FARLEY BRAIDER 

The braiding surface of the Farley braider consists of an array of turntables that can be made 
to oscillate in 90° steps. Yarn ends are transported about the surface by motorized tractors which 
are controlled through an optical link with the turntables and powered through electrical contact 
with the turntables. The tractors are guided by track segments mounted on the turntables and are 
propelled by small motors driving through a gear train against a rack incorporated into the track 
segments. The necessary relative motions are produced by a series of linear tractor moves combined 
with a scries of tin ntablc rotations. With the turntables oriented in one direction, the track segments 
become aligned and form a continuous path in one direction for the tractors. When oriented in the 
other direction (rotated 90°), the track segments again form a continuous path but in a direction 
perpendiculai to the first orientation. A proper combination of turntable rotations accompanied by 
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appropriate tractor moves along the paths formed by the track segments weaves the yarn ends into 
the desired pattern. Stationary axial yarns, if needed, would pass through the braider surface in 
the space between the turntables. The basic concept is illustrated in Figure I. 



Figure 1-- Prototype Farley Braider Showing Basic Components of the 
System 

As originally conceived, the Farley braider would consist of a large array of independently 
controlled rotating turntables. Such independent control of the turntables is desirable but would 
require an immense number of actively controlled devices when implemented on a piactical scale. 
Consequently, in the prototypes, and likely in a practical implementation, the turntables were made 
to operate in unison and were actuated by a single actively controlled prime mover. The switching 
action of the turntable array is controlled by a computer with each rotation occurring after a complete 
set of tractor moves along a given axis. For example, with the turntables set in the X-axis, the 
tractors are moved as necessary in the X-dircclion. Alter each tractor reaches its cunent destination, 
the turntables are switched to the Y -axis. The next set of moves of the tractors, all in the Y-dircction, 
then take place. The turntables are then returned to the X-axis orientation, and another set ol tiactoi 
moves occurs. The switching back and forth of the turntables continues in this alternating manner 
until the entire braiding program has been executed. 

Each of the tractors incorporates an electronic control circuit and a small d.c. motor and gear 
train. Power is conveyed to the tractors through contact with electrically isolated conductors 
incorporated into the turntables. Control signals arc transmitted by frequency modulated optical 
signals through emitter-detector pairs mounted in the turntables and in the tractors. Signals are 
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directed to specific turntables in the array. Directional start signals are transmitted to the locations 
occupied by the tractors. Stop signals are erected at the destinations of the tractors. The tractors, 
once set in motion, continue in motion until they encounter the stop signals. When all tractors have 
completed moving, the turntables arc commanded to rotate a quarter-turn to align with the opposite 
coordinate axis. In the current prototype, this rotation is accomplished via solenoid controlled 
valves and pneumatic cylinders. After the rotation has occurred, the next set of instructions is sent 
to the tractors and the appropriate moves are accomplished as before. The sequence of operations 
continues, alternating between tractor moves and turntable rotations, until the desired braided shape 
is completed. 

THE SHUTTLE PLATE BRAIDER 

The shuttle plate braider consists of a braiding surface formed by an array of stationary square 
sections, each separated from its neighbors by a gap (Figure 2). A fiat plate beneath this surface is 
caused to reciprocate in two perpendicular directions, first in on e dire ction and then in the other. 
This movement is made possible by openings in the plate that clear short columns supporting the 
surface segments. Yam ends are moved about the surface and interwoven by shuttles which engage 
the reciprocating plate as needed to yield the desired movements. An operating sequence is 
illustrated schematically in Figure 3. 



Figure 2-- Prototype Shuttle Plate Braider 


In the current version, both power and control signals are transmitted to the shuttles through 
electrical contact with the braiding surface. The shuttle plate is a passively driven prime mover that 
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supplies the power to move all shuttles. The shuttles themselves are very simple devices that employ 



Direction of Shuttle Elate Movement 



Figure 3-- Schematic Representation of Shuttle Plate Operating Sequence. 

only solenoids to extend a plunger that engages the shuttle plate on command. Each shuttle is 
assigned a unique identity and is controlled independently by directing commands to particular 
addresses. 
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The shuttle plate approach was originated after identifying two additional attributes likely 
to be found in an ideal braiding machine. These attributes were as follows: 

• The braiding action should require the minimum of actively and independently controlled 
devices required be reduced to an absolute minimum. Assuming simultaneous and 
independent control of all braiding yarns, this minimum would equal the number of 
braiding yarn ends. 

* Actively controlled actions should be mechanically uncomplicated. 

The shuttle plate device possesses both of these attributes. The move commands are transmitted 
directly to the shuttles and the controlled action is a simple on/off command to actuate a solenoid. I 

Such simplicity is in stark contrast to other methods that require control of actuators, direction ! 

control devices, and the like at each point on a braiding surface that could be occupied by a yarn | 

end. For example, the Bluck f 11 and Fakula [5] braiders and the AYPEX [8] process, as originally 
proposed, require an x-y grid of actuators, all independently controlled and quite complicated in 
their function. A 1 00*100 braiding grid would require ten thousand such actuators, even when only 
a few hundred or perhaps a few dozen yarn ends are being conUolIed. With the shuttle plate 
approach, the size of the braiding grid has no effect on the number of required controlled devices. 

It is possible to make the shuttle plate itself a completely passive device by driving it j 

alternately in one direction, then the other at a constant frequency. However, the braiding process 
can be sped up by independently driving the plate in the two axes in a fashion Uncontrolled eliminate 
wasted moves when possible. Such control adds one clement to the number of controlled devices 
and promises substantial speed increases for certain braid patterns. 

THE PROTOTYPE MACHINES 

Both braiding approaches have been reduced to practice in the form of small devices consisting 
of a 5*5 braiding grid with three shuttles. Both work well. However, the greater mechanical 
complexity of the Farley braider rendered it much more temperamental and difficult to make reliable. 

The shuttle plate braider works with hardly a hitch. The prototypes arc shown in the photographs 
of Figures 4-6. 

COMPARISON OF THE TWO APPROACHES 

The two braiders discussed both accomplish generalized braiding, both in theory and as 
reduced to practice, in that they arc both capable of moving any yarn end from any endpoint to any 
other endpoint by any path specified by the programmer. To the investigators’ knowledge, this has 
not been practically achieved before. 
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Figure 4-- Photograph of Prototype Farley Braider Showing Two Yarn 
Carriers and with Track Segments Oriented along the X-Axis 



Figure 5— Photograph of Prototype Farley Braider with Track Segments 
Rotated to Align with the Y-Axis 
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Figure 6 Photographs Showing Prototype shuttle plate Braider 
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Comparing the two braiders against each other, as opposed to comparing against other 
braiding techniques, the following advantages and disadvantages have been determined through 
operation of the two prototypes in the laboratory. The shuttle plate braider is a very simple design 


from a mechanical viewpoint, and its control requirements are as simple as they can be made, since 
all that is required are simple on/off commands. Further, all the power needed to move the shuttles 
is derived from the shuttle plate, and thus little power is needed for the shuttles themselves. The 
modified Farley braider does not have this simplicity, but it docs have the advantage of speed for 
braiding patterns which require numerous long moves of the yarn carriers. In addition, while at any 
given lime all the yarn carriers of the modified Farley braider must move along a given axis, some 
can be moving in the forward dir ection while others are moving in the reverse direction. Of course, 
this speed advantage diminishes as the average move length of a yarn end becomes shorter in 
complex patterns. 

For the Farley braider the timing and synchronization of moves between yarn carriers become 
an issue of concern, especially as the number of carriers increases. This concern could force the 
use of more complicated devices, such as stepper motors, and neighbor proximity detectors. The 
shuttle plate braider does not have this timing difficulty, since all shuttle moves are automatically 
synchronized by the driving plate. 

Although both braiders transmit power to the yarn carriers via the braiding surface, the need 
for such power is significantly different in kind. The shuttle plate braider needs power on the surface 
to engage the solenoid in each shuttle. As currently implemented, this power is held continually to 
maintain engagement. If several solenoids are activated at the same time, this would require high 
currents on the surface. However, there are numerous ways to overcome this difficulty in a scaled 
up version of the shuttle plate braider. These include such options as using mechanical latching and 
momentary currents to engage the latch. For the Farley braider, the motors must be powered 
continually. Thus the current, of necessity, must increase as the number of moving yarn carriers 
increases. There is no simple solution to this dilemma. Finally, as the size of the braiders is scaled 
up to practical applications, addition difference would be evident. The shuttle plate braider scales 
up readily, since the control problem remains the same no matter the size of the braider. The Far ley 
braider might more easily be implemented on an upwardly curved surface. Use of such a surface 
would reduce the size of the braiding surface needed to control braid angles. However, such an 
approach would complicate the design significantly. For example, the turntables of the modified 
Farley braider would have to be of unequal size or rotate through unequal angles, depending upon 
location on the curved surface. Finally, set up and operation of the shuttle plate braider is much 
easier and more reliable, as discovered in operations to date. 

SUMMARY 

A successful attempt to develop and implement generalized, three dimensional braiding has 
been accomplished and two practical schemes for implementation have been designed, built, and 
tested. Both schemes, as implemented, produce the motions necessary to generate any desired 
braiding pattern. Each scheme has its advantages and disadvantages. However, the shuttle plate 
braider offers the greater immediate promise because of its mechanical simplicity and ease of 
control, especially when scaled up to practical dimensions. 
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Over the past few years, Multiaxial Warp Knit (MWK) 
fabrics have made significant Inroads into the 
industrial composites arena. This paper will examine 
the use of MWK fabrics in industrial composite appli- 
cations (1-10). Although the focus will be on current 
applications of MWK fabrics in composites, this paper 
will also discuss the physical properties, advantages 
and disadvantages of MWK fabrics. The author will also 
offer poss ibilities for the future of MWK fabrics in the 
industrial composites arena. 


According to a Frost and Sullivan study, the demand for high 
performance composite materials In the USA should more than triple 
between 1987 and 1992 (3). A yearly growth rate of high performance 
composite materials for this same period has been forecast to be at 
least 25%. Others predict that the worldwide market for composite 
materials, currently at $2 billion, will grow to $20 billion by the 
year 2000 ( 8 ) . One of the biggest reasons for this projected growth 
is the ever-expanding use of composite structures in industrial 
applications . 

Advanced composite structures, high strength materials combined 
with a resin system, are being used in the most diverse sectors of 
the industrial arena. One fabric that is playing an increasingly 


77 



important role in this arena is the Multiaxial Warp Knit (MWK) . MWK 
fabrics are ideally suited for this type of end use because of their 
flexibility and engineerabil i ty . MWK fabrics, produced in a one-step 
process , have properties similar to those of quasi - isotropic layup 
structures. Besides good handleabili ty , MWK fabric preforms are 
extremely conformable. 

This paper will examine the use of MWK fabrics in industrial 
composite applications. The focus will be on the MWK fabric itself, 
and why it is being used as a preform in industrial composites. 
Current end use applications for the MWK fabric composites will be 
discussed and poss ibili t res for the future of MWK fabrics in the 
industrial composites arena will also be offered. 

What Is An MWK Fabric Composite 

The basic principle of fabric composites is to combine the singular 
properties of more than one substrate and/or medium to create a 
single structure that performs much better than any of the individual 
components. An advanced composite is based on the same principle as 
using steel reinforced concrete to build a bridge, or straw and mud 
to make a brick for building a fireplace. In advanced composite 
structures, reinforcing fibers or fabrics are embedded in resin 
systems such as epoxy or other binders. Multiaxial warp knitting is 
a means of manufacturing stable reinforcing fabrics from fibers such 
as graphite, glass, Kevlar, ceramic, or other textile fibers which 
have the ability to resist stress and strain from any direction. 

MWK fabrics are unique structures which are produced by warp 
knitting techniques. With these techniques, straight ends of 
parallel and uncrimped yarns are inlaid into the knitted structure at 
virtually any desirable angle. Because of the versatility of this 
process, fabric characteristics can be engineered into the structure 
to give the ideal combination of mechanical properties at a favorable 
production cost. This produces MWK fabrics, with the combined 
advantages of design flexibility, performance, productivity and 
availability, and the potential to become a major fabric preform 
for industrial' composites . 

An Industrial composite can be defined as any composite structure 
that is not being used as a household good such as kitchen utensils, 
furniture, clothing, jewelry, and/or sporting goods. 

THE MWK FABRIC 

MWK fabrics generally possess up to four different load bearing yarn 
systems arranged so that each can take on stress and strain in 
virtually all directions. Since these load bearing yarns lie 
straight In the fabric, with no crimp, the physical parameters of the 
individual yarn system are fully utilized. 


From a structural geometry viewpoint (see Figure 1), MWK fabrics 
consist of warp (0 degrees) t weft (90 degrees) , and bias (+/- various 
degrees) yarns that are stitched together during the warp knitting 
process by a fifth yarn system through the thickness of the fabric 
structure (4) . The warp and weft yarns stabilize the fabric in the 
machine and cross machine directions while the diagonally arranged or 
biased yarns absorb tension from any required angle. The fifth yarn 
precisely binds together all of the load bearing yarn systems. The 
bias or diagonal yarns in the fabric can be inlaid at any angle along 
the plane of the machine direction, the most common of which is +/- 
45 degrees. It should be noted that all four load bearing yarn 
systems do not have to be used in a MWK fabric construction. Also, 
different yarn types and counts can be used in each of the yarn 
systems. MWK fabrics allow the designer freedom of choice in the 
arrangement of load bearing yarns in the structure. This provides 
the fabric designer with a tremendous amount of design flexibility. 

Unlike the crimp inserted into the yarn in woven fabrics during 
weaving, the load bearing yarn in MWK fabrics lies straight and 
parallel to other yarns in its yarn system. This characteristic of 
MWK fabrics allows for the yarn properties to be more fully utilized 
in withstanding in-plane forces. Fabric design is made easier 
because the designer can more accurately calculate the tensile load 
of a MWK fabric with a much higher degree of confidence than had been 
previously attainable with woven fabrics. Testing is still required 
to verify the estimates, but the starting points can be calculated 
instead of guessed at. In addition to design flexibility, isotropic 
stress and strain resistance, excellent tear resistance, and improved 
conformability to complex shapes are also characteristics of the MWK 
family of fabrics. 

MWK fabrics are capable of withstanding stresses and strains in 
an optimum fashion. This is due largely to the parallel and straight 
arrangement of the load bearing yarns in the MWK fabric. Because the 
load bearing yarns lie straight in the fabric, their tensile 
properties are fully utilized and are able to absorb tension without 
the elasticity that occurs when the yarns are crimped or in a 
wavelike form, such as in a woven fabric. Because of the arrangement 
of the yarn systems at various angles (see Figure 1), the MWK fabric 
is able to withstand shear forces from various angles 1.7 1 . This 
Isotropic ability is particularly important in many of the advanced 
composite structures and a primary reason for the use of MWK fabrics. 

Due to the parallel nature of the load bearing yarns In the MWK 
fabric, excellent tear propagation resistance is achieved. If a tear 
were to occur in the MWK, the yarn layers would shift slightly under 
the force of the tear and bunch together. The inherent movement of 
the load bearing yarns and the diagonal element would act to 
reinforce the area and prevent further tearing (7 1 . This resistance 
to tear propagation becomes increasingly important when a MWK 
structure Is damaged while In use (such as in a sail, inflatable 


structure, or aircraft skin). The damage is minimized by the 
resistance to further tearing ( 2 ) . 

Because the yarn systems are not interwoven, but rather lie 
directly on top of each other and are held together by the fifth yarn 
system, conf ormabili ty of the fabric is greatly improved. This 
allows the MWK fabric preforms to conform to many complex geometrical 
shapes and still maximize the translation of fiber mechanical 
properties to the composite structure. The conf ormabil i ty of the 
uncured MWK fabric also provides good shape retention during the 
laying up and curing process. 

By combining a web or nonwoven fabric (usually nylon, polyester 
or fiberglass) to the MWK fabric during the knitting process, it is 
possible to control many other physical aspects of the fabric 
structure. Both the MWK fabric and the web can demonstrate their 
specific advantages. Because the web is fed into the knitting 
machine during production of the MWK fabric, it is linked to the load 
bearing yarns by the stitch yarn, rather than being rigidly bonded. 
This allows for a certain amount of give to the MWK and web 
structure. 

Adding a web to the MWK structure allows the designer even 
greater design flexibility. Addition of a web can further control 
the strength and elongation of the MWK fabric, and at the same time 
provide variations of fabric cover and density, water and air 
permeability, stiffness, thickness, initial tear resistance, 
increased tear propagation resistance, and yarn slippage resistance; 
all of which can be tailored by the designer to meet the end use 
requirements of the MWK fabric. The web also provides greater 

stability of the yarn layers during the further processing stages of 
creating an industrial composite. It has been shown that the 
addition of a web to the MWK preform increases the flow of resin 
during the resinating process in making the advanced composite 
structure. This results in better processing times and lower 
production costs ( 5) . 

Few, if any other fabric production techniques, offer such a wide 
range of properties with such versatility In relation to type and 
count as MWK fabrics. Although a comprehensive data base is not yet 
available for MWK fabrics due to their relatively recent introduction 
into the industrial composites arena, several studies have been done 
and are being carried out to assess their potential (6) . The obvious 
properties that MWK fabric composites offer the industrial composites 
arena are their incredible design flexibility, isotropic stress and 
strain resistance, tear propagation resistance, and conformability 

m. 

End Use Applications of Industrial Composites Using MWK Fabrics 

Because the MWK fabrics are relatively new to the industrial 
composites arena, the number of recognizable end use applications is 
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relatively small when compared to those of traditional woven fabrics. 
But when given the head start that woven fabrics have had, some 
several hundred years, the progress made by MWK fabrics in the 
industrial composites arena is nothing less than phenomenal. 

The majority of current end uses for industrial composites made 
from MWK fabrics can be separated into two different industries, 
marine and aerospace. Probably 65 percent of all MWK fabrics 
currently made are used in marine composite applications, while 
another 20 percent are used in the aerospace industry. The remaining 
15 percent encompass all of the varied end use applications being 
evaluated with MWK fabric composites. Table I divides many of the 
current end uses into the three categories. 

MWK fabrics are becoming the fabric preform of choice in the 
marine industry, especially in yachts, sailboats, and high speed 
racing boats. The MWK fabric composite is generally used in these 
vessels for the hulls, deck supers tructures and subs true tures , and 
motor bays . 

Because of the isotropic properties of the MWK fabric structure, 
boat designers are finding that they can use less MWK fabric in the 
composite structure and still maintain, or often improve upon, the 
structural integrity and torsional stiffness of the boat. This also 
means that boat hulls made of MWK fabric composites can withstand 
greater stresses and strains with less overall weight. Less overall 
weight obviously requires less energy to power the boat, which 
translates into fuel savings and/or faster boats. The improved 
structural integrity makes the boat safer at higher speeds . 

MWK fabric composites are being used in most of the fastest ocean 
going racing boats and yachts because of their increased stability 
and weight savings. The hulls and masts of several of the sailboats 
used in the America's Cup competition were made of MWK fabric 
composites because of the performance edge experienced by using these 
composite structures. Partially as a result of using MWK fabric 
composite structures, speed boats and racing boats are achieving 
speeds previously thought to be unreachable with any degree of 
safety. 

Another example of improved performance as a result of using MWK 
composite structures in marine applications was seen in a new 
generation of racing shells used by some of the top rowing teams in 
the country. The shells (long narrow row boats, usually powered by 
8 rowers) were found to give a greater translation of power into 
speed because of the improved torsional stiffness. This allowed 
energy to be translated more directly into speed rather than being 
absorbed by the shell when flexing. 

MWK fabrics are also being looked at for applications in sails , 
In this case, however, the composite structure is the MWK fabric 
combined with a plastic film, usually through laminating. Because of 
the lack of crimp in the yarn in the MWK composite structure, the 
force of the wind is immediately translated into power and not 



absorbed at all by the crimp deformation associated with woven fabric 
structures . 

Fabric composites of all types are being used in the marine 
industry because of their inherent resistance to corrosion. This 
saves on the manufacturing cost because expensive metal treatments 
and repeated paintings are not needed to protect the craft from the 
corrosive nature of salt water. 

In the aerospace industry, which includes the military, aerospace 
and commercial aircraft industries, MWK fabric composites are being 
used more and more. Relatively speaking, they are being used In this 
industry for many of the same reasons as in the marine industry, 
namely reduced weight, and increased strength and integrity. Because 
of the flexibility and tailorability of mechanical and physical 
properties, MWK fabric composites can be customized for the 
application and specific properties can be emphasized to suit the 
particular need. 

Currently, there are very few aircraft that use MWK fabric 
composites in critical structures such as the fuselage or wings. 
Most current applications center around the skin of the aircraft. 
Other areas of use are in the top and side tail units, fuselage 
paneling, leading edges on side rudders, and engine paneling. MWK 
fabric composites are also being evaluated for rotor blades, outer 
skin, and ballistic protection for helicopters. It is thought that 
the use of MWK composites is also being evaluated in the new military 
plane/helicopter, the V-22 Osprey, and the all-composite Beech 
Starship business plane (8) . 

The lower weight achieved through the use of MWK composite 
structures means that less fuel is consumed by the aircraft, which 
translates into significant energy savings for the user. Also, 
because of the improved structural integrity offered by the MWK 
fabric composite, it is believed that safety Is enhanced. 

Other various applications for MWK fabric composites can be found 
in the industrial composites arena. In Europe, MWK fabric composites 
are being used for flooring in sports halls where the combination of 
multi - directional force distribution and excellent tear resistance 
are beneficial. The MWK fabric composite also helps to improve the 
sound damping characteristics of the flooring. MWK fabrics coated 
with rubber are also being used in the Industrial roofing industry. 

The Future of MWK Fabric Composites In Industrial Applic ations 

The main obstacle in gaining acceptance for MWK fabric composites in 
the industrial composites arena has been the lack of confidence 
derived from inexperience and a lack of sustained performance data 
for MWK composites. With the advent of time, these obstacles will 
undoubtedly be overcome and MWK fabric composite usage will 
dramatically increase ( 1) . 
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Numerous end use possibilities exist for MWK fabric composites. 
They can not only be used to replace traditional materials, but also 
to improve the performance of many new industrial composites which 
seem to have reached their performance limit. MWK fabric composites 
can be used to replace traditional structural materials such as 
concrete, wood, and steel, thus creating new possibilities in various 
industries and end uses. 

The advantages of using MWK fabrics in composite structures are 
clearly the flexibility and freedom of choice in the desired 

properties in all directions which can be matched to individual 
needs. As a result, MWK fabric composites are uniquely suited to a 
wide range of industrial applications. 

MWK fabric composites incorporating a nonwoven structure are 
ideally suited for many high strength geotextile applications, where 
isotropic strength, resistance to tear and tear propagation ,; good 
water permeability, low creep, and good fabric/soil interaction are 
required. 

With the flexibility of fiber placement and potentially high 

productivity, MWK fabric composites are ideally suited for many 

structural load bearing applications in the automotive and aerospace 
industries. MWK fabrics, because of their structural makeup, have 
good flexibility which allows them to be formed during molding into 
virtually any desired shape. The through thickness reinforcement 
provided by the stitching process helps to reduce the possibility of 
delamination of layers in the composite structure.. 

Numerous applications for MWK fabric composites also include 

protective helmets and armored protection of vehicles, buildings, and 
people. Various drive belts, V-belts, fan belts, and conveyor belts 
will benefit from the availability of diagonal load bearing yarns in 
the composite structure. Inflatable rafts, cushions, balloons, and 
fuel cells are ideal applications due to the MWK fabric composite s 
isotropic strength and tear resistance. 

There are probably hundreds or even thousands of other areas that 
could benefit via the use of MWK fabric composites. The. flexibility 
of the MWK fabric system provides endless potential end use 
applications. All that is needed is for the designers and engineers 
of the world to dare to improve upon what they already have and open 
their minds to the future. 

Conclusions 

The MWK process offers limitless possibilities for the formation of, 
new fabric preforms for the composites industry. By varying the 
angle of the load bearing yarn systems, and the type and count of 
yarn used, the strength in any fabric direction can be tailored to 
the requirement, rather than the requirement being tailored to the 
fabric available. Because of the complete versatility of the MWK 
fabric, design flexibility is virtually limitless. 


MWK fabric composites offer many never before realized advantages 
to the industrial composites arena, giving the designer flexibility 
that he or she has never before had. The ability to design a 
composite structure with the load bearing systems aligned precisely 
where they are needed provides opportunities to make better, more 
cost efficient structures. 

The use of MWK fabric composites in many of the areas mentioned 
earlier shows that their development is worthwhile. New end uses are 
being developed daily and in many cases MWK fabric composite 
structures are being used to replace expensive, heavier or 
technically inferior constructions produced from other materials. 
Many times these replacements are converted directly into cost 
savings. Also, with each new application, comes more experience and 
an addition to that ever - expanding data base, which so many engineers 
and designers require. All that is needed is for them to open their 
minds and give MWK fabric composites a try. 
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Properties of Textile Grade Ceramic Fibers 
Eric Pudnos 

Dow Corning Corporation 

The availability of textile grade ceramic fibers has sparked great 
interest for applications in composite reinforcement and high temperature 
insulation. This paper summarizes the properties of various small diameter 
textile grade ceramic fibers currently available. Room temperature mechanical 
and electrical properties of the fibers are discussed for three cases: ambient 
conditions, after heat aging in argon, and after heat aging in wet air. ow 
Corning (R) HPZ Ceramic Fiber, a silicon nitride type fiber, is shown to have 
improved retention of mechanical and electrical properties above 1200C. 

Introduction k 

The utility of ceramic fibers for reinforcement of ceramic matrix 
composites (CMC) is well established. The extensive research on a variety of 
composite systems is worldwide in scope and has produced some outstanding high 
temperature composite materials. The high strength and excellent toughness of 
these systems is a direct result of the mechanical properties of the 
reinforcing fibers. Commercially available fibers such as Ceramic Grade (CG; 
Nicalon (R) fiber, Nextel (R) 480, and SCS-6 offer a composite fabricator a 
range of desirable properties to work with. It is also well recognized that 
in many of these systems, the reinforcing fiber is a limitation due to thermal 
decomposition and/or loss of mechanical properties at the elevated 
temperatures that are encountered difring composite fabrication or end-use 
conditions . 

The purpose of this paper is to review the properties of some commercial 
fiber materials and compare these to a new fiber recently introduced by Dow 
Corning. The paper will be broken down into two sections: first, some genera 
information about HPZ fiber, and then some more specific mechanical and 
electrical properties of various textile grade ceramic fibers. This 
discussion will be limited to the class of small diameter, textile grade 
ceramic fibers generally considered to be most desirable for CMC 
reinforcement. The data presented is limited to direct fiber tests and is not 
composite data. The author feels that these comparisons are useful for 
initial evaluation of fiber properties and material selection. We also 
recognize that such comparisons do not predict composite performance. Ongoing 
studies of fiber-matrix compatibility and fiber interfaces are necessary to 
fully determine performance of composites- 

Gp.neral Information about H PZ Ceramic Fiber 

HPZ Ceramic Fiber is made via a polymer precursor route. This unique, 
amorphous silicon nitride-type fiber has been developed as part of a DARPA/Air 
Force contract. The processing chemistry and characterization of this 
material have been the subject of several publications and will not be covered 
here. Silane monomers manufactured by Dow Corning (Figure I) are reacte to 
form a preceramic polymer which is then extruded through a spinnerette at 
temperatures ranging from 100-300C. The resulting uncured polymer fiber 
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(green fiber) is then cured to a fully crosslinked, non-melting form. The 
cured fiber is then pyrolyzed to yield the ceramic fiber. During pyrolysis, 
various flue gases (CO, CH 4 , SiO, etc.) are given off. 

HPZ was officially introduced to the marketplace at the SAMPE show in 
April. It is available in developmental quantities and has a useful balance 
of mechanical and electrical properties. Typical properties are shown in 
Table I. The filament diameter has an oval cross-section that is calculated by 
taking the geometric average of the major and minor axes. As a developmental 
product, properties typically quoted cover a range of values. As you can see 
from Figure II, over the last two years we have been able to almost double the 
average tensile strength of the fiber through a series of process and 
procedural improvements. While no guarantees are made with regard to 
continuing this rate of improvement, suffice it to say that we are constantly 
working on improving fiber properties and that the values achieved today may 
be improved upon in the future. 

Figure III is a typical range of tensile strength values from one 
particular lot of HPZ fiber. Each data point represents an individual 
filament tensile test (recall that HPZ is currently supplied in 500 filament 
tows). The range of values falls between a low of over 200 ksi to a high of 
about 650 ksi with an average tensile strength of 417 ksi. Such a spread of 
values is typical for ceramic fibers and other relatively brittle ceramic 
material . 

As we'll see in the following section, HPZ maintains a large percentage 
of its mechanical and electrical properties at elevated temperatures. It has 
been proven elsewhere that this behavior can be attributed to the stability of 
the Si— N— C chemistry and that the oxygen in the fiber primarily resides near 
the surface of the fiber. A composition depth profile (Figure IV) produced 
with a scanning Auger microscope reveals that at the surface there is roughly 
37X oxygen and 40X carbon, but at a depth of less than 700 nanometers, these 
concentrations have dropped to 5Z and 12% respectively. 

HPZ Ceramic Fiber is available in both lOOg and 250g spools. The 250g 
spool contains nearly 1.5 miles of continuous tow (Figure V). The fiber is 
also available in several weave styles (8HS and Plain weave). Weavers who 
have handled HPZ for us have remarked that the fiber handles at least as well 
as mid-modulus carbon fibers. The fiber is supplied with a PVA sizing that 
can be removed by either a double hot water wash, or a 30 minute/600C/air 
exposure. Both of these methods will effectively remove the sizing without 
any effect on the mechanical or electrical properties. 

As a class, non-oxide fibers such as HPZ have a higher strength than 
oxide fibers. Since non-oxide fibers typically have a lower density, they are 
considerably stronger than oxide fibers on a specific strength basis. Figure 
VI illustrates this nicely. The small circled area shows the range of 
properties currently possible with laboratory prepared HPZ. The excellent 
combination of specific strength and modulus of T-300 carbon fiber (as an 
example of carbon fibers) is shown by its location in the upper right quadrant 
of the chart. Clearly, we would like to improve the properties of ceramic 
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fibers to equal the room temperature properties of the graphite fibers, and 
have the properties retained at very high temperatures in air. 

Room Temperature Prope rties. 

Table II illustrates the mechanical and electrical properties of several 
textile grade ceramic fibers considered in this paper. Structural and 
compositional information is presented in Table III. Of the fibers in this 
list, only Dow Corning (R) HPZ Ceramic Fiber is likely to be new to this 
audience. The reader should note that unless otherwise specified, fiber 
tensile data is for single filament testing at 1 M length, tensile modulus 
values are calculated using machine compliance corrections, and that 
electrical measurements are made at 10 Ghz (x-band) at room temperature. 

Obviously, it is the elevated temperature performance of ceramic fibers 
that is critical for the majority of applications where they are considered 
for use. Although some workers only consider oxidizing conditions to be of 
interest, we felt that inert, non-oxidizing conditions are useful since they 
represent conditions encountered during some composite fabrication processes 
and some use conditions (For example, the interior of a fully dense. CMC may 
contain a very low partial pressure of oxygen). A two hour, 1400C flowing 
argon heat soak became a somewhat standard test condition during early studies 
within our program. This condition was initially chosen because Nicalon fiber 
was almost completely destroyed by such treatment. As shown in Figure VII, 
both Nicalon and Tyranno almost totally lose their reinforcing properties 
after such a thermal treatment. These studies, as well as work by others, 
show that this is a result of massive decomposition that occurs with loss of 
CO and SiO. As is also reported elsewhere, the oxide fibers generally weaken 
somewhat due to grain growth. In the case of HPZ fiber, it has been 
determined that the chemistry and the amorphous structure are very stable at 
1400C. However the strength of the fiber is sensitive to minute levels of 
impurities. Recently we have tested a number of samples that retain fully 
100% strength and maintain >300 ksi after a 2 hour soak at 1400C. No other 
fiber tested has shown such excellent inert atmosphere stability. 


Similar argon agings have been carried out at lower temperatures. 
Generally, these studies show trends that are consistent with studies reported 
in the literature. Figure VIII shows the large weight losses that occur with 
Nicalon and Tyranno as CO evolves at higher temperatures. Note however, that 
the High Volume Resistivity (HVR) Nicalon fiber does remain reasonably stable 
for 65 hours at 1100C. As expected, the more stable chemistry of Nextel fiber 
shows only the weight loss associated with organic sizing present on the tows. 
Similarly, HPZ fiber shows the sizing burnoff and finally the loss of some 
surface oxygen at 1300C. 

The electrical properties of the fibers essentially reflect the changes 
in bulk chemistry occurring. Figures IX and X show that the dielectric 
behavior of Nicalon and Tyranno change significantly as the chemistry changes 
due to loss of CO. In sharp contrast, HPZ and Nextel are electrically very 
stable as a result of the much more stable chemistry of these fibers. 

Since certain applications require that fibers are exposed to air at 
elevated temperatures, it is important to examine heat aging under oxidizing 
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conditions. We selected 100 hours at 1000C as a representative condition. 
Flowing wet air (humidified at RT) was selected for its aggressiveness. In 
the study illustrated in Figure XI, we found that SiC-type fibers CG Nicalon, 
HVR Nicalon, and Tyranno lost strength in sequence as expected based on our 
findings in the inert atmosphere. Thus, the retained strength was (in order) CG 
Nicalon, HVR Nicalon and Tyranno, Once again, the titanium present in Tyranno 
did not aid stability under these conditions. Another surprise was the 
significant strength loss of the oxide fibers under oxidizing conditions. 
Apparently, some grain growth occurs even at 1000C in air causing a strength 
loss. Of the commercial fibers, only CG Nicalon retains >200 ksi strength. 

Although it was expected that a silicon nitride fiber such as HPZ would 
show good oxidative stability, the results were still somewhat surprising. 
Within the normal error of the tensile test, no change in strength was 
observed during this experiment. Other work at Dow Corning at 1100C has also 
shown excellent strength retention for periods up to at least 15 hours. 

As illustrated by the data, HPZ ceramic fiber strength is unaffected by 
a 100 hour oxidative treatment. Multiple batches and repeat runs show no loss 
of strength although some evidence of surface oxidation is observed. At 
higher temperatures, for instance, this surface oxidation can lead to fusion of 
the individual filaments within the tow bundle. Such fusion creates a very 
brittle, weak tow when treated at room temperature. Paul Sawko (NASA Ames) 
has studied fiber tow oxidation over a range of temperatures as illustrated in 
Figure XII. This data was generated by tensile testing fiber tows at room 
temperature after a 2 hour soak in air. This data nicely shows that while 
only HPZ fiber retains strength well in 1200C heat soaks, the strength 
reduction caused by fiber fusion does degrade strength at higher temperatures. 
Thus, in tow form, the usable temperature of HPZ fiber in free air appears to 
be 10Q-200C higher than for other non-oxide materials, but is still limited to 
about 12Q0-1300C. Note that the HPZ sample used in this experiment was from 
the earlier stages of our program. The denier here is only 400, while today's 
fiber has 1000 d. Also the starting fiber strength is about 200 ksi, and 
current values are around 400 ksi. Newer material would start at about 8 
kg/denier and be expected to follow the same trend. 

Although room temperature tests after heat soaks can be useful in judg- 
ing the performance of a fiber, actual properties at various temperatures are 
necessary to effectively design composites for high temperature use. Workers 
at Penn State University have developed equipment and procedures for elevated 
temperature tensile testing of single filaments in air. One uniqueness of 
this technique is the capability for testing at greater than 1400C. Figure 
XIII shows the results of one test series. As reported previously, the oxide 
based fibers are generally lower strength at room temperature and maintain 
strength reasonably well to about 800C. Above 80QC the strength degrades 
rapidly until at 1200C there is sufficient deformation to prevent testing at 
higher temperatures. In contrast, non-oxide fibers are considerably stronger 
at room temperature and generally maintain strength well up to about 1200C in 
fast fracture tests of this type. Above 1200C, Tyranno fiber loses strength 
rapidly, presumably due to its high oxygen content, while only Nicalon and HPZ 
hold useable strength levels at 1400C. Even at 145QC, the HPZ fiber has 
greater than 100 ksi strength. Elevated temperature modulus data (Figure XIV) 
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generally parallels Che same behavior, although here the excellent stability 
of HPZ's microstructure and fiber chemistry is even more apparent as evidenced 
by the superior modulus retention. 

Conclusions 

1) Non-oxide fibers have good room temperature strength and modulus, generally 
much higher than for oxide fibers, 

2) CG Nicalon and HPZ are the only fibers to retain useable tensile strength 
and modulus above 1300C, 

3) HPZ fiber exhibits good retention of mechanical properties after oxidation 
at 1000C for 100 hours. This is the only material showing no strength loss 
under such conditions. 

4) The handling properties of HPZ are similar to mid-modulus carbon fibers. 
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DOW CORNING® HPZ CERAMIC FIBER 

Typical Properties 

i 

• Composition (Wt. %) 

Si 

67% 


N 

28% 


C 

10% 


O 

5% ! 

• Amorphous Structure 



• Diameter 

10-12 microns j 

• Tensile Strength 

420 ksl (2 

9 GPa) ! 

• Tensile Modulus 

28 Ms! (193 QPa) \ 

• Density 

2.4 g/cc 

1 

• DK 

6.1 


• LF 

0.06 


• Heat Capacity 

0.68 J/g C 



Coeff. of Thermal Ex 
(20-1 100°C) 


4*1 ppm/C 


Table I 


HPZ Ceramic Fiber - Typical Properties 


i 


FIBER PROPERTIES 

\ 


Tensile 


Fiber 

Diameter 

Strength 

Modulus 

DK 

LF 1 

CG NICALON 

15 

400 

28-30 

9.5 

1.0 I 

HVR NICALON 

15 

350 

26-28 

6.4 

0.05 i 

i 

LVR NICALON 

15 

425 

26-28 

>20 

>10 1 

C-Coated NICALON 

15 

425 

28-30 

- 

| 

TYRANNO 

8 

450 

22 

7.8 

i.o i 

i 

HPZ 

10-12 

420 

28-32 

6.1 

0.06 1 

NEXTEL 480 

8-10 

220 

30-32 

6.4 

cos 1 


Table II : Comparable Fiber Properties 
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FIBER CHARACTERISTICS 


Fiber Composition Density (g/cc) Structure 


CQ NICALON™ 

Si-C-O 

2.55 

1 

HVR NICALON 

Si-C-O 

2.35 1 

i 

f 

l nanocrystalllne SIC/ 

* amorphous 

LVR NICALON 

Si-C-O 

2.35 \ 

i 

C-Coated NICALON 

Si-C-O 

2.55 


TYRANNO® 

Si-O-C-Ti 

2.40 

amorphous 

HPZ 

Si-N-C-0 

2.40 

amorphous 

NEXTEL® 480 

Al O /SiO _/B_0_ 

2 3 2 2 3 

3.10 

crystalline, mulllte 


Table III 


Fiber Characteristics 
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Figure III : Failure Probabilities 
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Figure IV : Compositional Depth Profile 
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Figure VII : Inert Aging Comparisons 
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Figure VIII: Aging Effect on Fiber Weight 
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Figure IX : Aging Effect on Fiber DK 
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Figure X : Aging Effect on Fiber LF 
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100 Hours / 1000°C in Air 
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Figure XI : Fiber Oxidation Results 
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Figure XII : Effect of Conditioning Temperature on Break Strength 










N93- 14752 


THE OPTIMAL FIBER VOLUME FRACTION AND 
FIBER-MATRIX PROPERTY COMPATIBILITY IN FIBER 
REINFORCED COMPOSITES 


Ning Pan 

Division of Textiles and Clothing 
University of California 
Davis, CA 


Abstract 

Although the question of minimum or critical fiber volume fraction beyond which a 
composite can then be strengthened due to addition of fibers has been dealt with by sev- 
eral investigators for both continuous and short fiber composites, a study of maximum 
or optimal fiber volume fraction at which the composite reaches its highest strength has 
not been reported yet. The present analysis has investigated this issue for short fiber 
case based on the well-known shear lag (the elastic stress transfer) theory as the first 
step. Using the relationships obtained, the minimum spacing between fibers is deter- 
mined upon which the maximum fiber volume fraction can be calculated, depending on 
the fiber packing forms within the composites. The effects on the value of this maximum 
fiber volume fraction due to such factors as fiber and matrix properties, fiber aspect 
ratio and fiber packing forms are discussed. Furthermore, combined with the previous 
analysis on the minimum fiber volume fraction, this maximum fiber volume fraction 
can be used to examine the property compatibility of fiber and matrix in forming a 
composite. This is deemed to be useful for composite design. Finally some examples 
are provided to illustrate the results [1-14]. 


1 INTRODUCTION 

Adding fibers to strengthen materials is a technique which has been used since ancient times. 
It is applied mainly to materials which are much weaker in tension than in compression so 
that by adding fibers into them the superior tensile property of fibers can be fully utilized 
and stronger new materials arc obtained. 

Yet it is understandable that if very few fibers are added to a matrix, the mateiial 
is weakened rather than strengthened. Therefore there must be a minimum critical fiber 
volume fraction V m i„, only exceeding that with which the fiber reinforcing function can be 
realized. There havebeen several studies which addressed thisproblcm and derived thespecific 
values of V min for various cases and under different conditions [1,8,11]. On the other hand, 
however, as the fiber amount in the system is being increased, the tensile strength of the 
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composite will increase to a certain point where, upon further increasing of fiber amount, 
the bonding between the fibers and the matrix will start to deteriorate as the fibers become 
too close to each other. As a result, for a short-fiber composite, the tensile strength of the 
composite will decline due to the bond failure of the system caused by the excessive fibers. 
Therefore there will be a maximum value of fiber volume fraction as well, that being the 
upper limit of fiber amount allowable in the system for reinforcement. 

There have been no reported studies on this issue as far as the present author is aware. 
This may be due mainly to the fact that, in most cases, the fiber amount which can be 
incorporated into a matrix system is limited by the processing technology [4] so that tech- 
nically it may be difficult for the fiber volume fraction to reach this maximum allowable 
value. Hence the maximum fiber volume problem may not be as significant as the minimum 
one for practical applications. Nevertheless, study on this issue is still desirable partly due 
to its theoretical significance, and more importantly, because the investigation of this 
issue as presented in this article cannot only provide the maximum value of fiber volume 
fraction, but also determines the propertycompatibilityof various fiber and matrix materials 
for a composite so as to guide the design procedures in achieving the optimum composite 
strength and full material usage. 

The present study deals with this problem based on a shear strength criterion between 
fibers within the composite. The effects of fiber length and fiber misalignment are also 
investigated. 

2 THE MINIMUM ALLOWABLE SPACING BETWEEN 
FIBERS IN A SHORT FIBER COMPOSITE 

It has been a well known mechanism that when a fiber composite is under a uniaxial 
tension, the axial displacements in the fiber and in the matrix will be different because 
of the differences in tensile properties of these two components. As a result, shear strains 
will be created on all planes parallel to the axes of the fibers. The shear strain and the 
resulting shear stress arc the primary means by which load is transferred to fibers ( for a 
short fiber composite ), or distributed between and supported by the two components of 
composites. It is through this interaction between fibers and matrix that a fiber reinforcing 
function is realized. There have been several theories trying to explain this fiber-matrix 
interaction. The first one was entirely based on the elastic mechanism by Cox [3] in 1952, and 
is now referred to as the shear lag theory, and another similar version was later proposed 
by Rosen [12], Since then, a number of new theories were suggested such as the slip theory 
(7] to account for matrix plasticity at the fiber surface near the fiber ends, applicable to 
well-bonded reinforced metals, and the theory of frictional sliding [7] to reinforced polymers 
and ceramics. However for the present study, the model of the elastic stress transfer will be 
used as the main theoretical basis. It will be shown that, although this theory basically 
only explains the behaviorof composites at low stress, it still provides adequate information 
in determining the maximum fiber volume fraction for design purposes. Furthermore the 
analysis will surely be helpful in the attempt to look into the case of the inelastic 
interaction as well. 
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Assumptions made in this analysis include : 

1. Since the elastic model is used here, conclusions from the present analysis are valid 
only if the original assumptions associated with this model hold. 

2. The composite consists of many short fibers each with constant length /, circular 
cross-section area Af of uniform radius r and tensile modulus Ej. 

3. All fibers are distributed uniformly along the length of the composite so that the fiber 
area fractions on all the cross sections of the composite are identical. 

4. Both fibers and matrix behave elastically, and the interface transfers the stress between 
fibers and matrix without yielding or slipping. 

5. Fiber ends are all normal so that the shape effect of fiber end on the stress transfei 
[6] is excluded in this analysis. 

6. Furthermore, the fiber-fiber interaction within the composite and the effect of matrix 
property change as a result of the fiber interfering with dislocation motion in the 
matrix are also ignored. 


We take the mean fiber center to center spacing normal to their length to be 2R (see Figure 
1). Assume the composite as a whole is subject to a strain e c which will cause a strain ey 
in a fiber. If P is the load in the fiber at a distance x from the fiber end, then according to 
Cox [3], the distribution of tensile stress in this arbitrary fiber is 


Af 




( 1 ) 


where 

fl - 1 /Cm f 2 I ) 

P r\) E f K HR/r) 

and G m is the shear modulus of the matrix. Note that a = 0 at x = 0, and l. 
The maximum stress occurs at the middle of position x — 1/2 where 


( 2 ) 


&max — 


cosll/? i 


( 3 ) 


It can be seen from Equation 3 that, in order to fully make use of the tensile strength of 
the fibers, i.e. to make o max = a b j, the fiber- fiber spacing R is the key factor for given fiber 
strain, and fiber and matrix properties. 

If r is the shear stress in the direction of the fiber axis, on planes parallel to this axis, 
then at the fiber surface we have 

— — = —2TTTT ( 4 ) 

ax 

Equation 1 and 4 give the expression for the shear stress distribution 


r 


E f e f 1 


2Ef \n(R/r) 


siiih/?(^ — x) 
cosh /? ^ 


( 3 ) 



The maximum value of r occurs at the fiber ends, i.e. at x ~ 0 and /, 


Tmaz ~ Efe/ y 2 E f HR/r) tanh/? I (6) 

and it is zero at the middle of the fiber. Both of these stress distributions are shown in 
Figure 2. The ratio of the maximum value of shear stress to the maximum tensile stress in 
the fiber is 

T rnax I Grn 

max \j 2Ef\n(R/r) C ° 1 4 

This ratio is of great importance as it represents fiber and matrix properties as well as fiber 
spacing 72 within the matrix. The validity of equations 1-7 has been verified by several 
experimental studies [G,7 ]. 

In reinforcing the composite to its maximum tensile strength, the tensile strength of 
fibers has to be utilized to the fullest. In other words, a stress equal to the tensile breaking 
stress of the fibers must be reached at the middle of a fiber, i.e. u max = So the 
above equation can be rearranged as 


Tmax — &bf 1 


coth p- 
2Ef\n(R/r) 4 


This equation gives the relationship between fiber spacing R ( or the spacing ratio R/r) 
and the maximum value of shear stress. When the spacing between fibers R decreases, 


the value of r wax 
determined when 


will increase as shown by Figure 3. The minimum spacing R min is thus 
r max has been increased to such a large value that it reaches the shear 


strength of the matrix adjacent to the interface or the shear strength of the fiber/matrix 
interface, whichever is less, designated as r,. Because of the clastic assumption where the 
matrix can not deform plastically, this will cause either the fiber /matrix interface or the 
matrix to fail in shear. 

Furthermore replacing r nuix by r s and rearranging Equation 8 give the final relationship 
between the minimum spacing ratio R m in/ r anc l the strength ratio the fiber aspect 

J y ~ f 3 

ratio - as well as the modulus ratio 


In (Rmin/r) = (— coth/?b 2 ^r- = — cothfif . / ] 2 (9) 

r, 4 2 E f r s ir]j E f ln{R mijl /r)‘ 2 E f 

This is a transcendental equation for and its solution can only be calculated 

numerically. 

However if the fiber length is relatively long so that coth — » 1, we have an explicit 
relationship between the fiber spacing ratio and the fiber-matrix properties 

HR M = ( 10 ) 

r . 9 


(ZL)* Hfli. 

Rmin/r = e " 



In this analysis, the effect of stress transfer across the fiber ends is neglected which will 
cause an extra load on both the fiber and the matrix in this region. However this effect is 
considered insignificant [11] as long as the fiber aspect ratio Ijr > 10. Also the influence of 
stress concentration across the fiber ends, which will lead to a greater shear stress [8] and 
will affect the slip behavior of the fiber ends, is ignored. 

In addition, in the present analysis, fiber and the matrix are assumed to be completely 
elastic. This is of course an ideal case, and only valid in practice to brittle materials. For 
some cases where plastic deformation of the matrix does exist, the conclusion drawn from 
this study will be a conservative one and certain modification may be needed, since the 
plastic deformation of the matrix will alleviate the shear stress. However a different model 
of the spacing/stress relationship is desirable for a matrix which is significantly plastic and 
flows under loading, or for structures where the effect of frictional sliding between fiber and 
matrix during the stress transfer is not negligible. 


3 THE MAXIMUM FIBER VOLUME FRACTION IN COM- 
POSITES 


As indicated above since there is a minimum spacing R m in/r between fibers within 
a composite below which the structure will collapse due to shear failure, correspondingly 
this minimum spacing will define an upper limit of fiber amount which is allowed to be 
incorporated into a given matrix. The composite will reach its highest strength at this 
maximum fiber volume fraction V max , as there will be a maximum amount of fibers in the 
composite and each is fully utilized. In other words, this maximum fiber volume fraction 
is also the optimal value for maximizing the composite tensile strength. Obviously the 
specific value of V max is dependent on the forms of fiber arrangement within a composite 
as well. Moreover, for the short fiber case where fiber ends don’t meet, the maximum fiber 
volume fraction also varies with the distance between fiber ends. Let us assume this distance 
between the ends of two fibers is 25 f as shown in Figure 4(a). 

The following are the two cases most often encountered in a composite [11]. 


3.1 Hexagonally Packed Fibers 


The fiber arrangement of this type is schematically shown in Figure l-(c). Suppose there are 
totally N fibers within the composite. According to the definition of fiber volume fraction 
of a composite, we have 


Vf = 


F/jfccr 

Vtotal 


(12) 


Considering the hexagonal area enclosed by the dotted line in Figure l-(c), the maximum 
fiber volume fraction in this case is 


Vjj n h ~ 


37r r 2 / 


7T 


) 2 (: 


a#(2ilm, •„)*(/ + 26 f) 2n/3 Ryix i n 1 + 28 j /l 


77) 


(13) 
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When all fibers are packed so densely that they are actually contacting each other in full 
with fiber ends also connected, that is 


Rmin — T j — 0 

the maximum fiber volume fraction becomes 


(14) 


(15) 


Also for the case when the fiber length is so long as / >> 8f that the fiber end effect can 
be neglected 


Vfmh — 


2\/3 R 


(TT ^) 2 


(16) 


Because of this direct relationship between the maximum fiber volume fraction and the 
minimum fiber spacing, it is equivalent in the later analysis to refer to either of them. 


3.2 Square-Packed Fibers 


The fiber arrangement in this case is shown in Figure l-(d) and we will have 


v, = 111 = n (- r *2, 1 

frn 3 (2R min ) 2 (l + 2S f ) VR m J 4 + 2 Sf/l 

For the longer fiber case we have 

* Jtmin 


) 


(17) 


(18) 


In the extreme case when fibers are closely packed to each other so that i2„ ll7l = r, there 
will be 

Vjms = J (19) 

In either of two packing forms, the value of maximum fiber volume fraction monotonically 
increases as the fiber spacing decreases. The relationship of or difference between the 
maximum fiber volume fractions of these two packing forms is given by 


Vjms _ V^3 
Vfmh 2 


( 20 ) 


That is, the maximum possible fiber volume fraction for square-packed fibers is less than 
that of a hcxagonally Packed case. Again because of this direct relationship between the 
two fiber packing forms, for briefness, only the Square-Packed form is used in the following 
analysis. 

Note that when there is fiber misalignment existing in the composite, the fiber arrange- 
ment may not be as regular as the two examples shown here. Consequently the value of the 
actual maximum fiber volume fraction may be lower than the present results. 



4 THE MINIMUM FIBER VOLUME FRACTION IN COM- 
POSITES 


In the next paragraph, e&/, fbm a-nd e ytn represent the fiber breaking strain, the matrix 
breaking strain and the matrix yield strain. Although there may be three cases [11] 

1 * ^ 6 / ^ £ymi 

2. €ym ^ ^6/ ^ ^tmi 

3 . €bm ^ ^/>/* 

existing in composite, each of which will lead to different failure behavior of the composite, 
usually the breaking strain of the fiber is less than the yield strain of the matrix € ym so 
that only the first case is considered here. The treatments of two other cases, however, 
are in principle the same. 

If the variations of fiber tensile properties are ignored, according to the Law of Mixture, 
the breaking strength of the composite a^ c is of the contributions from both fibers and the 
matrix, and can be expressed as [1] 

&bc = 077.1(1 - Vf) + wioVbjVf, Vf > V lllin ( 21 ) 

where a^j is the breaking strength of the fiber, and a j ln is the stress on the matrix at the 
breaking tensile strain of the fiber. The factors i]i and i)o account for the effects of limited 
fiber length for the short fiber case, and of fiber misalignment, and are often called the 
length efficiency and fiber orientation efficiency factors. V„ un is the minimum value of the 
fiber volume fraction which must be exceeded if the strength of the composite is to be given 
by the Law of Mixture. The value of Km',, can be determined analytically, according to 
Kelly [8], as shown below. 

If the amount of fiber added into the composite is very small, it will actually weaken 
the composite so that the strength of the composite becomes 

^6c = - Vf) (22) 

where <7{, m is the breaking strength of the matrix. 

Inserting this relation into Equation 21 gives: 

cq,m(l - Vf) = °7m( 1 - Vf) + vrnovbfVf ( 23 ) 

The minimum fiber volume fraction can be derived from this equation, i.e. 

r„„„ = "‘T "'Z — (24) 

VlVOVbf d* G 6 m & fm 


For a continuous fiber composite where all fibers arc aligned in the loading direction, 


there will be 


VI ~ l ’Vo = 1 
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and 


1 /- &bm G fm 

*min — " 

&bf » &bm & fm 


( 25 ) 


5 DETERMINATION OF THE FIBER LENGTH EFFI- 
CIENCY FACTOR 


It has been claimed [9] that in the post-cracking stage the combined efficiency factors due 
to both length and orientation can not be simply calculated as the product of the length 
efficiency factor and the orientation efficiency factor because the orientation efficiency factor 
is also a function of the fiber length in the case of short fibers. For the elastic model of 
pre-cracking stage as in the present case however, these two factors can be considered to be 
independent of each other and therefore can be determined separately. 

The length efficiency factor, specifying the effect of a definite length of fibers in ashort 
fiber composite, has two expressions, depending on the stress interaction mechanisms [9]. 
For the inelastic case, the most common version of this fiber length efficiency factor is 
expressed in terms of critical fiber length [1,8, 9]. For the clastic case, it can be easily 
determined based on the tensile stress distribution. From Equation 1, the average tensile 
stress over the length of this short fiber can be calculated as 


Gjl = Eftf[ 1 - 


tanh/?i 

pI " 


While for continuous fibers, there is 


(26) 


a ~f> = a fi = Ef e f 


(27) 


Therefore the fiber length efficiency factor can be defined as 


VI — 1 - 


t.anh/?2 


(28) 


When / — * oo, ?// = 1. 

This expression shows that, compared to the continuous fibers, the tensile stress on a 
short fiber is discounted by a factor i ]i due to limited fiber length. 

It is easy to prove that when l — * 0, ?// = 0. 


6 DETERMINATION OF FIBER ORIENTATION EFFI- 
CIENCY FACTOR 


In most published studies where the effect of fiber misalignment was considered, fibers were 
assumed either all aligned in the same direction at a fixed angle with respect to the axis 
of the composite loading direc:tion[l,8,ll], i.e., there is no variation or spread existing in 
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fiber orientations, or distributed totally in random [3, 9]. Although it is usually desirable 
to orient the fibers to enhance stiffness and strength properties, in short-fiber composites, 
it is normally very difficult, if not impossible, to achieve perfect alignment or completely 
random distribution of short fibers. The orientation distributions of fibers in a composite 
are determined by the processing conditions. Partial fiber alignment is typical in injection 
and transfer moulded composites while planar partial random orientation is typical in sheet 
moulding compounds [2]. Therefore for most cases, variation of fiber orientation distribution 
has to be included in the study. There have been several reports [2,5,13,14] dealing with 
fiber orientation spreading. A more explicit form of the expression of this fiber orientation 
efficiency factor is obtained in the present study. 


G.l Form of Fiber Orientation Density Function 


Since it is impractical to deal with fibers of different orientations individually, a statistical 
approach is usually a better, or the only, alternative. To do this, a known form of the 
function to describe the fiber orientation probability density is the premise. 

Two cases of the fiber orientation distribution are of practical importance [2]. In the 
case of injection moulded objects, fiber orientation distribution is independent of the base 
angle (f> if the direction of flow is along the ^ 3 ( 2 ) axis. In sheet moulding compounds it is 
reasonable to assume that the short fibers all lie within a plane and the problem is reduced 
to a two-dimensional one. In either case, by properly arranging the coordinate system, the 
fiber orientation density function can be expressed as 


Q{0) = 


b 0<9<a 
0 a < 0 < it/ 2 


where 9 is the polar angle of a fiber with respect to the a: 3 axis (the loading direction), and 
a is the limit of 9 . 


6.2 Relationship Between Strains of Composite, Matrix and Fiber 

Assume the composite as a whole is subject to a strain e c which will cause strain ey in 
the fiber and e m in the matrix. It has been widely accepted that the elastic stress transfer 
mechanism is dominant at the pre-cracking stage and therefore the longitudinal displace- 
ments of the fiber and matrix interface are considered geometrically compatible. In other 
words, the matrix strain will be the same as the composite strain before cracking. The fiber 
strain however is dependent on the fiber orientation with respect to the loading direction. 

There are several approaches in finding the relationship between the composite strain 
and fiber strain, such as the tensor transformation method[ 10 ] and the affine deformation 
model [5]. However a few simple differentiation operations as shown below can also derive 
the same result. 

Let us consider a cylinder of matrix material with height H and radius R c . Inside 
the matrix there is a fiber with length l and orientation 6 (see Figure 4-(b)). We have a 
relationship between the three variables 

l 2 = R 2 + H 2 (29) 
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Differentiating both sides gives 


It can be further expressed as 


Bringing 


2 Idl = 2 R c dR c + 2 II dll 


dl _ R 2 dRc H 2 dH 
l ~ i 2 R c + l' 2 H 


dH 

H 


= 


and 


into it yields 


dRc 

Rc 


— —V m € e 


dl 


ej = — = e c (cos^ 0 — u m siir 0) 


(30) 

(31) 

(32) 

(33) 

(34) 


where u m represents the matrix Poisson’s ratio. Note that, similar to previous analyses 
[5, 8], the effect of the fiber Poisson’s ratio has been excluded in equation 33. Equation 
34 lias been found to be consistent with both the experimental data and the results based 
on other more sophisticated analytical analysis [5]. However, since the change of the fiber 
orientation during composite deformation is neglected in the above analysis, it is preferable 
to apply equation 34 to the small strain case. 

6*3 The Result of Fiber Orientation Efficiency Factor 

Once we have the relationship between fiber strain and the overall composite strain, the 
average strain on an arbitrary fiber due to its misalignment can then be calculated as 


q = / € c (cos 2 0 — v m sin 2 9)Q(0)dd 

Jo 


Bringing the distribution function into the above equation gives 

€/ = ~[2a(l - i/ m ) + (1 + v m ) sin 2a] 

4 a 

The overall average tensile stress on this fiber thus becomes 

c 7} = r hEftf = - u m ) + (1 + ^ m )siii2a] 


(35) 


(36) 


(37) 


Furthermore because of fiber misalignment, the contribution of this fiber toward the 
composite strength will be discounted according to the equation in [9] 


M = f T][aj ) 


(38) 
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where [T] is the transformation matrix 


m = 


cos 2 a sin 2 a —2 sin a cos a 

sin 2 a cos 2 a 2 sin a cos a 

2 • 2 

sin a cos a — smacosa cos or — sm a 


( 39 ) 


and [aj] and [<ry c ] are the actual fiber stress tensor and the fiber stress tensor in the orthog- 
onal directions with respect to loading direction. For the present uniaxially loading case, 
the above equation reduces to 

Ofc = of cos 2 9 — 2 fj sin 9 cos 0 (40) 


It can be easily proven from Equation 5 that 


Tf = 0 


So we have 


j COS 


2 o 


(41) 

(42) 


The average value of this stress with respect to fiber orientation is 

rot J 

oje = j Of cos 2 9Q(9)d9 = 7 jiEfCc — 2 [2a(l — i/ m ) + (1 + v m ) sin 2a](2a + sin 2a) (43) 
The fiber orientation efficiency factor is thus derived as 


r l° = [2a(l ~ + (1 + ^m) sin 2a] (2a 4- sin 2a) (44) 

It can be proved that when a — ► 0, ijo = 1. The minimum value of i)o — • £ n> ^ is achieved 

when a — ► 7r/2 


7 THE PROPERTY COMPATIBILITY OF FIBER AND 
MATRIX IN COMPOSITES 

Now that we have determined the maximum allowable fiber volume fraction and the min- 
imum necessary fiber volume fraction, we can use these two values to examine the fiber- 
matrix property compatibility. 

Obviously, for composite design with any possible combinations of fiber and matrix, the 
criterion 

Vmax ^ Klin (4^) 

has to be satisfied. As these two values arc determined by the properties of the fiber 
and matrix as well as the spacing between fibers, Equation 45 actually provides the inter- 
relationships between all these parameters in a composite, and can hence be used to study 
the fiber-matrix property compatibility and to select proper materials for a composite. 

The easier way of using this equation is to study the boundaries enclosed by the maxi- 
mum and the minimum fiber volume fraction curves for a given property. Some examples 
will be shown in the next section. 
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8 CALCULATION AND DISCUSSION 


First of all, since we have had all the equations describing the relationships between the 
composite structural parameters and the fiber and matrix properties, a parametric study 
becomes possible to show the effects of these properties on a composite structure. The 
data used for calculation are listed in Table 1. For generality, the ratios of fiber matrix 
properties are used wherever possible. When the efFcct of a specific parameter in Table 1 
is investigated over the given range, other parameters will take the typical values provided. 
The results are illustrated in Figures 5 - 12. 


Table I. The Fiber Matrix Properties Used for Calculation 


Item 

Range 

Typical 

Strength Ratio ^ 

Modulus Ratio 

f 1 

Fiber Aspect Ratio - 
Fiber Breaking Strength 
Matrix Stress Difference ~~ &fm 

Fiber Orientation Range a 
Matrix Poisson’s Ratio v m 
Fiber End spacing Length ratio -f- 

1.5 - 5.0 [7] 

0.02 - 0.3 [7] 

5 - 200 [assumed] 

4 - 20 Gpa [9] 

0.2 - 4.0 Gpa [assumed] 

°-S 

0.2 - 0.5 [assumed] 
0.0008 - 0.05 [assumed] 

2.4 

0.03 

GO 

8 Gpa 
2 Gpa 

7T 

G 

0.3 

0.003 


Figure 5 shows the effect of the Strength Ratio ^ on the values of V max . As the 
strength ratio increases, meaning stronger fibers are used, or a weakerbonding shear strength 
between the matrix and fibers, V max is decreasing, a greater spacing between 
fibers is required in order to maintain a stable structure. Note that fiber length does not 
have significant effect on the result. 

The effect of the modulus ratio on the V max value is illustrated in Figure 6 . It is also a 
monotonically decreasing relationship. This means that a matrix with higher shear modulus 
or a less tough fiber will result in a smaller V max value, or allow greater spacing between 
fibers. In other words, fewer fibers will be needed in the structure. Again there is no 
noticeable difference for different fiber length cases. 

Figure 7 and 8 show the relationships between the fiber length efficiency factor 77 / and 
the fiber aspect ratio and between the fiber orientation range a and the fiber orientation 
efficiency factor rjo respectively. As shown in the results, increasing the fiber aspect ratio ( a 
longer or thinner fiber ) will raise the fiber length efficiency factor, and a wider spread ( a 
greater a value ) of fiber orientation will lower the fiber orientation efficiency factor. 

Figures 9 and 10 on the other hand indicate the effects of the matrix Poisson’s ratio v m and 
the fiber orientation range a on the value of minimum fiber volume fraction V m i n . When a 
becomes larger, the value of V m { n will increase as shown in Figure 10, indicating that the fiber 
reinforcing function is hampered due to fiber misalignment so that more fibers are needed. 
A similar trend is found between v m and V 77n *„ in Figure 9 except that the relationship 
appears to be linear. 

The curves in Figures 11 and 12 can be used to test the property compatibility between 
the fibers and matrix. First of all, Figure 11 shows the effects of the modulus ratio on the 
values of V mnx ( the same curve as the short fiber case in Figure 6 ) and U min . Unlike V maxy 
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Knin decreases very slightly when ^ is increasing. Based on Equation 45, only those fiber 
and matrix types whose ^ values are greater than the critical ^ value are compatible for 
being selected to form a properly functioning composite. Likewise in Figure 12, first, the 
effects of £ on both volume fraction values V max and can be seen, showing different 

trends but both gradually approaching its own asymptote as ~ increases. On the other 
hand, as shown in the figure, there is a critical £ value above which a feasible structure can 
then be made. 


9 CONCLUSIONS 

The stress transfer between matrix and fibers in a composite is not only determined by the 
intrinsic properties of fiber and matrix, but also affected by the geometric parameters of 
fiber arrangement within the matrix such as the spacing between fibers and the orientation 
of fibers. Consequently the shear strength of the interface between fibers and the matrix can 
be used as a criterion to determine the spacing between fibers in a composite. 

For a composite made of given fiber and matrix materials, there is an optimal spac- 
ing between fibers at which the fiber tensile strength will be fully exploited. Moreover 
this optimal spacing is also the minimum allowable spacing between fibers below which 
the structure will start to disintegrate under loading before the fiber tensile failure. This 
minimum spacing then defines a maximum fiber volume fraction allowable for a composite. 

The maximum fiber volume fraction combined with the minimum fiber volume fraction 
studied previously can be used for composite design. Both volume fractions are found 
dependent on such parameters as fiber modulus Ef, fiber tensile strength <7^/, fiber aspect 
ratio ~ and fiber orientation range a, the matrix properties as Poisson’s Ratio i/ m , shear 
modulus G m , and the bonding shear strength r s between fiber and matrix, as shown in 
this study. Consequently, these two values of fiber volume fraction V max and V m { n can 
be applied to define the boundaries in determining the property compatibility of various 
combinations of fiber and matrix types for a particular application so as to optimize the 
result of composite design. 
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(d) Piber with nearest square-packed neighbours 


Figure 1 Fiber embedded in a matrix and the fiber packing forms 
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Figure 5 Maximum Volume Fraction vs. Strength Ratio 
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Figure 6 Fiber Volume Fraction vs. Gm/Ef 
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ABSTRACT 

An experimental investigation was conducted to identify the key 
geometrical parameters and quantify their influence on the mechanical 
response of through-the-thickness (TTT) reinforced composite materials. 
Composite laminates with TTT reinforcement fibers were fabricated using 
different TTT reinforcement materials and reinforcement methods and 
laminates were also fabricated of similar construction but without TTT 
reinforcement fibers. Coupon specimens were machined from these 
laminates and were destructively tested. 

Through-the-thickness (TTT) reinforcement yarns enhance damage 
tolerance and improve interlaminar strength. Thick-layer composites 
with TTT reinforcement yarns have equal or superior mechanical 
properties to thin-layer composites without TTT reinforcement yarns. A 
significant potential exists for fabrication cost reduction by using 
thick-layer composites with TTT reinforcement yarns. Removal of the 
surface loop of the TTT reinforcement improves compression strength. 
Stitching provides somewhat higher mechanical properties than integral 
weaving . 
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OBJECTIVE AND APPROACH 


The objectives of the present invest igation, as presented in Fig. 1, are 
to identify the key geometrical parameters that inf luence ^ the mechanical 
response of through-the-thickness (TTT) reinforced composite materials 
and, where possible, to quantify the influence of these key parameters 
on the mechanical response. The approach used to identify and quantify 
the effects of key geometrical parameters was based upon observations 
and destructive testing of coupon specimens, respect ively . . Composite 
laminates with TTT reinforcement fibers were fabricated using different 
TTT reinforcement materials and reinforcement methods. These laminates 
were microscopically examined to identify potential geometrical features 
that influence mechanical response. Coupon specimens were machined from 
these laminates and were destructively tested. . Test results were 
compared with test results from materials of similar construction but 
without TTT reinforcement fibers. 


Objective: 

• Identify key geometrical parameters and quantify 
their influence on the mechanical response of 
through-the-thickness (TTT) reinforced composite 
materials. 

Approach: 

• Investigate specimens having different TTT 
reinforcement materials and reinforcement methods. 

• Identify potential material and geometrical features 
that influences mechanical response. 

.Investigate the influence of TTT reinforcement 
of materials having thick ply layers. 


Figure 1 


FEATURES OF TTT REINFORCEMENT THAT INFLUENCE MECHANICAL RESPONSE 


There are four geometrical features of TTT reinforced composite 
materials that significantly influences their mechanical response, as 
shown in Fig. 2. These geometrical features are resin rich regions, in- 
plane fiber waviness, surface loop of continuous TTT reinforcement 
fibers and breakage of in-plane fibers. 


• Resin rich regions 

• In-plane fiber waviness 

• Surface loop of continuous through-the- 

thickness reinforcement 

• Breakage of in-plane fibers 

Figure 2 
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PHOTOMICROGRAPHS OF STITCHED AND INTEGRALLY 
WOVEN COMPOSITE MATERIAL 

The resin rich regions identified in Fig. 3 were created by the 
inclusion of TTT reinforcement fibers. The stitch yarn creates a small 
inclusion around each TTT reinforcement penetration. The TTT 
reinforcement fibers of the integrally woven material create resin 
channels because the integrally woven TTT reinforcement fibers separate 
adjacent in-plane yarns . It is well understood that changes in fiber 
volume fraction can change the mechanical properties of composite 
materials (ref .1) . Therefore, it is reasonable to believe that resin 
rich regions created as an artifact of inserting TTT reinforcement 
fibers could influence both laminate stiffness and strength. 

Fiber waviness in composite laminates has also been identified as a 
mechanism that can adversely influence compression strength (ref. 2). 

The vertical portion of the TTT reinforcement fibers creates in-plane 

waviness of the in— plane reinforcement fibers, as observed m Fig. 
3, for both stitched and integrally woven TTT reinforced laminates. The 
amplitude of the waviness is approximately equal to the diameter of the 
TTT reinforcement fibers which is generally greater than the fiber 
waviness found in non-textile composite laminates. 

When TTT reinforcement fibers penetrate in-plane yarns, such as in the 
stitched laminate depicted in Fig. 3, the penetration by the needle and 
stitching yarn can result in breakage of in-plane fibers. With proper 
processing the majority of the in-plane fibers at stitch penetration 
are pushed aside as the needle penetrates the preform. Generally, less 
than 5 percent of the in-plane fibers at the penetration site are broken 
(ref. 3). No significant breakage of in-plane fibers occurs for 
integrally woven TTT reinforced materials because^ the TTT reinforcement 
does not penetrate the in-plane yarns and hence little or no breakage of 
in-plane reinforcement fibers occurs in the weaving process. 


(Figure 3 appears on the next page.) 
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Figure 


KINKED IN-PLANE SURFACE YARN 


Both the stitched and integrally woven materials studied in this 

investigation use continuous TTT reinforcement fibers. Between 

successive TTT fiber penetrations, a part of the TTT reinforcement 

fibers extends across the laminate's surface. The portion of the TTT 

reinforcement fibers that lies on the surface is referred to as the 

surface loop, see Fig. 4. When the preform is Infiltrated and processed 

into a composite part, the surface loop of the TTT reinforcement fibers | 

is forced into the surface - oT" the laminate causing a kinking or bending j 

of the laminate's in-plane fibers near the laminate surface. The 5 

magnitude of the kinking or bending of the in-plane fibers is a function § 

of the diameter of the surface loop yarn and the compaction pressure j 

used during consolidation. Only the layers of in-plane fibers near the I 

surface of the laminate are influenced by the surface loop. That is, s 

the layers near the center of the laminate are not kinked or bent by the 

surface loop. I 




(Figure 4 appears on the next page.) 
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KINKED IN-PLANE SURFACE YARN 
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Photomicrograph of laminate cross section 


MATERIALS INVESTIGATED 


All laminates studied in this investigation were fabricated using 
Hercules, Inc. AS4 carbon fiber with Hercules, Inc. 3501-6 epoxy resin 
and were approximately 0.64 cm thick. See Fig. 5- The stacking sequence 
for all laminates consisted of variations of a 0/90 cross-ply stacking 
sequence. However, the per ply thickness, stacking sequences, inclusion 
of TTT reinforcements, TTT reinforcement method and TTT reinforcement 
material differed between laminates. The 0/90 stacking sequence was 
chosen so that a direct comparison between materials with stitched and 
integrally woven TTT reinforcement fibers could be made. Three 
different laminates were fabricated without TTT reinforcement fibers to 
quantify the influence of in-plane fiber architecture and TTT 
reinforcement fibers on the mechanical response. Two of these laminates 
were fabricated from unidirectional prepreg material. Their stacking 
sequences are [ (0/90) 2 /0/ (0/90) 5 /0/ (0/90) 3 /0] s and [ (0 5 /90 5 ) 2 /O 3 ] s and 
they are referred to herein as the thin-layer and thick-layer material, 
respectively. The third material without TTT reinforcement, referred to 
as the uniweave material, was fabricated from a nine layer stack of dry 
uniweave fabric that was infiltrated and cured to form the composite 
laminate. This uniweave fabric has approximately 99 percent of its 
carbon fibers oriented in the fabric's warp direction. The remaining 
one percent of the fibers are fine denier glass yarns used to hold the 
carbon warp yarns together. The glass yarns are oriented in the 
fabric's fill direction. Each layer of the uniweave fabric is 
approximately equal to five layers of unidirectional prepreg tape 
material. The uniweave fabric is composed of 21K-filament carbon yarns 
with a yarn spacing of approximately 5 per cm. The thick-layer and 
uniweave materials are similar in layer thickness and construction. 

All laminates with TTT reinforcement fibers have the same ply 
orientation as the uniweave material, that is [0/90/0/90/0/90/0/90/0], 

Two different types of TTT reinforcement yarn were used, namely carbon 
and Kevlar yarns. The carbon yarn is a Toray T-900-1000A stitching yarn 
and the Kevlar yarn is a 1100 denier Kevlar-49 yarn. Both the stitched 
and integrally woven materials were made with the same number and 
spacing of TTT reinforcement yarns. The stitched preforms were produced 
using uniweave fabric. All laminates fabricated with TTT reinforcement 
fibers were fabricated by stitching or integral weaving. Both stitching 
and integral weaving are referred to herein as continuous TTT 
reinforcement processes. Simulated tufted TTT reinforcement was 
produced by machining the surfaces of stitched or integrally woven 
laminates to remove the surface loop of the TTT reinforcement fibers. 

The simulated tufting TTT reinforcement is referred to herein as a 
discontinuous TTT reinforcement process. 

AS4-3501-6, 0.64 cm thick 


Laminates without through-the-thickness (TTT) 
reinforcement 

Thin layer - [(0/90y0/(0/90) 5 /0/(0/90y0] s 
Thick layer - [(ty9^y<y s 
Uniweave - [0/90/0/90/0/90/0/90/0] 

Laminates with TTT reinforcement 
[0/90/0/90/0/90/0/90/0] 

TTT fiber 

Carbon 
Kevlar 


Figure 5 


TTT method 

Continuous 

Stitched 

Woven 

Discontinuous 
Simulated tufting 



PREFORM ARCHITECTURE 


Three different TTT reinforced preform architectures were used in this 
study, as seen in Fig. 6. The stitched preforms were created by 
stacking nine layers of uniweave fabric in a 0/90 orientation and the 
stack was stitched using a modified through-the-thickness lock stitch in 
orthogonal rows and columns oriented parallel and perpendicular to the 0 
degree direction of the preform, respectively. The integrally woven 
preforms were produced as a single unit on a loom. The in-plane warp 
and fill yarns were not interwoven. The nine layers of fabric were held 
together by the TTT reinforcement yarns woven in both the warp and fill 
directions. The TTT reinforcement yarns were inserted and wrapped 
around a catcher yarn in the center of the preform and returned to the 
surface. Therefore, the TTT reinforcement yarn in the integrally woven 
material was actually composed of two yarns, a yarn originating from the 
upper surface and one originating from the lower surface. The simulated 
tufting material was produced by machining away one half of the upper 
and lower outer surface layers from the stitched and integrally woven 
cured composite laminates. 



No-crimp stitched Integrally woven 3-D weave Simulated tufting 


Figure 6 


131 



liiHllllllliillll I il 


TEST SPECIMENS 


Three different test specimens were used to evaluate the mechanical 
response of the different materials. They were the compression-after- 
impact (CAI) specimen, the short-block-compression (SBC) specimen and 
the multi-span-beam (MSB) specimen, as depicted in Fig. 7. The CAI and 
SBC specimens were machined from laminates such that the direction of 
the applied load coincided with the 0 degree direction of the laminate 
stacking sequence. The MSB specimens were machined from laminates such 
that the 90 degree direction of the laminate stacking sequence was 
parallel to the length of the specimen. 

The CAI specimen was used to evaluate the damage resistance and residual 
compression strength of a damaged panel. The CAI specimen was 12.7 cm 
wide by 25.4 cm long. A 0 . 127-cm-diameter aluminum ball was shot at 
the center of each CAI specimen with approximately 42J of energy. Back- 
to-back strain gages were bonded to each CAI specimen. Each CAI 
specimen was then mounted into a side-supported compression test fixture 
and statically compressed to failure. 

The SBC specimen was used to evaluate the undamaged compression strength 
of the material. The SBC is 3.81 cm wide by 4.45 cm long. Each 
specimen was instrumented with back-to-back strain gages. The SBC 
specimens were mounted in a compression test fixture without side 
supports and a static compression load was applied until the specimen 
failed. 

The MSB specimens were used to evaluate, in a qualitative manner, the 
interlaminar strength of the material. The MSB test consists of an 
upper and lower load introduction structure which makes contact with the 
specimen at five locations across the width of the specimen. The width 
of the specimen is 2.54 cm and the length is 12.7 cm. A static 
compression load is applied to the specimen until failure. 
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Figure 7 


EFFECTS OF TTT REINFORCEMENT ON COMPRESSION STRENGTH 

Short -block-compression tests were performed on the thin-layer, thick- 
layer, uniweave materials as well as the carbon and Kevlar TTT 
reinforced stitched and integrally woven materials. The highest 
compression strength was obtained from the thin-layer material followed 
bv the uniweave, thick-layer materials and the materials with TTT 
reinforcement yarns. See Fig. 8. The difference in strength between the 
three materials without TTT reinforcement yarns is related to the 
architecture of the materials. Increasing layer thickness produces a 
corresponding increase in the magnitude of interlaminar stresses, 
therefore the thin-layer material should have higher compression 
strength than the thick-layer material. The thick-layer and uniweave 
materials have similar layer thickness, however the uniweave had higher 
strength. This difference in strength is attributed to the glass fill 
yarn which partially intertwines each warp yarn in the uniweave material 
and produces some intralaminar reinforcement. 

All of the specimens with TTT reinforcement yarns have significantly 
lower undamaged compression strength than those materials without TTT 
reinforcement yarns. The stitched material exhibited somewhat _ higher 
undamaged compression strength than the integrally woven materials. The 
strength difference is attributed to the reduced amount of kinking or 
bending of the in-plane fibers adjacent to the surface loop formed by 
the needle yarn in the stitched material. The needle yarn used in the 
stitched material was a small denier yarn which produced less kinking or 
bendinq of the in-plane yarns when the preform was processed into a 
composite part. The type of reinforcement yarn in the TTT reinforced 
laminates had no influence on compression strength. 
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Thin layer 
Thick layer 
Uniweave 


Compressive 

strength, 
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Carbon TTT reinforcement 
Kevlar TTT reinforcement 



reinlorcemenl 


Figure 8 



FAILURE MODES OF SHORT -BLOCK-COMPRESS ION SPECIMENS 


The compression failure mode for the thin-layer, thick-layer and 
uniweave materials was primarily delamination which precipitated a 
transverse shear failure, as shown in Fig. 9. This failure mode was 
caused by interlaminar stresses. All materials with TTT reinforcement 
yarns exhibited a transverse shear failure mechanism because the TTT 
reinforcement yarns increase the interlaminar strength of the laminate 
which suppresses the delamiriation failure mechanism. 


(Figure 9 appears on the next page.) 
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Figure 





INFLUENCE OF TTT SURFACE LOOP ON COMPRESSION STRENGTH 

Short -block -compress ion tests were conducted on stitched and integrally 
woven composite laminates. One-half of the specimens had one-half of 
their surface layers machined away to simulate tufting. The simulated 
tufted specimens are referred to herein as those "without TTT surface 
loop”. In all cases removal of the TTT reinforcement surface loop 
and the associated kinked or bent in-plane fibers resulted in higher 
compression strength.^ See Fig. 10 . Compression strength was based upon 
the actual cross-sectional area of the specimen, that is the reduced 
area in the case of machined specimens. The axial-load-carrying 
capability is primarily due to the 0 degree fibers. Approximately 12 
percent of the 0 degree fibers were removed in the machining process 
because the outer layers were oriented in the 0 degree direction. Since 
the machining produced a net increase in strength, the 0 degree fibers 
that were machined away were ineffective in carrying a compressive load. 

The stitched material with the TTT surface loop had higher compressive 
strength than the integrally woven material. However, after the 
surfaces were machined there was little or no difference in compression 
strength between the stitched and integrally woven materials. The 
stitched materials had the surface loop on the bobbin yarn surface 
created by the carbon or Kevlar TTT reinforcement yarns, whereas the 
surface loop produced on the needle yarn surface was made with a much 
smaller denier yarn. The surface loops produced in the integrally woven 
materials were produced using the same type carbon or Kevlar yarn on the 
outer surfaces . Therefore, more in-plane fibers in the integrally woven 
materials were kinked or bent than in the stitched materials, hence the 
lower compression strength for the integrally woven materials. 
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FAILURE MODES OF SBC AND CAI SPECIMENS 

The failure modes of SBC and CAI specimens having TTT reinforcement 
yarns with and without the surface loop of the TTT reinforcement yarn 
are shown in Fig. 11. Both the unmachined (with surface loop) and 
machined (without surface loop) specimens have similar transverse 
shearing failure modes. 



Figure 11 
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EFFECTS OF TTT REINFORCEMENT 
ON COMPRESSION STRENGTH OF IMPACTED PANELS 

The materials with TTT reinforcement yarns exhibited nearly twice the 
CAI strength as materials without TTT reinforcement yarns, as shown in 
Fig. 12. The typical failure mechanism of the CAI panels is controlled 
by the interlaminar strength of the material. Inclusion of TTT 
reinforcement yarns improves the interlaminar normal and shear strengths 
of the material and produces a higher CAI strength. Stitched materials 
have slightly higher CAI strengths than the integrally woven materials. 

This trend is consistent with the undamaged compression strength of 
these materials. The thin-layer material exhibited higher CAI strength 
than either the uniweave or the thick-layer materials. The relative 
magnitude of the CAI strength of materials without TTT reinforcement 
yarns was consistent with the undamaged compression strength for these 
materials . 
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EFFECTS OF TTT SURFACE LOOP 
ON COMPRESSION STRENGTH OF IMPACTED PANELS 

Neither the CAI strength nor the impact-induced damage area were 
affected by the removal of the TTT surface loop, as shown in Fig. 13. All 
panels were impacted at approximately 42J of energy,’ however the panels 
with the surface loop removed were 11 percent thinner. It was Initially 
expected that the thinner panel would have a greater damage area due to 
impact but the results indicate that this was not the case. These 
results suggest that the surface loop has no positive influence on the 
CAI strength. Futhermore, the surface loop has no influence on the 
damage containment because the damage area of the specimens with and 
without a surface loop was approximately the same. 
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EFFECTS OF TTT REINFORCEMENT ON MULTI-SPAN-BEAM STRENGTH 

The MSB strength of the thin-layer material is considerably higher than 
the MSB strength of either the thick-layer or the uniweave materials as 
shown in Fig. 14. This ordering of strength is consistent with the 
compression and CAI strength for these materials without TTT 
reinforcement yarns because the failure mechanism is induced by 
interlaminar stresses. However, the MSB strength of the thick-layer 
material with carbon TTT reinforcement yarns was 60 percent higher than 
the thick-layer and uniweave materials and was 90 percent of the MSB 
strength of the thin-layer material. These results suggest that when 
properly designed, thick— layer material with TTT reinforcement yarns can 
be used in lieu of thin-layer materials. A significant cost savings ^ 
(material and processing) can be realized by using thick-layer material 
with TTT reinforcement yarns. 
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INFLUENCE OF TTT REINFORCEMENT SURFACE LOOP 
ON MULTI-SPAN-BEAM STRENGTH 

The surface loop had no effect on MSB strength for either type of TTT 
reinforcement or reinforcement method as seen in Fig. 15. This result 
is reasonable because the region of high interlaminar shear and the 
initial failure site are in the center of the laminate. The surface loop 
influences the mechanical properties of the in-plane fibers along the 
outer surface of the specimen. Therefore, it is reasonable to expect 
that the surface loop has no influence on the MSB strength. 

The MSB strength of the stitched material is consistently higher than 
the MSB strength of the integrally woven material. The MSB failures 
always initiated in the interior of the beam in a 90 degree ply or in 
one of the resin pockets. The integrally woven materials have larger 
resin pockets than the stitched material and have TTT yarn loops around 
a catcher yarn along the specimen's centerline. It is suspected that the 
larger resin pockets and the presence of the catcher yarns are 
responsible for the lower MSB strength of the integrally woven 
materials . 
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CONCLUDING REMARKS 


An experimental investigation was conducted to identify the key 
geometrical parameters and to quantify their influence on the mechanical 
response of through— the— thickness (TTT) reinforced composite materials. 
Composite laminates with TTT reinforcement fibers were fabricated using 
different TTT reinforcement materials and reinforcement methods. These 
laminates were microscopically examined to identify potential 
geometrical features that influence mechanical response. Coupon 
specimens were machined from these laminates and were tested to failure. 
Test results were compared with test results from materials of similar 
construction but without TTT reinforcement fibers. The concluding 
remarks are summarized in Fig. 16 and discussed further below. 

Through-the-thickness reinforcement yarns enhance the damage tolerance 
and improve interlaminar strength of composite materials. However, TTT 
reinforcement yarns in composite materials cause a reduction in 
undamaged compression strength. The reduction of undamaged compression 
strength due to incorporating TTT reinforcement yarns in composite 
materials is a result of local resin rich regions, in-plane fiber 
waviness, surface loop of the TTT reinforcement ^ yarn kinking or bending 
the in-plane fibers, and breakage of in-plane fibers. 

The thick-layer composites studied in this investigation with TTT 
reinforcement yarns have mechanical properties that are equal to or 
superior to thin-layer composites without TTT reinforcement. Thick- 
layer materials can be fabricated from large filament count yarns at a 
lower cost than with small filament count yarns. Cost of large filament 
count reinforcement yarns is substantially less than the cost of small 
filament count reinforcement yarns. Fewer thick layers of^ material are 
required to produce a part than are parts produced from thin layer 
materials. The reduced number of layers of material reduces part layup 
cost. Therefore, there is a significant potential cost savings when 
using thick-layer composites and large filament count yarns. 


The stitched materials evaluated in this investigation had different 
denier yarn for the needle and bobbin yarns. The needle yarn was a 
smaller denier TTT reinforcement yarn than the bobbin yarn. The 
smaller denier needle yarn kinked or bent the inplane fiber less than 
the larger denier bobbin yarn. Since the inplane surface yarns adjacent 
to the needle yarn are kinked and bent less, then their compression 
strengths are degraded less. The TTT reinforcement yarns in the 
integrally woven materials were of the same denier as the larger denier 
TTT reinforcement yarn in the stitched material . The outer surfaces of 
the integrally woven laminate were equally influenced by the surface 
loop of the TTT reinforcement yarns. Therefore, stitched laminates 
provide somewhat higher mechanical properties than integral weaving. 


• Through-the-thickness reinforcement enhances damage 

tolerance and improves interlaminar strength. 

• Thick layer composites with TTT reinforcement have 

equal or superior mechanical properties to thin layer 
composites without TTT reinforcement. 

(Potential exists for fabrication cost reductions.) 

® Removal of surface loop improves compression strength. 

o Stitching provides somewhat higher mechanical properties 
than integral weaving. 


Figure 16 
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Abstract 


The importance of textile composite materials in aerospace structural ap 
plications has been gaining momentum in recent years. With a view to better 
understand the suitability of these materials in aerospace applications, an ex- 
perimental program was undertaken to assess the mechanical properties of 
these materials. Specifically, the braided textile preforms were infiltrated with 
suitable polymeric matrices leading to the fabrication of composite test 
coupons. Evaluation of the tensile properties and the analyses of the results 
in the form of strength, moduli, Poisson’s ratio, etc., for the braided compos- 
ites are presented. Based on our past experience with the textile coupons, the 
fabrication techniques have been modified (by incorporating glass microbal- 
loons in the matrix and / or by stabilizing the braid angle along the length of 
the specimen with axial fibers) to achieve enhanced mechanical properties of 
the textile composites. This paper outlines the preliminary experimental re- 
sults obtained from testing these composites. 


Introduction 


The advantages in using the engineered laminated composites haying 
high specific strength and specific modulus in the design of payload sensitive 
spacecraft have been known for some time. As the range of applications using 
these laminated composites is increasing, other design limiting problems such 
as the interlaminar strength, material toughness, fabrication and tooling costs 
need to be addressed critically. By virtue of the inherent geometrical layups of 
laminated composite structural components, non-uniform stresses in the vari- 
ous layers and directions take place under loads. These stresses in turn induce 
interlaminar shear and normal stresses at the ’free’ boundaries or edges [1] of 
composite laminates resulting in damage and/or premature failures. While de- 
sign modifications or techniques can be incorporated to suppress the delami- 
nations, nevertheless newer approaches in developing alternate and cost-effec- 
tive technologies are needed. As a result, the existing technologies from the 
textile industry are being explored and developed [2,3] in the design of 
aerospace composite components. 
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Many of the textile manufacturing processes such as braiding, knitting, 
and weaving will reinforce through the “thickness” of the composite structural 
components, thereby virtually eliminating or minimizing the effect of delamina- 
tions as well as giving rise to additional strength, if necessary, between the lay- 
ers. Further, it's believed that the textile fabrication technologies could con- 
tribute to lower costs. As a result, the research and developmental activi- 
ties in the use of textile technologies for structural composites are gaining mo- 
mentum in developing multiaxial fabrics, integrally woven structures with stiff- 
eners, and near -net-shape preforms [4,5]. These textiles preforms assume the 
near net structural shapes after impregnating and curing with suitable matrix 
with very little, if any, further processing operations. In addition, the textile 
preforms in their near -net-shape may be dry, prepregged or impregnated with 
(matrix) powder or commingled with matrix filaments. The latter two tech- 
niques by-pass the (liquid) resin transfer molding operation. One of the attrac- 
tive features in using the net-shape preform technologies is its amenability for 
automation such as pultrusion, where production quantities justify, thereby 
maintaining quality control and lowering the cost of fabrication. 


Materials, Processing, and Fabrication 


The graphite textile preforms were braided with tows of Celion G30-500 
(BASF Structural Materials, Inc.) fibers. These fibers measure 7 microns in di- 
ameter and have a modulus of 30 Msi. The preforms were braided by a 4-step 
process where each tow was interlocked with the other tows to form a true 3-D 
structure. The braid angle was varied between 17 and 30 degrees. 

The preforms were fabricated using a vacuum/compression molding con- 
solidation process. The matrix system was comprised of two components- 
Epon 828 (Shell Chemical Co.) and Jeffamine T-403 (Texaco Chemical Co); 100 
parts-by -weight to 42 parts-by -weight, respectively. This resin system was cho- 
sen because it had a relatively long pot life and low viscosity which was neces- 
sary for this method of molding. The two components were heated separately 
at 60° C for ten minutes to reduce the viscosity which in turn reduces the time 
required to evacuate the catalyzed resin. The two components were then com- 
bined and evacuated for ten minutes after the vacuum pressure reached one 
Torr. The bottom of the stainless steel mold cavity was completely covered to a 
depth of about 1/16” with the evacuated resin. The textile preform was then 
placed into the mold cavity and was covered with the remaining resin, and 
again evacuated for ten minutes after the vacuum pressure reached one Torr. 
Once the evacuation process was complete, the plunger was placed in the mold 
cavity and the mold was placed in a heated press where the excess resin was 
squeezed out in order to obtain a desired composite thickness. The nominal 
specimen dimensions are 0.100 x 0.750 x 10 inches. The composites were 
cured for three hours at 100° C. The fiber volume of the specimens was found 
to vary between 45% and 55%. The tensile coupons were tabbed using 1 /8” 
fiberglass tabs. The tabs were attached with an adhesive film, FM 123-5 
(American Cyanamid), which required a cure cycle of 90 minutes at 100° C and 
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50 psi. The geometry and dimensions of four different types of braided test 
specimens are shown in Figure 1 . 


In an effort to lower the density of the composite specimens, glass mi- 
croballoons were introduced into the matrix of some of the test coupons. 
Glass microballoons can be thought of as controlled voids which will displace 
some of the resin. Thus, the idea was to use controlled voids to reduce the 
density of the samples without any significant loss in many of the properties. 
The glass microballoons used in the study had a diameter range from 7 to 
microns. The resulting density reduction was found to be about 2 /o. The spec 
™en without microbailoons had a density of 1.51 g/cm> and with nucrobal- 
loons it was 1 48 g/cm 3 . The test results indicate that the resin displaced by 
the glass microballoons (5% by weight) did not significantly affect the tensile 
strength of the composite specimen. However, a difference m the stom behav- 
ior was observed. Further studies are suggested to optimize the strength and 
density with respect to the percentage of microbailoons that can be incorporat- 
ed in the composite materials. 


Experimental Results and Analyses 


Numerous factors affect the performance of braided textile composites. It 
is not uncommon to observe significant variations in the ^mechanical Properties 
of braided composites. At the laboratory coupon-level fabrication, the quality 
control of textile composites becomes even more critical if one is assessing 
relative merits of these composites for their mechanical properties. Several fac- 
tors such as the quality of braided preform, stabilizing the braid angle throug 
the molding process, extent of curing, compaction of the preform maintaining 
precise fiber volume ratio, size effect, uniform density, location of the 
cut from a finite length of the preform, complete removal of voids and tiny 
pockets of trapped air bubbles, uniform distribution of fiber tows and the sur- 
rounding matrix throughout the length of the test specimen, ^nmetry of the 
composite, dimensional variations in the as-molded specimen, etc., n ay 
adversely affect the mechanical properties. 

Preliminary experimental data of the tensile properties of four -different 
types of braided composites are shown in Table I. In j^ ese nting thii s lnitia 
data no attempt has yet been made to censor the data based on the factors 
mentioned in the preceding paragraph that may adversely affect die ^™^amca 
properties Further analyses of the raw data are underway to refine and com 
pare the properties of the four types of specimens with the factors that influ- 
ence the results in the background. 

The objective in introducing the glass microbailoons in the composite 
was to not only lighten the braid but to reduce the formation of resin pockets 
Tat may develop at the crossings of large-sized (12k) tows. A cursory exami- 
nation of the resulting data as shown appears to indicate the nominal 

tensile properties were not affected except for the transverse strain Based on 
density and fiber volume ratio measurements to be performed shortly, the eval 
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GLASS / EPOXY 
COMPOSITE TAB 


TYPICAL BRAIDED SPECIMEN GEOMETRY 

TYPICAL DIMENSIONS OF FOUR TYPES OF SPECIMENS* 



a 

b 

c 

d 

e 

f 

g 

Plain braids 

0.750 

0.352 

0.102 

2.75 

4.5 

0.25 

2.5 

Braids with 
microballoons 

0.753 

0.354 

0.104 

2.75 

4.5 

0.25 

2.5 

Braids with 
axial fibers 

0.877 

0.354 

0.104 

2.75 

4.5 

0.25 

2.5 

Braids with 
axial fibers & 
microballoons 

♦ All : 

0.878 

0.362 

T \ • 1 

0.1 12 
ill 

2.75 

4.5 

0.25 

2.5 


* All dimensions in inches. Braids with axials have 37% (26 tows out of 70) 
axial fibers. 


Figure 1. Tensile Specimen Geometiy and Dimensions 
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TABLE-1 : TENSILE PROPERTIES OF BRAIDS* 



Braid 

Load 

Stress 

Strain 

Modulus 

Trans. Strain 

Poisson’s 


angle 

Pm ax 

(7 max 

£ 1 1 max 

El i 

£12 

ratio 


(deg) 

(kip) 

(ksi) 

(%) 

(Msi) 

(%) 

V 12 

Plain braids 

23/30 

13.22 ±1.7 

173.5+23 

1.47+0.09 

11.8 

2.31±0.36 

1.57 

Braids with 
microballoons 

21/30 

13. 44 ±2.3 

172.1+29 

1.52±0. 1 1 

11.3 

2.69±0.40 

1.77 

Braids with 
axial fibers 

- 

15.81+1.0 

179.0+14 

1.34±0.07 

13.4 

1 .19+0.27 

0.89 

Braids with 
axial fibers & 


13.14+2.3 

135.3±27 

1.26+0.10 

10.7 

1.39+0.33 

1.10 


microballoons 


» The experimental data shown here in tabular form is for convenience only. Comparison of 
data among the various cases shown here is not recommended. Extensive data reduction and 
detailed analyses are planned for the purpose of comparing the properties on a rational basis. 
Factors such as braid angle, fiber volume ratio, failure zones, size effect, density, specimen 
imperfections, if any, etc., are expected to influence the preliminary data shown above. 
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uation of specific mechanical properties will be conducted. If the succeeding 
analyses prove to be encouraging, a corollary that needs to be addressed is the 
percentage of microballoons that could be introduced in a unit volume of the 
braid resulting in one or more than one designated optimum property or proper- 
ties. Further, what is the effect of introducing other grades of microballoons 
as to size and specimen thickness on these properties ? 

A nonuniform drift in the braid angle was noticed which may be at- 
tributed to the compression molding process. In order to stabilize the braid 
angle as well as to further improve the mechanical properties of the braids si- 
multaneously, axial fibers(tows) were introduced during the preform fabrica- 
tion. As expected, preliminary data reveal an improvement in the nominal 
tensile strength and modulus, and a decrease in the strains. The percentage of 
axial fibers selected in this study was random. Optimization techniques may 
be applied to determine the percent of axial fibers that can be embedded in 
braided specimens. 

The experimental results with a matrix containing microballoons and 
braids containing axials were obtained very recently and are being analyzed 
critically to understand their significance. 


Conclusions 

1. Four different types of braided specimens were studied to assess their 
tensile mechanical properties. Factors that affect the properties are 
identified for further data reduction and analyses. 

2. A preliminary evaluation on the effect of microballoons in the plain 
braided composites on the mechanical properties is presented. Many 
of the properties appear to remain unchanged. Further analyses based 
on density measurements may show an improvement in some of the 
specific properties. 

3. Axial fibers were introduced in the braids to improve the mechanical 
properties further, and the preliminary results are encouraging. 

4. Through these studies, it has become clear that many factors affect the 
mechanical properties ranging from preform and composite fabrication 
through testing and rational data reduction. 
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MODELING OF EUCLIDEAN BRAIDED FIBER ARCHITECTURES TO OPTIMIZE COMPOSITE PROPERTIES 
E. Armstrong-Carroll *, C. Pastore^, and F. K. Ko^ 


ABSTRACT 

Three-dimensional braided fiber reinforcements are a 
very effective toughening mechanism for composite 
materials. The integral yarn path inherent to this fiber 
architecture allows for effective multidirectional dispersion 
of strain energy and negates delamination problems. In 
this paper a geometric model of Euclidean braid fiber 
architectures is presented. This information is used to 
determine the degree of geometric isotropy in the braids. 
This information, when combined with candidate material 
properties, can be used to quickly generate an estimate of 
the available load-carrying capacity of Euclidean braids at 
any arbitrary angle. 

INTRODUCTION 

Three-dimensional braided fiber architectures generate 
thick fabrics with yarns traveling diagonally through the 
fabric thickness. The outstanding performance features of 
these braids are high damage tolerance and delamination 
resistance. The mechanical properties of these materials 
are documented in previous works |1, 2|. 

Correlations between braid fiber architecture and 
resultant mechanical properties have been made [2|. 
Further efforts have been made to develop a processing 
science model of Euclidean braid fiber architectures [3|. 
The processing model incorporates yarn geometry and other 
fabrication variables in a geometric cross-sectional slice 
model. The processing model correlates fabric design 
inputs parameters with a geometric model in order to more 
accurately determine the resultant yarn orientations and 
fabric dimensions. 

This paper presents a unit cell model of Euclidean braids 
in the close-packed condition. The dependence of fiber 
volume fraction on braid angle is formulated. A geometric 
isotropy model [4] is applied to the unit cell model. This 
application generates plots of the distribution of effective 
fibrous reinforcement in Euclidean braids. These plots 
provide a quick and effective guideline for selecting the 
proper braid angle for a specific application. 

FABRIC FORMATION PROCESS 

Braids are formed by the intertwining of yarns. This 
intertwining is accomplished by the crossing of yarns on 
individual yarn carriers. Three-dimensional Euclidean 
braids are formed in a four-step braiding sequence. In step 
one columns of yarn carriers are moved up and down 
relative to each other. In step two tracks of yarn carriers 
are moved back and forth relative to each other. In step 
three columns of yarn carriers are moved up and down 
relative to each other. In step four tracks of yarn carriers 
are moved back and forth relative to each other. The 
correspondence between track and column loom movements 
and yarn movements in forming a fabric are shown in 
Figure 1 . 



Step 

Four 


Step 

Five 



Figure 1. Track/column loom motion used to form an 
Euclidean braid and the initial unit cell 
produced by this process 


After each set of track and column movements, the 
yarns are compacted. In this process body diagonal yarn 
pairs resulting from a track/column movement are 
compacted against body diagonal pairs arising from the 
previous track/column motion. The compacting motion 
intertwines the yarns. 

The yarns transverse through the thickness of the fabric 
in a zig-zag motion as the braiding sequence is repeated. 
Figures 2 and 3 project the resultant yarn paths into the loom 
plane for rectangular and circular looms. The three 
dimensional path of one yarn in a Euclidean braid is depicted 
in Figure 4. In this figure a discrete lattice is used to locate 
the yarn in the braid. The presence of this lattice has 
generated the nomenclature, Euclidean braiding, to describe 
the track/column braiding process. 


SBBBBBBBB 


U 


FOOTNOTES: 

* Naval Air Development Center, Warminster, PA 
t North Carolina State University, Raleigh, NC 
tt Drexel University, Philadelphia, PA 


Figure 2. Yarn path projected onto the braiding plane 
of a rectangular loom |3| 


preceding page blank not filmed 
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Figure 3. Yarn path projected onto the braiding plane 
in a circular loom [3| 



Figure 4. Isolated path of a single yarn in a 3D 
braided fabric [3 J 


UNIT CELL MODEL 

This model describes the unit cell of a three-dimensional 
braid when the yarns are close-packed in the xy fabric 
formation plane. Several simplifying yarn geometry 
assumptions are made for this model. The yarns are 
assumed to be incompressible, identical, and circular in 
cross-section. 

The loom and compaction motions involved in the 
Euclidean braid fabric formation process cause the diagonal 
movement of the yarns in space. The path of these yarns 
through a unit volume element is shown in Figure 5. The 
volume is tetragonal since the yarns possess a circular cross- 
section and the yarn carriers on the loom are equally 
spaced in the x and y directions. 



Figure 5. A schematic of yarn paths in an Euclidean i 
braid j 

Figure 6 is a schematic of through-thickness cross- 
sections of an Euclidean braid. These cross-sections 
correspond to the xy planes of the unit cell. The yarns 
appear elliptical due to the braiding angle. The squares 
superimposed on the yarns represent the unit cell 
boundaries in the xy planes. The length of the cell in the 
x direction is a. The length of the cell in the y direction 
is b. The unfilled yarns pass through the fringes of the unit 
cell. The yarns progress through the unit cell as a result I 
of the loom motions. For example, as the columns move, * 

the shaded yarn in the lower left-hand corner of the 
bottom cell plane zigs halfway up the b length. The 
subsequent track movement zag and compaction place the 
shaded yarn in the central region of the unit cell. 


Bottom Cell Plane Central Ceil Plane Top Cell Plane 

I 

Figure 6. Close-packed through-thickness planes of an 1 
Euclidean braid I 


Thelengths of the a and b parameters are equivalent. i 
These lengths are a function of the braid angle and the yarn | 
diameter. Figure 7 defines the braid angle as it appears | 
in the unit cell. The length of the a and b parameters can | 
be defined as: I 

i 

| 
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Figure 7. The relationship of the braid angle to the 
unit cell 



a = b = (D/sin0) 4- 2D)/SQRT(2) (1) Figure 8. The effect of varying braid angle on fiber 

volume fraction 

The height of the unit cell is measured by the 

parameter c. The length of c can be defined as: DISTRIBUTION OF FIBROUS REINFORCEMENT 


c = a SQRT(2)/tan0 


( 2 ) 


Since there must be enough room in the unit cell for 
yarns to form the V-shaped and X-shaped cross-overs in the 
cell, there exists a minimum value for c. This minimum 
can be defined as: 


c min 


2D/si n<j) 


( 3 ) 


The volume of the unit cell is the product of these three 
lengths. 

The unit cell contains the six yarns depicted in Figure 5. 
The sum of the lengths of all these yarns is equivalent to 
four times the body diagonal length. The volume of yarn 
in the unit cell is equal to the yarn length times the yarn 
cross-sectional area. This relationship can be expressed 
as: 

V y = 7rD 2 \/2o 2 + c 2 (4) 


The fiber volume fraction in the unit cell is equal to the 
yarn volume divided by the unit cell volume times the yarn 
packing factor. The yarn packing factor accounts for 
interstitial space within the yarn bundle. The value used 
for this work is 74%. The value results from 
experimental work performed by C. Pastore at 
NASA Langley. 

MODELING THE EFFECT OF FABRICATION VARIABLES 

Since the yarns are close-packed in this model, the fiber 
volume fraction is solely a function of the braid angle. The 
effect of varying this angle is plotted in Figure 8. Fiber 
volume fraction increases with braid angle. The maximum 
braid angle, 52°, yields a fiber volume fraction of 0.79. The 
maximum angle corresponds with a minimum in the value 
of c. 


A fiber architecture possesses geometric isotropy if the 
variation in effective fiber volume fraction directly 
contributing to a loading direction is constant for any angle. 
The effective fiber volume fraction of a fabric is defined as 
the fraction of fibers aligned in the proper direction in order 
for the applied load to be transferred to the fibers. For 
this model of geometric isotropy the load bearing capacity 
of a fiber in the transverse direction is assumed to be 
zero. This capacity is assumed to be one in the longitudinal 
direction. These assumptions can be applied since the load 
bearing capacity of a yarn in the transverse direction is many 
orders of magnitude lower than that in the longitudinal 
direction. 

The unit cell is composed of six yarns which transverse 
the cell in four distinct yarn orientations. These yarn 
orientations are orthogonal to each other. The yarn 
orientations are translated to correspond with 0, 90, 180, 
and 270 degrees. 

The distribution of fibrous reinforcement is described 
by plotting the effective fiber volume fraction at an arbitrary 
angle in the xy, xz, and yz planes. The braid angle is 
projected onto a given plane. Then the effective fiber 
volume is calculated for any angle in this plane. The fibrous 
reinforcement distribution function in the xy plane is: 

V(S) = 0.25V f sin <J> (lcos(0)1 4 !sin(r?)l) (5) 

The fibrous reinforcement distribution function in the 
xz and yz planes is: 

V(i?) = 0.25Vf cos</> (lcos(^)l 4 lsin(r?)l) (6) 
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Figure 9 plots fibrous reinforcement distributions in the 
xy plane for two different braid angles. The xy fabric 
formation plane corresponds with the through-thickness 
direction in the finished fabric. Minima in the plot occur 
at the principal yarn orientation directions. In these 
directions the other orthogonal yarn pair contributes 
nothing. Maxima occur 45° from the minima. At this 
location all yarns are contributing to load-carrying ability. 
The effective load-carrying ability varies 33% as a function 
of the arbitrary angle. 


90 



Braid angle 

— — 52 degrees 

— 20 degrees 


Figure 9. The distribution of fibrous reinforcement in 
the through-thickness plane 


Figure 10 plots fibrous reinforcement distributions in 
three orthogonal fabric planes for 52° braid angle. 
Euclidean braid fiber architectures with this maximum braid 
angle possess the highest amount of through-thickness load 
carrying ability. As the angle decreases, overall fiber 
volume decreases and the portion of the fiber volume 
contributing to through-thickness load-carrying ability 
decreases. 


90 



0 


Braid angle is 52 degrees 

— _ XY plane 
« ■■■ — XZ — YZ plane 


Figure 10. Fibrous reinforcement distribution in the 
three fabric planes of an Euclidean braid 


SUMMARY 

Euclidean braid fiber architectures do possess a close- 
packed fiber plane. With the close-packed condition, fiber 
volume fraction is solely a function of the braid angle. The 
fiber volume fraction increases with the braid angle. 

The fibrous reinforcement distributions for Euclidean 
braids closest approach isotropic conditions when the load- 
carrying ability in each plane is equivalent. Within a plane 
the fibrous reinforcement distribution is similar to a four- 
leaf clover. Maxima and minima differences correspond 
with a 33% difference in load-carrying ability. Due to the 
degree of isotropy present, no arbitrary angle contributes 
a large percentage of reinforcement. 

In general, geometric isotropy plots are a useful means 
to determine which fiber architecture is most suited for the 
desired loading conditions for a particular application. 
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ABSTRACT 

Explication of the fracture processes of braided composites is 
needed for modeling their behavior. Described is a systematic exploration 
of the relationship between microstructure, loading mode, and micro- 
failure mechanisms in carbon/epoxy braided composites. The study 
involved compression, and fracture toughness tests and optical and 
scanning electron fractography, including dynamic in-situ testing. 
Principal failure mechanisms of tow sliding, buckling, and unstable crack 
growth are correlated to microstructural parameters and loading modes; 
these are used for defining those microstructural conditions which are 
strength limiting. 

INTRODUCTION 

The integrated nature of multidimensional braided composites holds 
promise for overcoming some of the shortcomings of conventional 
laminate composites, while introducing compromises of their own. While 
their general behavioral trends are known there are still many knowledge 
gaps and obstacles to their use. In particular, the high anisotropy and 
interaction of several structural levels have made selection of 
representative volume elements and homogenization procedures difficult. 
At the present, no superior method of modeling and analyzing these 
materials has emerged. 

This paper gives an overview of several projects focused on the 
deformation and failure of braided materials. The approach is 
experimental micromechanics-based, relying on testing of small-scale 
specimens and optical and SEM fractography. There is a prohibitively 
large number of material and braid geometry combinations to be 
exhaustively tested; therefore the focus has been on determining the 
failure sequence and failure controlling microstructural features. The 
goal is the establishment of a heuristic, rule-based description of the 
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failure process/microstructural interaction, and of the definition of a 
mechanical properties unit cell which can be used as a bridge to macro- 
scale behavior. 

MATERIALS AND PROCESSES 

The textile preforms were braided with 12K tows of Celion G30-500 
(BASF Structural Materials, Inc.) graphite fibers. These fibers measure 7 
microns in diameter and have a modulus of 30 Msi. The preforms were 
braided by a 4-step process in a 3 X 14 design where each of the tows is 
interlocked with the other tows to form a true 3-D structure. Reference 
[1] contains a discussion of the 4-step braiding process. 

Composites are fabricated from these textile preforms using a 
vacuum/compression consolidation process. The matrix is comprised of 
two components, Epon 828 (Shell Chemical Company) and Jeffamine T-403 
(Texaco Chemical Company), in a ratio of 100:42 parts-by-weight 
respectively. After mixing, the catalyzed resin is vacuum degassed for 10 
minutes at < 1 Torr. The preform is placed in a mold and submerged in the 
liquid resin and the evacuation step is repeated. The mold is closed and 
the composite is press cured for 3 hours at 100 C. A mechanical stop is 
used to control the thickness of the composite. Typical dimensions for 
the composite samples are 0.100 X 0.750 X 10.0 inches. More details of 
the fabrication aspects are given in [2], The fiber volume fraction was 
typically 55%. 

Composite tubes were also used in these experiments. They were 
fabricated from a textile preform braided in much the same manner as the 
flat braid but the machine bed was circular instead of rectangular. The 
tube preforms were impregnated and consolidated by resin transfer 
molding (RTM) with Tactic 123 (Dow Chemical Co.) and catalyzed with 
Millamine 5260 (Milliken Chemical Co.) in a ratio of 100:17 parts-by- 
weight respectively. This matrix system requires a cure of 1 hour at 70°C 
followed by 2 hours at 170°C. 

STRUCTURAL DESCRIPTION 

A structural description of the composite is necessary for any 
consideration of failure modes. The preform structure has been defined 
from a textile point of view [1], Such a description, however, is only 
partly descriptive of the final component. While computational unit cells 
have been proposed for braids and weaves [3,4], care must be taken in 
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extracting them from the actual structure. A different approach will be 
taken here. 

Figures la and 1b show the textile preform and the subsequent 
composite coupon fabricated from it. In this case, a white tracer tow has 
been added indicating the location of a single tow. Figures 2a and 2b show 
the top and side views of a single tow extracted from the composite. 
Based on these, a schematic of the tow location, along with identification 
of the important microstructural descriptors and the surface/volume 
mapping can be constructed (Figure 3). 

Note that a fiber tow lies on a diagonal-type plane and sinuously 
spans the thickness of the specimen. This fiber does not exactly match the 
apparent surface pattern. The true cycle length, or the return of the tow to 
its starting position may be quite long and depends on the braiding cycle. 
For the purpose here, it is sufficient to define one cycle as the return to 
the equivalent point on the diagonal plane. Figure 4 is a further 
idealization and contains only two representative crossing tows. The 
diagonal length of this cell is one half the cycle length previously defined. 

Figure 4 can be used for creating an efficient representation of the 
tow structure. Consider a unit cell as commonly defined in x-ray 
crystallography [5]. This type of unit cell allows the use of Miller Indice 
notation, which describes directions, planar orientations, and directional 
relationships in compact form. The X, Y, and Z axes become [100], [010], 
and [001] respectively. The diagonal plane on which the tows lie are of 
the form (110), and the crossing points or intersection of the tows are 
generally on the (lll)-type planes. The cell here is, of the the 
orthorhombic type, and actual cell parameters are established by 
measurement from the fabricated composite. 

TEST RESULTS 

Table 1 contains selected test results from small coupons, along 
with the external load conditions and specimen details. The matrix (neat 

resin) specimen failed along the expected 45° shear lines. The axial 
specimens initially failed in shear along the apparent braid angle. SEM 
fractographs are shown in Figure 5. Internally, the failure was of a mixed 
mode on the (110) and (lll)-type planes. Note the rapid drop in axial 
strength with the increase of braid angle. The overall trend with fiber volume 
fraction is shown in Figure 6. In later stages of the failure sequence fiber 
buckling occurred. For transverse compressive loading, the failure was 
tows sliding over each other on the (lll)-type planes. 
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The notched specimens were based on ASTM E399 compact fracture 
toughness specimens [6]. This test was developed for metals but attempts 
have been made to use it for laminate composites. For braided materials 
this test appears useful only as a worst-case indicator. If the notch is 
along the braid direction, the specimen tends to fail in cleavage at values 
below that of the matrix. The crack front moves in general along the 
(HO)-type planes. For notches perpendicular to the braid direction the 

failure load is considerably higher but the crack length, area and failure 
mode are difficult to determine. Figure 6 is an SEM of this case. 
Considerable additional work is required to remove uncertainties from the 
testing results. 

Tests were also conducted on components. Axial compression of one 
inch diameter braided tubes and T-sections followed the same trends as 
the coupons. Figures 7 and 8 are micrographs of these specimens. 

The loads and stresses given in Table 1 are those for yielding. There 
is much post-yield phenomena in these materials, with higher values 

possible with continued testing, along with increasing damage. At some 
point the damage becomes catastrophic and the specimen fails. 

DISCUSSION 

The previous results represent only selective, illustrative data. Our 
experience with these materials runs to hundreds of samples and a variety 
of tests, and these will be amplified in other publications. The common 
trends and generalizations that have emerged will be outlined here. 

The dominant descriptive parameter is the braid angle. The 

differentiation between the apparent and internal braid angles is 

approximately 3° for most braids, and the surface angle can be used as an 
approximation in most cases. 

In the initial stages of the failure sequence, the tows act as units, 

giant pseudo-fibers. It is useful to partition the matrix into inter- and 

intra- tow components. The inter-tow matrix and fibers form the tows 
which are surrounded by the intra-tow matrix. If the total fiber fraction 
is V tota |, and the fiber fraction in the tow is V^ ow , then the inter-tow 

matrix fraction is 


' ^tota/ ^tow 


( 1 ) 
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This assumes the tows are straight. 

Most tows are not circular; combined with tow curvature this has 
the effect that in later stages of failure the individual fibers in a tow do 
not carry the load evenly and failure occurs by sequential tearing. This is 
a strength lowering but energy absorbing mechanism. 

Fiber-matrix interface debonding is frequently the initiating event 
in failure, allowing either tow sliding or the formation of cracks at resin- 
rich areas. It appears that the accumulation and growth of these cracks is 
encouraged by the resin-rich areas. 

There is also a stochastic aspect introduced by structural 
variability. Most braided materials contain variation in the braiding 
angle, and this is linked in a complicated way to the local density and 
fiber volume. The variation in even high quality braids can be as much as 
2 degrees. Considering the effect of braid angle on compressive strength 
given previously, there is bound to be an uncertainty in strength values of 
components. Further, in components, cracks in tension have a probability 
of being in the proper orientation for runaway crack growth. 

Do these observations point the way to any solutions of the 
shortcomings? Stronger interfaces are one possibility. To reduce the 
effects of matrix-rich areas fewer fibers per tow or stress dispersants 
such as microballoons could be used. But the toughness/energy absorption 
vs. strength tradeoffs are not fully known, and much work remains in 
understanding these materials. 

As a closure, the question of a mechanical properties unit cell will be 
considered. Based on the foregoing experimental work, it appears that 
this could be derived from the unit cell defined previously. This unit cell 
would reflect the fracture-strength anisotropy, and would contain scaled 
equivalent volumes of matrix and fiber. These would be corrected to 
reflect Equation 1, and the interface angles would be based on observed 
failure processes.. Variability can be introduced by considering a 
distribution of cells. We have used such an approach by assigning failure 
criteria and using Weibull statistics, and have been able to model some 
macro failure processes. This work is continuing. 
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Table 1 Test Results 


Braid 



Angle 

Axial Yield Stre 

3 

110,000 psi 

13 

66,000 

1 8 

44,000 

23 

16,000 

Matrix 

13,000 


Fiber Volume - 55% approx. 
Thickness - 0.100 approx. 

Angle 

Transverse Stress 

21 

5,000 psi 

Angle 

Notched Strength 

21 

140 pounds 

Angle 

Notched Strength 

21 

52 pounds 

Matrix 

60 pounds 
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Figure 1 a - Braided Textile Preform, b. - Molded Composite. Note the 
tracer tow. The composite is 3/4 inch wide. 
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b. 


Figure 2. a - Side View of Single Tow. Note the “diamond” where the 
tow intersects the surface. The lines indicate the apparent 
braid angle. 

b - Top View of Same Tow. Lines indicate the width of the 
specimen 
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Figure 3. Surface/Volume Mapping of Braid. 










Figure 6. Notched E399-Type Specimen, Notched Perpendicular to 
Principal Braiding Direction. 



Figure 7 Axially Compressed Tube 
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Introduction 

Inhomogeneous finite elements are an attractive alternative to homogeneous elements in 
the mechanical analysis of fabric reinforced composites (Figure 1). These elements greatly 
simplify the mesh generation problem created by the complex reinforcing geometry. 
However, this advantage has some drawbacks associated with it. Convergence, with 
diminishing element size, becomes less certain. Also, the computation of stresses within the 
various constituent materials of an element becomes a problem. This paper addresses both 
these concerns. 

The convergence can be improved by replacing the inhomogeneous elements with 
special homogeneous elements whose properties are chosen to match the inhomogeneous 
element response to simple average strain states. One and two dimensional examples are 
considered. The three dimensional application is discussed without examples. The analysis 
also provides the basis for an approximate solution for the average stresses within each 
constituent material of each inhomogeneous element. This makes an approximate 
stress/failure analysis possible. 



PLAIN WEAVE 



2X2 BRAID WITH STUFFER 

FIGURE 1. FABRIC REINFORCED COMPOSITE MICROGEOMETRIES 
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The Analysis 

The analysis uses the unit cell concept to create a boundary value problem that fully 
characterizes the reinforcing microgeometry. Division into subcells then establishes a fmite 
element model of the unit cell and reduces the reinforcing complexity in each subcell to the 
point where average, subcell, constituent stress levels are meaningful. For example, 
consider the plain weave unit cell in Figure 2. Subcells that were adequate for stiffness 
analysis must be reduced in size to yield meaningful detailed stress information. 

The individual subcell stiffness matrix, [k], can be obtained by numerical integration 
of the general finite element energy formula [k] = fff B T DB dvol where the matrix D 
contains only material property distribution functions and the matrix B contains only 
displacement mode shape derivatives (Ref. 1). Once a library of different subcell stiffness 
matrices has been created it remains to transform them into the global coordinates and 
assemble them into an overall stiffness matrix for the unit cell. The surface nodal forces 
and the average strains in each subcell, corresponding to each of the six independent strain 
states of 3-D elasticity, may then be solved. The surface nodal forces give the average 
stresses on the unit cell surfaces. The stiffness coefficients of the composite may then be 
computed. The convergence of this process is considered first. The subsequent 
computation of the constituent stresses is taken up later. 




SUBCELL FOR 
STIFFNESS ANALYSIS 


0 


SUBCELL FOR 


STRESS ANALYSIS 


FIGURE 2. DIVISIONS OF UNIT CELL 
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One Dimensional Example 


The tension bar made from dissimilar materials (Figure 3) illustrates the convergence 
problems associated with inhomogenous elements and simple displacement mode shapes. If 
a finite element node coincides with the point of material discontinuity then the elements 
become homogeneous, the strain in each element becomes uniform, and the analysis 
converges abruptly to the true displacements. But, if the material discontinuity is always 
contained within some element, as in Figure 3.., then the solution is approximate and the 
accuracy and convergence rate depend on the choice of assumed displacement mode 
shapes. For example, consider a linearly varying displacement mode shape within each 
element and an internal node placement at the 1/3 and 2/3 points along the bar. Each 
subsequent refinement divides each element into three equal segments. The middle element 
remains inhomogeneous as element size decreases. Figure 3 is a plot of bar elongation 
error as the element size diminishes. The load point displacement approaches the exact 
solution monotonically. Although rapid, this convergence is less satisfactory than that of 
the homogeneous element solution. 
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FIGURE 3. TENSION BAR PROBLEM 


Improved Convergence 


The convergence rate for the previous problem can be improved by using higher 
order displacement mode shapes (Ref. 2). It can also be improved by making a 
modification that overrides the source of error. The error source is the inability of the 
energy formulation, in combination with the linear displacement assumption, to distinguish 
between dissimilar material stiffnesses in series and parallel (Figure 4). The use of low order 
displacement modes is a presumption of parallel response when in reality the stiffnesses 
are functioning in series, in this particular application at least. The use of higher order 
modes would permit the analysis to make such a distinction. The same refinement can be 
attained by intervention of the analyst. Since the two materials in the center element are 

truly arranged in series, an effective modulus, E, of this element can be computed from the 

elementary series formulation 2 _ 1_+ 1 where subscripts L, R designate left and right 

E E l Er 

halves of the bar. The inhomogeneous center element can then be replaced by a 
homogeneous one with the modulus E. This correction leads to the immediate convergence 
of the deflection analysis. The logic in this argument seems trivial but as the dimension 
increases to two and three it becomes more complicated. 



ELEMENTS IN PARALLEL 



ELEMENTS IN SERIES 


FIGURE 4. ELEMENTS IN SERIES AND PARALLEL 
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2-D Example 


Figure 5 contains an example of a bimetallic composite in the form of bonded sheets 
of dissimilar materials. In the natural coordinates of the material the composite stiffnesses 
and internal stresses can be established from elementary mechanics. The results can be 
transformed into any other coordinates. In the coordinate system of Figure 5 A the material 
can also be analyzed using inhomogeneous finite elements. A unit cell and a subcell 
division are shown in Figure 5B. Using the displacement modes usually associated with 8 
node- isoparametric-brick elements ( Ref. 1), generalized plane strain analysis, and the 32 
node finite elements grid shown; the analysis overestimates the x and y moduli by over 
20%. Refinement of the finite element grid leads to the composite moduli estimates of 
Figure 5C. 

The convergence may be improved by recognizing how the reinforcing layers provide 
stiffness and then applying the appropriate rule of mixtures formulae while substituting 
homogeneous orthotropic elements in place of the inhomogeneous ones. Stiffnesses in 
parallel apply in the two principal reinforcing directions. Stiffnesses in series apply normal 
to the material boundaries. Figure 5C shows that this leads to abrupt convergence of the 
finite element sequence. The Appendix contains the relevant equations for this analysis. 
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FIGURE 5. LAMINATED 2-D COMPOSITE MATERIAL 



2-D Generalization 


Consider generalizing the 2-D case to include more than two materials while 
restricting the material distribution within the subcell. Assume the material boundaries can 
be represented by a number of straight lines radiating from a common point within the 
subcell (Figure 6A). If additional lines are drawn from that point to each subcell comer 
node then each material will be contained in one or more homogeneous triangular Finite 
elements. Comer displacements of the rectangular subcell are Fixed by the applied strain 
case. Edge conditions must be approximated. The stresses within each material of the 
subcell can now be computed corresponding to each independent unit strain case. The 
nodal forces along the sides of the subcell are also available for calculating the mean 
stress/strain relations for the subcell. The subcell can then be considered homogeneous and 
anisotropic in the unit cell analysis. 

If the common point lies outside the subcell boundaries then a similar analysis can 
be performed using some trapezoidal elements (Figure 6B). In either case full displacement 
continuity is not preserved between adjacent rectangular subcells. 



FIGURE 6. INHOMOGENEOUS/HOMOGENEOUS ELEMENT REPLACEMENT 



Edge Boundary Conditions 


The aforementioned boundary conditions on nodes created by the intersection of 
material boundaries and subcell edges are key points (Figure 7A). The assumptions 
used here are motivated by consideration of the special case where a material boundary 
and subcell edge cross at right angles (Figure 7B) and the surrounding material/loading 
arrangement is symmetric across the same subcell edge. Then the response will be 
symmetric. This implies zero tangential displacement along the material boundary 
(except for rigid motion) and zero nodal force normal to the material boundary. 

The unit shear strain case (Figure 7C), relative to the same axes, represents 
antisymmetric loading with zero displacement normal to the material boundary (except for 
rigid motion) and zero nodal force component parallel to the material boundary. 

For material boundaries which intersect subcell edges at angles other than 90° the 
true boundary conditions cannot be rationalized so easily. Nevertheless, it is assumed that 
the same displacement and force conditions remain as suitable approximations to the true 
conditions of equilibrium and compatibility across adjacent subcell boundaries. 
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FIGURE 7. INTERFACE BOUNDARY CONDITIONS 
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3-D Analysis 


To extend the analysis into three dimensions it is assumed that within any subcell of 
the rectangular 3-D array of subcells representing the unit cell, the dissimilar material 
boundaries consist of planes which intersect at a single straight line within the subcell 
(Figure 8). •Call this line the "common material axis". Consider the smallest sphere 
containing the subcell. A major circumference of the sphere exists which defines a plane 
which is normal to the common material axis and passes through the subcell centriod. By 
assuming the stresses do not vary along the common material axis the analysis reduces to a 
series of 2-D problems in the circumferential plane. Each 2-D problem resolves the 
displacements as a result of an applied unit strain with reference to the circumferential 
plane. The stresses within each wedge of dissimilar material may then be solved and the 
average stresses over any planar area within the subcell computed. These average stresses 
over the subcell surfaces permit the calculation of the equivalent anisotropic constants of 
the subcell. The boundary conditions at node points formed by the intersection of 
dissimilar material planes and subcell boundaries are established in the same way as in the 
2-D case. Thus, inhomogeneous subcells are replaced by homogeneous ones in the unit cell 
analysis. 



CIRCUMFERENTIAL 

PLANE 

FIGURE 8. 3-D INHOMOGENEOUS SUBCELL GEOMETRY 


178 


i ii i iimniumiui . i uiiiii i 1 1 muiuwMMijiiiii aim i i,, iiiiii ju, n nijiiniHiumiumtaMflM&M n iiiiiiiiiiiiiiiiuiniiiiaiii, iuiiii.ii nil ifia iik, mm 


Global-Local Analysis 


It is likely that critical regions of high stress within a unit cell can be identified 
beforehand, making a complete stress analysis unnecessary. Certain subcells or portions of 
subcells will be of interest and can be isolated for detailed study through global-local 
analysis. This reduces the computational effort. For example, a crude model of a unit cell, 
based on a few subcells, can be used to obtain composite stiffnesses and define the 
displacements on the boundaries of a smaller volume of microstructure that contains a 
single subcell of further interest. A more refined grid may then be superposed on the 
smaller volume and more detailed stresses or displacements obtained for either a failure 
analysis or a more refined analysis of some smaller portion of that volume (Figure 9). Such 
a sequence of grid refinement may be expected to yield detailed finite element average 
stresses and strains in regions of tow bypass, tow contract, and sharp tow curvature. The 
problem of resolving these average element strains into average constituent strains and 
stresses will be considered next. 
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FIGURE 9. TYPICAL SEQUENCE OF INCREASING ANALYSIS DETAIL 


179 




Constituent Stresses 


Returning to the one dimensional example, the determination of the stresses in the 
truly homogeneous elements of the tension bar is the same as in any finite element analysis. 
The strain is first computed from the displacements of the element nodes. The 
homogeneous, uniaxial stress/strain law of the material then leads to a stress calculation. 
For an inhomogeneous element the average strains are computed in the same way. 
Recourse to the stiffness in the series model that was used to compute the equivalent stiffness, 
E, also provides the basis for the constituent material stress calculation. The equations of 
strain compatibility and stress equilibrium plus the individual constituent material 
stress/strain relations are adequate in number for the stress calculations (Ref. 2). 

In the general 2-D and 3-D cases the same mini-finite-element models that form a 
basis for equivalent stiffness calculations also provide a mechanism for backing out the 
constituent stresses (Figure 10). First, the element comer displacements, in conjunction 
with the displacement mode shapes establish internal displacements. Derivatives of the 
mode shapes establish detailed strains. The volume averages of these detailed strains 
establish the subcell average strains. Each average strain component, along with the mini- 
finite element model yields average strains in each constituent material. Constituent 
stress/strain laws then yield average stresses at the constituent material level. 

Select a rectangular unit cell 

SUPERPOSE RECTANGULAR FINITE ELEMENT GRID ON UNIT CELL 

REPLACE INHOMOGENEOUS ELEMENTS WITH MIMICKING HOMOGENEOUS ONES 
(VIA MINI FINITE ELEMENT MODEL) 

ANALYZE UNIT CELL FOR EACH OF THE SIX INDEPENDENT UNIT STRAIN CASES 

COMPUTE COMPOSITE MODULI 

COMPUTE UNIT CELL STRAINS FOR APPLIED LOAD CASES OF INTEREST 

RATIO CONSTITUENT STRESSES IN EACH ELEMENT FROM UNIT STRAIN CASES TO 
APPLIED LOAD CASES 

COMBINE CONSTITUENT STRESSES FOR THE APPLIED LOAD CASE STRAINS 
FIGURE 10. STEPS IN STRESS ANALYSIS 
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Concluding Remarks 


There is a need to do routine stress/failure analysis of fabric reinforced composite 
microstructures to provide additional confidence in critical applications and guide materials 
development. Conventional methods of 3-D stress analysis are time consuming to set up. 
run and interpret. A need exists for simpler methods of modeling these structures ana 
analyzing the models. The principal difficulty is the discrete element mesh generation 
problem. Inhomogeneous finite elements are worth investigating for application to these 
problems because they eliminate the mesh generation problem. However, there are 
penalties associated with these elements. Their convergence rates can be slow compared to 
homogeneous elements. Also, there is no accepted method for obtaining detailed stresses 
in the constituent materials of each element. This paper shows that the convergence rate 
can be significantly improved by a simple device which substitutes homogeneous elements 
for the inhomogeneous ones. The device is shown to work well in simple one and two 
dimensional problems. However, demonstration of the application to more complex two 
and three dimensional problems remains to be done. Work is also progressing toward more 
realistic fabric microstructural geometries. 

FABRIC MICROSTRUCTURAL ANALYSIS NEEDED 

CONVENTIONAL (HOMOGENEOUS) FINITE ELEMENTS DIFFICULT TO APPLY 
TO FABRIC 

SIMPLER METHODS/MODELS NEEDED 

INHOMOGENEOUS ELEMENTS ATTRACTIVE ALTERNATIVE 
SLOW CONVERGENCE PROBLEM CAN BE REMEDIED 
CONSTITUENT MATERIAL STRESSES OBTAINABLE 
ANALYSIS NOT PROVEN YET 
WORK REMAINS TO BE DONE 

FIGURE 11. CONCLUSIONS 
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Appendix 


The rectangular inhomogeneous element of Figure 12 consists of two homogeneous isotropic 
materials (A and B). Parallel to the material boundary plane (y = o) the average Young’s moduli (E,, Hj) 
and Poisson’s Ratios ( Ujj, u^) can be approximated by the parallel Rule of Mixtures: 

E, = Ej = E A v A + E B v B ? u ^ = Uyj= t) A v A +u B v B 

where v designates material volume fraction and subscripts A, B designate the material. The shear 
modulus (Gjj) can be approximated by the same rule 

g « =g a v a + g b v b 

Normal to the material boundary plane the average Young’s modulus (E-) and shear modulus 
(G^) are given by the series Rule of Mixtures: 

_l_=\^ + v B , 1_ = v A + v B _ 

Ey E a E b Gjj G a G b 

Considering y = o to be a plane of isotropy, the Poisson’s Ratio (ujj or u s ) in that plane can be 
approximated by (Ref. 3): 

V A 

where R = 

B 

These equations are sufficient to support the construction of a homogeneous finite element 
stiffness matrix after the moduli are transformed into the global x,y,z coordinates. 

x 
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Figure 12: 2-D Analysis for Two Materials 
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